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Abstract 



It is believed that a novel state of matter - Quark Gluon Plasma (QGP) will be 
transiently produced if normal hadronic matter is subjected to sufficiently high tem- 
perature and/or density. We have investigated the possibility of QGP formation in 
the ultra-rclativistic collisions of heavy ions through the electromagnetic probes - pho- 
tons and dileptons. The formulation of the real and virtual photon production rate 
from strongly interacting matter is studied in the framework of Thermal Field Theory. 
Since signals from the QGP will pick up large backgrounds from hadronic matter we 
have performed a detailed study of the changes in the hadronic properties induced by 
temperature within the ambit of the Quantum Hadrodynamic model, gauged linear 
and non-linear sigma models, hidden local symmetry approach and QCD sum rule ap- 
proach. The possibility of observing the direct thermal photons and lepton pairs from 
quark gluon plasma has been contrasted with that from hot hadronic matter with and 
without medium effects for various mass variation scenarios. The effects of medium 
induced modifications have also been incorporated in the evolution dynamics through 
the equation of state. We find that the in-medium effects on the hadronic properties in 
the framework of the Quantum Hadrodynamic model, Brown-Rho scaling and Nambu 
scaling scenarios are conspicuously visible through the low invariant mass distribution 
of dileptons and transverse momentum spectra of photons. We have compared our 
evaluation of the photon and dilepton spectra with experimental data obtained by the 
WA80, WA98 and CERES Collaborations in the heavy ion experiments performed at 
the CERN SPS. Predictions of electromagnetic spectra for RHIC energies have also 
been made. 
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Chapter 1 
Introduction 



As per contemporary wisdom, hadrons i. e. baryons and mesons are made up of quarks 
which interact by the exchange of gluons. This interaction is governed by a non-abehan 
gauge theory called Quantum Chromodynamics (QCD). Gluons, the gauge particles 
of QCD unlike photons of Quantum Electrodynamics, carry colour charge and arc self 
interacting. This endows the QCD interaction with remarkable properties as a function 
of the relative separation, or equivalently, the exchanged momenta. At short distances, 
or large momenta (q), the effective coupling constant Oisiq^) decreases logarithmically. 
This means that the quarks and gluons appear to be weakly coupled at very short 
distances, a behaviour referred to as asymptotic freedom [1, 2]. At large separations, the 
effective coupling, it is believed, progressively turns stronger resulting in the phenomena 
called quark confinement [3, 4, 5] which describe the observation that quarks do not 
occur isolated in nature, but only in colour singlet hadronic bound states as mesons 
(quark- antiquark bound states) and baryons (three quark bound states). In this case 
one also has chiral symmetry breaking [6] which expresses the fact that the quarks 
confined in hadrons do not appear as nearly massless constituents but are endowed 
with a dynamically generated mass of several hundred MeV. In nature there exist six 
quark flavours of three colours in addition to eight bi-coloured gluons. 

The QCD renormalization group calculation predicts that as the temperature and/or 
density of strongly interacting hadronic matter is increased, the interactions among 
quanta occur effectively at very short distances and is thus governed by weak cou- 
phng due to asymptotic freedom, while the long range interactions become dynamically 
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screened due to Debye screening of colour charge [7, 8, 9]. The interaction thus decreases 
both at small distances as well as for distances large in comparison to typical hadronic 
dimensions. As a consequence, nuclear matter at very high temperature/density ex- 
hibits neither confinement nor chiral symmetry breaking. This new phase of QCD 
where the bulk properties of strongly interacting matter arc governed by the funda- 
mental degrees of freedom - the quarks and gluons, in a finite volume is known as the 
state of quark gluon plasma (QGP). 

There are sufficient reasons to expect that the transition between the low and the 
high temperature manifestations of QCD is not smooth but exhibits a discontinuity, 
indicating the occurrence of a phase transition at some intermediate value of tempera- 
ture [10, 11]. In this context one refers to the possibihty of two different kinds of phase 
transitions, namely, the chiral symmetry restoring transition and the deconfinement 
transition. Now, it is necessary to realize that though QCD is firmly established as the 
theory of strong interactions, the long range behaviour of QCD is not well understood. 
This is because, at larger separations the strong coupling becomes large and invalidates 
perturbative approaches. As a result the study of long range phenomena such as phase 
transitions and hadronic properties and interactions using QCD is highly inhibited. One 
of the attempts to study such phenomena non-perturbatively on a discrete space-time 
lattice is the lattice gauge theory [12]. 

At low energies, the QCD vacuum is characterized by non-vanishing expectation 
values of certain composite operators, called vacuum condensates e.g. (qq) which de- 
scribe the non-perturbative physical properties of the QCD vacuum. They also act as 
the order parameter for the chiral phase transition. QCD calculations on the lattice 
suggest a first order (discontinuous) chiral transition for two massless quark flavours 
and a second order (smooth) transition for three flavours. The chiral phase transition is 
indicated by a strong drop in the value of the condensate (qq) at the transition tempera- 
ture, Tc believed to have values between 130-160 MeV. Lattice calculations also indicate 
that there is a sharp change in the energy and entropy densities at Tc [13]. This can 
be interpreted as an increase in the number of degrees of freedom, corresponding to a 
deconflnement transition. Thus restoration of chiral symmetry and deconflnement seem 
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Figure 1.1: Phase diagram of strongly interacting matter. 

to occur concomitantly. 

The quark-hadron phase transition is also expected to occur at high baryon density 
even at zero temperature. There is however a large uncertainty in the critical baryon 
density of transition. Model calculations predict critical densities to lie anywhere be- 
tween 4 to 10 times normal nuclear matter density. One expects a smooth connection 
between the high temperature and high baryon density phase transitions, giving rise 
to a continuous phase boundary Tc{pb)- Fig. (1.1) shows a typical phase diagram of 
strongly interacting matter. For T < Tc{pb), the effective description is in terms of 
hadronic degrees of freedom, whereas for T > Tc{pb) the degrees of freedom carry the 
quantum numbers of quarks and gluons. However, recently very interesting theoretical 
developments regarding the QCD phase diagram has taken place. At high densities, 
quarks may form Cooper pairs and a new colour superconducting phase may exist. The 
Pb — T phase diagram may have a critical or tricritical point somewhere along the phase 
transition line [14, 15, 16]. 

According to the standard big bang model, the universe has gone through the QCD 
phase transition and matter is supposed to have existed as QGP during the microsecond 
epoch after the big bang when the temperature was about 100 MeV. This has important 
consequences in cosmology [17]. The possible remnants that may have survived that 
primordial epoch till date can provide valuable clues about the nature of the phase 
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transition. A first order cosmological QCD phase transition scenario [18] could lead 
to the formation of quark nuggets made of u, d and s quarks at a density somewhat 
larger than normal nuclear matter density. Primordial quark nuggets with sufficiently 
large baryon number could survive even today and could be a possible candidate for the 
baryonic component of cosmological dark matter [19]. The QCD phase transition has 
its importance in astrophysics because cold deconfined matter is also believed to exist 
in the cores of neutron stars where the baryon density is supposed to be high enough for 
the occurrence of the phase transition. However, these cosmological sources of quark 
matter can only provide indirect and limited information. As a result, attempts to 
create conditions conducive for the production of QGP (energy density ~ 1 GeV/fm'^) 
in the laboratory through ultra-relativistic heavy ion collisions (URHICs) [20, 21, 22] 
have been a major sphere of activity of physicists over the last several years. A number 
of accelerators have been designed for this purpose. Experimental data from Pb-Pb 
collisions at 158 GeV/nucleon at the Super Proton Synchrotron (SPS) at CERN are in 
the last phases of analysis. The Relativistic Heavy Ion Collider (RHIC) at Brookhaven 
is about to start functioning and experimental data will be available in the very near 
future. Also, the Large Hadron CoUider (LHC) at CERN is expected to be ready in 
about five years from now. Heavy ion {e.g. ^°^Pb and ^^^Au) coUisions at the RHIC 
and LHC with centre of mass energies 200 AGeV and 5500 AGeV respectively are 
expected to produce extremely high energy densities. The question is whether one can 
produce a large enough and sufficiently long-lived composite system so that collective 
and statistical phenomena can occur. An even bigger challenge is to extract information 
on the dynamics and the properties of the earliest and the hottest stage of such collisions, 
because even if QGP is produced it would only have a very transient existence. Due 
to colour confinement quarks and gluons cannot escape from the collision and must 
combine to colour-neutral hadrons before travelling to the detectors. Hence, all signals 
emerging from the QCD plasma will receive substantial contribution from the hadronic 
phase. Therefore, a detailed study of the hadronic phase is important in order to 
disseminate the emissions from the QGP. 

It is necessary to understand that the main difference between the QCD phase 
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transition in the laboratory and in the early universe is that the effects of gravitation 
is very important for the later. The dynamics of the QCD phase transition in the early 
universe is therefore governed by the Einstein's equations in the Robertson- Walker 
space-time. The solution of Einstein's equation with equation of states for strongly 
interacting matter results in a characteristic time scale which is of the order of few 
micro seconds unlike URHlCs where time scales of the order of a few fm/c are involved. 
Consequently, in the early universe the evolution takes place in a very leisurely pace 
compared to the interaction time scales for the quarks and gluons which is ~ few fm/c. 
Therefore the QCD phase transition in the early universe occurs in an environment of 
complete thermal equilibrium which might not be satisfactorily realized in URHlCs. 

In the following Sections we will briefly discuss the the general features of the evo- 
lution of relativistic heavy ion collisions and the proposed signals of the quark-hadron 
phase transition. 

1.1 Formation and Evolution of QGP in URHlCs 

As the two heavy ions collide at very high energies, they deposit a substantial part 
of their kinetic energy into a small region of space. Depending on the energy density 
achieved, the initial state of the system will be either in the form of a QGP or a 
hot/dense hadronic gas. The evolution of an URHIC with an intermediate state of 
QGP can be described in terms of four regimes as discussed below. Since these vary 
widely in nature the physics governing these stages are quite different. The evolution 
appears as shown in Fig. (1.2) when projected in the plane of the longitudinal coordinate 
z and time t. 

a) Formation of QGP: Immediately after the collision of the two Lorentz con- 
tracted nuclei, the scattered partons decohere and free stream in the longitudinal di- 
rection. Rescatterings of these partons lead to locally thermalized QGP after a time 
generally taken to be ~1 fm/c. The value of the thermalization (realization of kinetic 
equilibrium) time is still rather uncertain. There are mainly two approaches to the 
microscopic description of the process by which the fully coherent parton wavefunc- 
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Figure 1.2: Space-time diagram of the collision in the t — z plane. 

tions of the two nuclei in their ground state just before the collision, evolve into a 
locally thermal distribution of partons in a QGP - the QCD string breaking and the 
partonic cascade. In the string breaking model colour strings formed between mutually 
separating partons after the collision fragment producing quark anti-quark pairs which 
eventually bring about thermalization. Some well known realizations of this approach 
are VENUS [23], FRITIOF [24] and RQMD [25]. The parton cascade model [26] on 
the other hand is based on the parton picture and renormalization group improved 
pcrturbative QCD. Whereas the string picture runs into conceptual difficulties at very 
high energy, the parton cascade becomes invalid at lower energies where most partonic 
scatterings are too soft to be described by perturbative QCD. 

b) Expansion of the QGP: Assuming that the interactions of quarks and gluons 
are sufficiently small at the temperatures achieved in URHICs, the energy density, 
pressure etc. can be calculated in QCD using thermal perturbation theory. Driven 
by the high internal pressure, the thermalized QGP expands according to the laws of 
relativistic hydrodynamics [27]. The hydrodynamic equations for an ultra-relativistic 
plasma admit a boost invariant solution describing a longitudinally expanding fireball 
with constant rapidity density [28, 29]. Relaxing the boost invariance by assuming 
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a Gaussian-like multiplicity density can lead to interesting features [30]. The most 
important question that arises during this part of the evolution is that of chemical 
equilibration of the partons. It is generally believed that gluons, because of their larger 
colour degeneracy equilibrate chemically much faster than the quarks. We have found 
that even light quark flavours fail to achieve chemical equilibrium during the lifetime 
of the plasma [31, 32]. 

c) Hadronization and the mixed phase: Expansion of the QGP proceeds till 
the critical temperature Tc is reached. At this instant, the phase transition to hadronic 
matter starts. Through the process of hadronization the coloured particles - quarks 
and gluons combine to form colour-neutral hadrons. The order of the transition is still 
somewhat controversial. Mostly, for the purpose of calculation of the signals of QGP 
it is assumed that the phase transition is of first order. The released latent heat main- 
tains the temperature of the system at Tc even though the system continues to expand. 
During this time a coexistence of quark and hadronic matter follows by a Maxwell con- 
struction. This mixed phase persists until all the matter has converted to the hadronic 
phase. Though not generally considered in a simple picture, one realizes that the large 
latent heat associated with a strong first order phase transition may lead to deflagra- 
tion waves during the expansion of the fireball along with possible superheating and 
supercooling. A detailed microscopic description of the hadronization process is yet to 
be achieved. 

d) The hadronic phase and freeze-out: The resulting hadronic matter expands 
and cools as long as the system can sustain interactions. As we will see later in this 
thesis, the properties of hadrons at high temperature are modified non-trivially due to 
interactions thereby changing the equation of state and consequently the cooling pro- 
cess. Once the mean free path of the hadrons becomes comparable to the dimensions of 
the system they decouple and free stream towards the detector. The temperature when 
this occurs is the freeze-out temperature Tf. The value of the freeze-out temperature 
is still an unsettled issue. The contribution to the signal output depends substantially 
on the freeze-out temperature, even more so when transverse expansion of the fireball 
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is considered. 

1.2 Signals of QGP 

The size of the plasma volume is expected to be at most a few fermis in diameter. It may 
live for a duration ~ 1-10 fm/c out of an overall freeze-out time of about 100 fm/c at the 
most. One realizes that once the system is produced, its space-time evolution cannot be 
controlled. In fact, the only experimentally controllable initial parameters are the mass 
numbers of the colliding nuclei and the collision energy. In addition, a handle on the 
impact parameter in each collision can be obtained by forming event classes of different 
multiplicities and transverse energies with a correlation to the energies observed in the 
zero-degree calorimeter. With these few controllable initial parameters, information 
of the whole space-time evolution of the system must be extracted from the various 
observables measured in the final state. Regardless of whether or not QGP is produced 
in the initial stages, the system turns into a system of hadrons. Hence, the particles that 
are detected are mostly hadrons along with photons and leptons. These hadrons, mostly 
light mesons like pions, kaons etc. make up the large multiplicity in relativistic heavy 
ion collisions. In fact, out of about 2500 particles created in central Pb-Pb collisions at 
the CERN SPS more than 99% have turned out to be pions. The hadrons are emitted 
predominantly from the freeze-out surface whereas photons and leptons are produced at 
all stages of the evolution as indicated in Fig. (1.2). Being strongly interacting particles, 
hadrons can provide only indirect evidence having undergone significant reinteractions 
between the early collision stages and their final observation. Leptons and photons in 
contrast are weakly interacting and are considered to be direct probes. 

The above discussions indicate the complexities involved in the identification and 
investigation of the QGP. However, various signals have been proposed which we will 
briefly discuss in the following. 

1.2.1 Probes of the Equation of State 

Since hadronic matter and the QGP is separated by a phase transition, most likely of 
first order one looks for modifications in the dependence of energy density e, pressure 
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P, and entropy density s of the evolving matter on the temperature T and/or the 
baryon chemical potential /ib- One searches for a rapid rise in the effective number of 
degrees of freedom, as expressed by the ratios e/T^ or s/T^ over a small temperature 
range. However,the thermodynamic quantities T, s and e are not directly measured in 
experiments. Usually, these arc identified with the average transverse momentum (pt), 
the rapidity distribution of the hadron multiplicity dN/dy and the transverse energy 
dEx/dy respectively. When (px) is plotted as a function of dN/dy or dEx/dy, one 
expects first a rise corresponding to the increase in the number of degrees of freedom in 
the hadronic phase, then a saturation during the persistence of the mixed phase followed 
by a second rise when the change from colour-singlet hadrons to coloured partonic 
objects is completed [33] . Another important feature which could indicate the collective 
nature of the evolution is the observation of transverse flow effects in the momentum 
spectra of final state particles. This is due to the fact that if the hfetime of the produced 
collective system is long enough, a strong collective flow is generated [27] and the heavier 
the particle is the more transverse momentum it gains from the transverse flow. Again, 
non-central collisions, i.e. collisions with non-zero impact parameter give rise to a 
different kind of flow - the asymmetric flow. In this case transverse flow is generated 
by the pressure gradients in the transverse plane leading to an azimuthally asymmetric 
flow [34] . This in turn causes the azimuthal angle distributions of flnal state hadrons to 
be asymmetric. Identical particle interferometry also provides an independent source 
of information regarding the space-time dynamics of the collision [35]. By studying 
measured two-particle {e.g. tttt, KK or NN) correlation functions in different directions 
of phase space, it is possible to estimate the transverse and longitudinal size, the lifetime 
and flow patterns of the hadronic flreball at the moment of freeze-out. The transverse 
sizes measured in URHlCs are found to be larger than the radii of the incident particles 
clearly indicating that the produced hadrons have rescattered before emission [36] . 

1.2.2 Signatures of Chiral Symmetry Restoration 

Strangeness enhancement and increase in antibaryon production relative to p-p or p-A 
collisions are some of the proposed signatures of chiral symmetry restoration. The basic 
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argument in both cases is the lowering of the threshold for production of strange hadrons 
and baryon-antibaryon pairs [37] . An optimal signal is obtained by considering strange 
antibaryons which combine both effects. The enhancement of strangeness in URHICs 
in simple terms is connected to the fact that in an ideal baryon-dense quark matter the 
production of ss is enhanced compared to that of light quark flavours u and d because 
the Fermi energy of the already present light quarks is higher that the strange quark 
mass. In Pb-Pb collisions at the CERN SPS, an enhancement in strangeness production 
by about a factor 2 has been observed from the {K + K) / {'k) ratio [38]. Also, a clear 
specific enhancement in the yield of Q." + W' per negative hadrons by a factor ~ 10 
have been observed by WA97 [39] relative to p-p and p-Be collisions. It is believed 
that as multistrange hadrons are difficult to produce due to high mass-thresholds, the 
strangeness increase could have an origin at the partonic level before hadronization. It 
is believed that domains of disoriented chiral condensate (DCC) could provide a more 
direct signal for the restoration of chiral symmetry in URHICs [40] . These correspond 
to isospin singlet, coherent excitations of the pion field which would decay into neutral 
and charged pions in such a way that there is a significant possibility of observing 
a large surplus of charged pions over neutral pions in certain regions of phase space. 
Restoration of spontaneously broken chiral symmetry is also reflected in the thermal 
modiflcation of the hadronic spectral function [41, 42, 43, 44, 45] particularly through 
the mass shift of the vector mesons in hot/dense medium. These modifications can be 
studied by analyzing photon, dilepton as well as hadronic spectra. 

1.2.3 Probes of Colour Deconfinement 

The suppression of J/^ production is considered as a direct probe of the colour decon- 
finement phase transition. The J/^ is a bound state of a cc pair dominantly produced 
by the fusion of gluons. In a deconfined environment like a QGP, the binding of a cc 
into a J/^ is suppressed due to the fact that the screening length is less than the bound 
state radius [46]. On the other hand, the J/^ may also be suppressed in a hadronic 
scenario due to nuclear absorption. This description is a probabilistic one where one 
assumes that the probability that a produced J/'^ escapes without making collisions is 
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oc exp(— L/A) where A = l/cXabs'^'O, is the nuclear density and Uabs is the cross section 
which controls the probability that the J/^ gets destroyed in a nuclear collision. The 
quantity L is the average length travelled by the J/^ in nuclear matter, is related to 
the transverse energy of the collision. This picture explains both p-A data as well as 
A-A data upto the S-U system. However, the Pb-Pb data from the CERN SPS has 
created a great deal of excitement [47]. This is due to the discontinuity observed in 
the ratio of J/^ to the Drell-Yan cross section at L ~ 8 f m (which corresponds to 
Et ~ 50 MeV and impact parameter 6 ~ 8 fm). Nuclear absorption apparently can 
not explain this discontinuity and one needs to invoke partonic degrees of freedom and 
colour confinement. It should be mentioned that there have been other attempts to 
explain this striking feature in a hadronic description. It has been argued that even if 
the J/^ escapes the nuclei, it can be destroyed at a later stage of the collision due to 
scattering on other produced particles, commonly referred to as comovers. A trustable 
estimate of such contributions is still lacking. 

Another possible way of probing the colour structure of the produced matter is by 
studying the energy loss of a fast parton, also known as jet quenching [48, 49]. The 
parton loses its energy either by excitation of the penetrated medium or by radiation. 
The magnitude of the energy loss is proportional to the strong couphng constant a^. It 
has been observed that the energy loss of a parton jet is greater in A-A coUisions than 
in p-p or p-A. Theoretical estimates have inferred that the energy loss of a parton in 
hadronic matter {—dE/dx ~1 GeV/fm) is much more than in QGP {—dE/dx ~0.1-0.2 
GeV/fm). In this respect, since jets in hadronic matter is suppressed more in hadronic 
matter than in QGP, jet "unquenching" is a signal of deconfinement. 

1.2.4 Electromagnetic Probes 

These include real (7) and virtual (7*) photons. Virtual photons decay producing pairs 
of leptons. Photons and dileptons are considered the cleanest signals of quark gluon 
plasma [50]. Because of the very nature of their interaction, they decouple immediately 
and leave the system without any distortion of their energy-momentum carrying with 
them the information from within the reaction zone. Hence, electromagnetic probes are 
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the only direct probes of QGP. They are emitted at all stages of the evolution but do 
not get masked by the details of the evolution process. Their production cross-section 
is strongly dependent on temperature and hence are copiously produced from the hot 
phase of the evolution. However, the photons and dileptons from the thermal phase of 
the evolution, referred to as thermal photons and dileptons have to compete with a large 
background [51, 52, 53, 54] due to production from many other mechanisms and this 
complicates the process of extraction of information about the basic thermal system 
we intend to study. The general features of the photon and dilepton spectra will be 
discussed in the following Section. 

1.3 Real and Virtual Photons: General Features 

Let us briefly discuss the various sources of production of photons and dileptons during 
various stages of an URHIC and the specific domains of phase space where they are 
known to dominate. Let us first consider real photons. In the low transverse momentum 
region the contribution from pseudoscalar meson decays involving 7r° and rj clearly dom- 
inate photon production. In fact they account for nearly 95% of the total photon yield 
in a collision. Photons emitted due to primary interactions among the partons of the 
colliding nuclei form the principal background in the large transverse momentum region. 
These are called prompt or QCD photons. The decay photons can be isolated experi- 
mentally by invariant mass analysis and the prompt ones can be accurately estimated 
by perturbative QCD calculations. These contributions are then subtracted out to get 
the thermal photon spectra. Thermal photons from the QGP phase arise mainly due 
to the QCD Compton {qg — > 57) and annihilation {qq — > 517) processes. The emission 
rate resulting from these reactions turn out to be infra-red divergent and can be evalu- 
ated [55, 56, 57] in the framework of Hard Thermal Loop (HTL) [58, 59] resummation 
in QCD. The hadronic matter is mainly composed of the pseudoscalar-isovector pion 
(tt) and the spin-isospin vector rho (p) mesons. One also includes the vector-isoscalar 
omega {cu), the pseudoscalar-isoscalar eta (77) and the axialvector-isovector oi mesons. 
Photons from the hadronic matter are emitted from reactions of the type hh — > hj 
(where h is one of the hadrons tt, p and 77) as well as from the decay of short-lived 
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>10 fm/c ~ 1 fm/c < 0.1 fm/c 

Figure 1.3: Schematic diagram showing the different contributions to the total photon 
yield in a ultra-relativistic heavy ion collision. 

hadrons [55] . In a phase transition scenario these are emitted during the later stages of 
the coUision when the system has cooled to temperatures below Tc. Hence the photon 
spectra due to the thermal hadronic phase is expected to have a steeper slope compared 
to the QGP. In Fig. (1.3) we have schematically shown the different contributions to the 
overall photon yield in a relativistic collision of heavy ions as a function of the trans- 
verse momentum of the emitted photons. One observes that the higher the transverse 
momentum of the photon the earlier they are produced in the collision process. 

The different contributions to the production of lepton pairs as a function of their 
invariant mass are shown schematically in Fig. (1.4). The principal source of thermal 
dileptons from the QGP and hadronic phases are the quark- antiquark and pion annihila- 
tion processes respectively. In the high mass region these compete with Drell-Yan pairs 
which are produced in primary interactions between incoming partons in the very early 
stages of the collision. The J/^ peak marks the cut-off scale for thermal pairs from the 
plasma. Around this region the decays of D and B mesons also become an important 
source. The vector mesons which decay both during the expansion and after freeze-out 
can be identified easily from their characteristic peaks in the spectrum. Below the vec- 
tor mesons, Dalitz decays of 7r°, rj, rjf and cu mesons provide the dominant source for 
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1 2 3 4 M [GeV] 
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Figure 1.4: Schematic diagram showing the different contributions to the total dilepton 
yield in a ultra- relativistic heavy ion collision [60]. 

dilepton production. The fact that higher mass lepton pairs are produced earlier in the 
collision process is evident from the time axis in Fig. (1.4). We emphasize that whereas 
the dileptons from hadronic matter have distinct features like the p, uo and peaks, 
the photon spectra is completely structureless. Hence it is very difficult to disentangle 
the thermal photons emitted from hadronic matter from those which have their origin 
in quark matter. In this case an accurate estimation of thermal photons originating 
from hadronic matter is of utmost importance in order to comment on the formation of 
QGP in relativistic heavy ion collisions. Such an estimation must incorporate medium 
modifications of hadronic properties in the evaluation of emission rates as well as in the 
equation of state (EOS) of the interacting hadronic matter. The change in the mass and 
decay width of a vector meson propagating in a medium occurs due to its interaction 
with the real and virtual excitations in the medium. Medium induced modifications of 
the properties of vector mesons, for example, the p and uj are likely to show up clearly 
in the invariant mass spectra of dileptons through the shifting and/or broadening of 
the respective peaks. In the following Section we will give a brief introduction to the 
study of hadronic properties in a thermal medium. 
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1.4 Medium Effects on Hadronic Properties 

In URHICs hadronic matter is expected to be formed after a phase transition from QGP 
when the plasma has cooled beyond the phase transition temperature Tc. Even if such 
a phase transition does not occur, realization of hadronic matter at high temperature 
(~ 150-200 MeV) and/or baryon density (a few times normal nuclear matter density) 
is inevitable. As a result the study of hadronic interactions at high temperature and 
density assumes great significance. However, progress in our understanding of hot 
and dense hadronic matter has been retarded since the underlying theory of strong 
interaction, QCD, is nonperturbative in the low energy regime. This severe constraint 
has lead to considerable amount of work on model building [61, 62, 63, 64] in order to 
study the low energy hadronic states. 

The principal contention of this thesis is to study the medium effects on thermal 
photon and dilepton spectra. We will be mainly concerned with the vector mesons 
p and CO. Various investigations have addressed the issue of temperature and density 
dependence of hadronic spectra within different models over the past several years. Ac- 
cording to Brown and Rho [65] the requirement of chiral symmetry (in particular the 
QCD trace anomaly) yields an approximate scaling relation between various effective 
hadronic masses, which implies that all hadronic masses decrease equally with tem- 
perature. The reduction in p meson mass has also been observed in the gauged linear 
sigma model [66] at low temperature; however, near the chiral transition point it shows 
an upward trend. The nonlinear sigma model claimed to be the closest low energy 
description of QCD shows the opposite trend, i.e. the effective p mass increases with 
temperature [67]. A similar qualitative behaviour of the p mass has been observed in 
the hidden local symmetry approach [68]. The relation between the self energy and 
the forward scattering amplitude has also been utilized to study the change of hadronic 
properties in the medium [69, 70], where the effects of non-zero temperature is rather 
small. 

In the Quantum Hadrodynamic Model (QHD) of Walecka [71, 72] scalar and vector 
condensates generated by the nucleon sources are themselves responsible for the modi- 
fication of the nucleon mass. The vector meson mass gets shifted due to the decrease of 
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the nucleoli mass which appears through thermal loops in the vector meson self energy 
[73], the imaginary part of which characterizes the response of the nuclear system to 
external (electromagnetic) probes. 

In-medium QCD sum rules are useful to make constraints on the hadronic spectral 
functions at finite temperature and density [74]. In the QCD sum rule approach the 
hadronic properties are related to the scalar and tensor condensates of quark and gluon 
fields. Due to lack of understanding of the behaviour of these condensates near the 
critical point, the hadronic properties at finite temperature in this approach is not 
firmly established. The spectral function of vector mesons can be parametrized in 
vacuum from the experimental data obtained in e+e~ hadrons (or from hadronic 
decays of the r lepton) for various isovector and isoscalar channels [75, 76]. The in- 
medium spectral function of the vector meson is then obtained by modifying the pole 
and the continuum structure as a result of its interaction with the constituents of the 
thermal bath. The vector meson masses and the continuum thresholds are taken to vary 
with temperature according to Brown-Rho (BR) and Nambu scaling [77] scenarios. 

So we see that a wide range of variation of hadronic properties with temperature 
are predicted by the models cited above. Other models e.g., those proposed by Rapp et 
al [78] and by Klingl et al [79], where the effects of non-zero baryon density (baryonic 
chemical potential) is dominant over non-zero temperature will not be discussed. This is 
because we intend to study the hot baryon free (central rapidity) region of URHICs. We 
shall consider various scenarios for the shift in the hadronic spectral function at finite 
temperature and evaluate its effects on the experimentally measurable quantities, the 
photon and dilepton spectra originating from a thermalized system formed in URHICs. 

1.5 Organization of the Thesis 

The thesis is organized as follows. In Chapter 2 we shall discuss the formalism of photon 
and dilepton production from a medium of interacting particles at finite temperature 
and/or density starting from first principles using perturbation theory. The rates of 
photon production from partonic interactions in the QGP as well as from hard primary 
interactions of partons will be considered. Chapter 3 is devoted to the study of spectral 
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modifications of vector mesons in tlie medium and tlie evaluation of static (fixed tem- 
perature) rates of photon and dilepton production from hadronic matter within various 
models. The medium masses of the vector mesons are calculated from the pole positions 
of the effective propagators which are obtained in terms of the self energies calculated 
in the framework of thermal field theory using a few well-known models. We have con- 
sidered the QHD model, the gauged linear sigma model, the gauged non-linear sigma 
model, and the hidden local symmetry and QCD sum rule approaches. The space-time 
evolution of these static rates using relativistic hydrodynamics is discussed in Chap- 
ter 4. The medium effects in the space-time evolution enters through the equation of 
state which is manifested in the cooling rate as well as in the estimation of the initial 
temperature of the produced matter. The final photon multiplicity with and without 
a QGP in the initial state is obtained. These are compared with data obtained by the 
WA80, WA98 and CERES experiments. Chapter 5 contains the thesis summary and 
related discussions. We have listed the invariant amplitudes of all the photon producing 
hadronic reactions and decays that we have used in the Appendix. 



Chapter 2 

Formulation of Electromagnetic 
Emission Rates 

Electromagnetically interacting particles - photons and dileptons are excellent probes of 
the thermodynamic state of evolving strongly interacting matter likely to be produced in 
ultra-relativistic nucleus-nucleus collisions. This is because electromagnetic interactions 
are strong enough to lead to a detectable signal and yet are weak enough to let the 
emitted photons and leptons escape from the finite nuclear system without further 
interactions. Hence the spectra of photons and dileptons can provide information of 
the properties of the constituents from which they were emitted. 

For most purposes the emission rates of photons and dileptons can be calculated 
in a classical framework. It was shown by Feinberg [80] that the emission rates can 
be related to the electromagnetic current correlation function in a thermalized system 
in a quantum picture and, more importantly, in a nonperturbative manner. Generally, 
the production rate of a particle which interacts weakly with the constituents of the 
thermal bath (the constituents may interact strongly among themselves, the explicit 
form of their coupling strength is not important) can always be expressed in terms of 
the discontinuities or imaginary parts of the self energies of that particle [81, 82]. In 
this Chapter we will make a detailed study of the connection between the emission 
rates of real and virtual photons and the spectral function of the photon which is con- 
nected with the discontinuities in self energies in a thermal system [83]. This in turn 
is connected to the electromagnetic current correlation function [50] through Maxwell 
equations. Alternatively, using the kinetic theory approach the photon emission rates 
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can be written in terms of the equilibration rate of photons in a thermal bath. The 
imaginary part of the photon self energy tensor in this case is related to the exclusive 
kinetic rates of emission and absorption in the system. We will discuss these facets 
regarding the formulation of the photon and dilepton emission rates from a thermal 
medium extensively in Section 2.2 after a brief review of propagators in thermal field 
theory in Section 2.1. The invariant rates of emission of both hard and soft thermal 
photons and dileptons from QGP will be dealt with in Section 2.3 where we will com- 
pare the static (fixed temperature) rates due to different processes contributing to the 
thermal yield. 

2.1 Review of Thermal Propagators 

As is well known, propagators play a central role in the description of the dynamics of 
systems of particles using quantum field theory. In this Section we will briefiy discuss 
the in-medium (thermal) propagators in Thermal Field Theory [82, 84, 85, 86] which 
will be used extensively in the thesis. We will begin by first defining the propagators 
in vacuum. 

The free propagator of a complex scalar field (j) propagating with a momentum p in 
vacuum is defined as 



iA%p) = J d'xe'P-^{T{<t>{x)<l>*{0)})o 



= ^y— - ;e^0+ (2.1) 

The operator T appearing within angular brackets ensures that the field operators are 
time-ordered and the subscript '0' indicates that there are no interactions. The vacuum 
propagator for fermions is defined as 

iGlpip) = J d'xe'^-^{T{M^)Mmo 

_ i(|^ + m)a(3 



p2 _ _|_ 



(2.2) 



where a and /3 denote the spinor indices of the fermion field ■0. In a similar way the 
free propagator of a vector field is defined as 

iD%ip) = I d'xe'P-^{T{A,ix)AM})o (2.3) 
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Depending on the nature of the quanta of the vector field A'^ the vector propagator can 
take the following forms. In the case of a vector particle of mass m we get 



= ''-f (2.4) 



which describes the free propagation of the uj meson, for example. For charged vector 
particles e.g. the p meson, the fields also carry isospin indices and we have 

^ p'' — + ie 

where i and j denote components in isospin space. In the case of massless vector fields 
corresponding to photons or gluons the propagator will contain a gauge parameter as a 
reminder of the arbitrariness of the gauge-fixing condition. The photon propagator in 
vacuum is obtained as 

Some well known choices are, ^ = 1 which is the Feynman (or Lorentz) gauge and ^ = 0, 
the Landau gauge. For gluons we need to add a Kronecker delta is colour space so that 
the free gluon propagator is 

,pOa6. ^ -t5-'[g^-^ + {^-l)p^p'^/p^] 

The gauge parameter ^ must be absent from any physical quantity we calculate. 

It must be noted that the fields 0, ip and ^4^ appearing above are free fields and the 
expectation values are calculated between noninteracting vacuum states. The propaga- 
tors defined through Eqs. (2.1), (2.2) and (2.3) are referred to as Feynman propagators. 

In the presence of interactions these propagators have to be redefined with inter- 
acting Heisenberg fields in place of the free fields and interacting vacua instead of the 
free vacua. The interacting propagator can be expressed in terms of the bare (non- 
interacting) propagator using perturbation theory. In the scalar case, for example, the 
exact propagator in the presence of interactions. A, is obtained as 



A = Ao AqIIA 



(2.8) 
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Where, IT is the self energy of the particle due to interactions. This equation is known 
as the Dyson- Schwinger equation for propagators. 

Let us now study the situation in a medium at finite temperature (and density). 
We will be interested in a system in thermal equilibrium. Hence we will assume that 
the interaction slowly switches off as wc go into the remote past and the fields become 
noninteracting fields satisfying the free equations of motion. These fields appear in the 
definition of the free propagators in the medium. The thermal propagator has more 
structure than the vacuum case as a result of different combinations of time-ordering on 
the real time contour [82, 84, 102]. In the real time formalism there are four non-trivial 
propagator structures possible which are collected in a 2x2 matrix [87, 88]. For scalars 
the free thermal propagator is defined as 



,Af{p) zAf(p) 

/ d^xe^P-''(T{^{x)q^*{0)})^ J d^xe^P-''{{(f)*(O)0{x)})^ 



(2.9) 



where the operator T denotes anti-time-ordered product. The subscript 'T' indicates 
that a thermal average is being performed. 

In order to obtain the thermal propagators in momentum space one follows the 
usual procedure of expanding the field operators in terms of the creation and annihila- 
tion operators and making use of the commutation relations between them. The four 
components are then obtained as 

Kip) = \^—--2T:i5{p''-m^)r]{p-u) 

p^ — + le 

Al^{p) = -2'Ki5{p^ - w?)[q{p ■ u) + e{-p ■ u)] 
Af ip) = -2mSip^ - m^) [r]{p ■ u) + Oip ■ u)] 
Afip) = , ~\ -2Tri5ip^ -m'')r)ip-u) 



. . . ... . (2-10) 

p^ — — le 

where ri{p ■ u) = 9{p ■ u)fBE{z) 9{-p ■ u)fBE{-z). fsE = [e'' - 1]"^ is the Bose 
distribution with z — {p -u — /i) /T, is the four velocity of the thermal bath and ^ is 
the chemical potential. We observe that the elements of the matrix propagator Aq are 
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not independent. From their definitions one can sec that Aq^ and Aq^ can be expressed 
in terms of Aq^ and Aq^. Also, the Kubo-Martin-Schwinger [89] periodicity condition 
yields Al\p)^e-^Al\p). 

Similarly, for fermions the corresponding thermal propagators are. 



0(12), 



J d^xe'P--{T{^a{x)MmT -I d'xe'P-^{MO)M^)})T 
Explicitly, the four components are 

1 



(2.11) 



GViP) ^ 



+ 27ri5{p — m )rj{p ■ u) 



Ga^^\p) = 27ri(;^ + m)a/3S{p^ - m'^)[ri{p ■ u) - e{p 



p2 _ ^2 _|_ 

27ri(jf^ + m)ap5{p'^ — m'^)[r]{p ■ u) — 6{—p ■ u)] 

2 _ ,^2^ 

-1 



:/3 



U) 

2 



p2 _ ^2 _ 



+ 27riS{p — m )ri{p ■ u) 



(2.12) 



where, r]{p • u) = e{p ■ u)fFD{z) + 9{-p ■ u)fFD{-z), fro = [e^ + 1] S the Fermi- 
Dirac distribution. For the fermions the KMS anti-periodicity condition leads to — 

Lastly, we define the finite temperature propagators for vector particles: 



/ rf4xe'P-(T{A^(a;)A,(0)})0 / d^x e'P-^{{A,{Q)A^{x)})% 
J d'xe^P-^{{A,{x)A,{0)})^T J d'xe'P-^f{A^{x)A,{0)})^j. 



.(2.13) 



As before, one will have additional indices corresponding to colour or electric charge; 
henceforth we will not mention them explicitly. Apart from the Lorentz indices, the 
explicit forms of the thermal propagators will be similar to the scalar case. For neutral 
particles, the chemical potential /i will be absent in the definition of the phase space 
factor riip ■ u). 
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It is important to note that the real time propagators as given by Eqs. (2.10) and 
(2.12) consist of two parts - one corresponding to the vacuum, describing the exchange 
of virtual particles and the other, the temperature dependent part, describing the par- 
ticipation of real (on-shell) particles present in the thermal bath in the emission and 
absorption processes. The tempcratm^c dependent part docs not change the ultra-violet 
behaviour of the theory as it contains on-shell contributions and has a natural cut-off 
due to the Boltzmann factor. Therefore, the zero temperature counter term is adequate 
for the renormalization of the theory. However, the infra-red problem becomes more 
severe at finite temperature [82, 90]. 

Thermal field theory in the real time approach can be reformulated by diagonalizing 
the 2x2 matrix propagators described above. A well-known possibility is to diago- 
nalize to a matrix constructed from the Feynman propagators [87]. The free thermal 
propagator defined by Eq. (2.9) can be written as 



Ao = U 



Ao 








U 



(2.14) 



where 



r] + 9{p- u) 



The exact propagators in the medium can be defined analogously as Eqs. (2.9), (2.11) 
and (2.13) with interacting Heisenberg fields instead of the free fields. In this case we 
write 

A 

-A* 



U 



u 



(2.15) 



where A is the matrix of interacting thermal propagators. Using thermal perturbation 
theory A can be expressed in terms of Aq. One obtains 



A = Ao AqIIA 



(2.16) 
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where 11 now is the self-energy matrix; 



n = 


nil 


ni2 




n 







nsi 


n22 







-n* 



Matching the elements appearing in the diagonal of Eq. (2.16) we have 

1 



(2.17) 



A = Ao + AonA 
Prom Eq. (2.17) it also follows that 



p2 _ _ n _|_ 



(2.18) 



Ren = ReHii 



(2.19) 



and 



Imll = e(po) tanh(po/2T) Imllii for bosons 

Imll = e(po) coth((po + A*)/2T) Imllii for fermions 



(2.20) 



In the following we will discuss the vector (spin 1) propagator in some detail. It 
is very similar to the scalar case except now one has to take into account the Lorentz 

structure of the propagator and the self energy. The exact propagator (matrix) D^j, 
can be diagonalized as above and the diagonal element satisfies Dyson equation 



D ^ D + D W^D 



(2.21) 



which gives 



D-l = {Dl^'' - n^., (2.22) 

where D^^^ is the vacuum propagator for vector particles. The quantity —ili^i, is the 
sum of all one particle irreducible (IPI) self energy insertions. In has a vacuum and a 
medium part so that 

(2.23) 



w = w + n 

vac ' 



1X1/ 

med' 



where 



vac 



(9 



p2 



)nvac(p'), 



(2.24) 



is the vacuum contribution to the self energy. 
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Naively, it would appear that finite temperature corrections to quantum field theory 
breaks Lorentz covariance since the rest frame of the heat bath selects out a specific 
frame of reference. However, by defining a fluid four-velocity with temperature 

— * 

defined in the fluid rest frame where it^ = (1, 0), a manifestly covariant formulation can 

be achieved [84, 91]. Using the techniques of tensor decomposition it can be shown that 
for a vector particle propagating with four-momentum = {po,p), 



(2.25) 



where nT,med and Hi mcd are Lorentz invariant self-energy functions which characterize 
the transverse and longitudinal modes. A^^'^ and B^'^ are the transverse and longitudinal 
projection tensors given by 

1 



and 



p2 -p2 



(P - Po)i9''" - uV) - pf^p" - p'ouV + poiu'^p" + pV) 



1 



which satisfy the following algebra: 



pIp'^p" + p'^u^u" - pop'^iu^p" + pV] 



, (2.26) 



(2.27) 



y -'-'til' 

y^M^ _|_ ^M'^ 






2 
1 



Pf'p" 

p2 



The two functions nr,med and n^^med are obtained by contraction: 

^2 



n 



L.med 



p 



\P\ 



riT.med — 



med 



n 



L.med 



)■ 



(2.28) 



(2.29) 
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Now, for massive vector particles, 



(2.30) 



Using Eqs. (2.22-2.28) the effective propagator becomes 



B 



fjLV 



p2 _ jyip. _|_ Ylj, 



where 



n 



T(L) 



n. 



T(L),med 



(2.31) 



(2.32) 



and m, we recall, is the bare mass of the particle. The real part of U.t(l) affects the 
dispersion relation of the particle in the medium. The displaced pole position of the 
effective propagator in the rest frame of the propagating particle {i.e. where the three 
momentum of the particle is zero) gives the effective mass of the particle in the medium. 
The imaginary part of IIt{l) is connected to the decay width. 

A different scheme in the formulation of finite temperature field theory known as the 
'R/A' formalism [92], is to diagonalize to a matrix composed of retarded and advanced 
propagators which are known to have better analyticity properties than the Feynman 
ones. In this case the analogue of Eq. (2.14) is 



V 









W. 



(2.33) 



The free retarded and advanced propagators are defined as 

I d^xe'P-^e{xo)mx),m])o 



p2 _ ^2 _|_ j^^p^ 



rf4^e*P-^(-a;o)([0(a;),0(O)])o 
i 



p2 _ ^2 _ ^gp^ 



(2.34) 



The matrices V and W depend on the momentum as well as the thermal factor con- 
taining the distribution functions. Their exact forms are given in Ref. [92]. For the 
case of massive vector particles one arrives at the following equation for the effective 
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retarded propagator at finite temperature: 



D 




(2.35) 



p2 _ _|_ JjR p2 _ _|_ JjR rn^p- 



i2' 



where 11^ and 11^ are respectively the retarded transverse and longitudinal components 
of the self energy. For photons we have [57] 



Before we end our discussion of finite temperature propagators let us briefly men- 
tion about the imaginary time formalism or Matsubara formalism which has been used 
extensively in the literature. In the imaginary time formalism, the form of the propa- 
gator at finite temperature is the same as that in vacuum but the time component of 
the four-momentum takes discrete values, i.e. po — 2mriT{— {2n+ l)Tr iT) for bosons 
(fermions) with n = — oo to -|-oo, the vertices are the same as the zero temperature the- 
ory and the loop integral / d^p/{27i)'^ is replaced by the sum iTJ^nl d^p/ {27r)^. There 
are standard methods to evaluate the sum over the frequencies [82]. The propagators 
in the imaginary time formalism can also be obtained by proper analytic continuation 
of the real time propagators [93, 94] . Another method, known as the SACLAY method 
has also been used extensively in the literature [58, 95]. This method uses the mixed 
representation of the propagator i. e. it depends on the three-momentum and Euclidean 
time. 



We note at the onset that the nature of emission of real and virtual photons depends 
crucially on the size of the hot thermal system from which they are emitted. If the 
system is large enough the photons will rescatter and thermalize and their momentum 
space distribution will be given by the Planck distribution. The corresponding emission 
rate will then be that of black body radiation which depends only on the temperature 
and the area of the emitting body but not on its microscopic properties. Since the 
typical size of systems produced in heavy ion collisions is much less than the mean free 
path of the photons, they are likely to escape the hot zone without rescattering and the 




(2.36) 



2.2 Thermal Emission Rates 
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emission rate in this case depends on the dynamics of the thermal constituents through 
the imaginary part of the photon self energy. We will begin by demonstrating this in 
a kinetic theory framework [60, 81] before going on to a more rigorous scheme. We 
will in general denote the four-momenta of the real and virtual photons by and 
respectively. 

2.2.1 Emission Rate as Equilibration Rate 

We know that the probabihty of a photon of 4-momentum p'* = [E, p) to be absorbed 
in matter is given by 

UE)-i^ E f dn^^uf}S\p+T.k^-T.f'f)\Mip+{^}^{f})\' n Ml±nf), 

{iUf} {iUf} 

(2.37) 

and the probability of emission is given by 

re(^)-^ E [ dn{iy,{fy5'{^kf-p-^h)\M{{i}^p+{f})\' n nf{l±n,). 
»,{/} «,{/} 

(2.38) 

Here, is the equilibrium distribution function, {i} and {/} denote the initial and 
final state particles, and denotes the phase space integration including spin 

and polarization sums. Now, 

where the squared matrix elements of the forward and backward processes have been 
taken to be equal on account of time revcrsality. If f{E, t) is the momentum distribution 
function of photons, the rate of decrease due to absorption is f{E, t)Ta{E) and the rate 
of increase due to emission is {1 -\- f{E,t))Te{E). So the time evolution equation for / 
is 

^ = -/r„ + (i + /)re. (2.40) 

The solution of this equation is 

f{E,t)= fBE{E)+c{E)e-^' (2.41) 

where, F = Fa — Fg is the equilibration rate, Jbe is the Bose distribution function and 
c{E) is a function which depends on the initial conditions. Since the system under 
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consideration is small enough for the photons to escape immediately on production, we 
must have f — 0- Consequently, Eq. (2.40) reduces to 



dl 

dt 



= Te = JbeT. (2.42) 



Now, the equilibration rate is related to the imaginary part of the photon self energy 
through Imn = EV{E) [96]. For a real photon = and we have from Eq. (2.29), 
Imn = ImHr = Inin^/2. Using 

dN=-^J{x,p)d^xd^p 

for a real photon we get from Eq. (2.42) 

The quantity hsIV^ contains information about the constituents of the thermal bath 
and thus is of great relevance. After this relatively simple but illuminating exercise, we 
will proceed to arrive at this result from more general considerations in the following. 

2.2.2 Emission Rate and Photon Spectral Function 

We begin our discussion with the dilepton production rate. Following Weldon [83] let 
us define A^^ as the exact Heisenberg photon field which is the source of the leptonic 

current J^. To lowest order in the electromagnetic coupling, the scattering matrix 
element, Shi, for the transition | /) —>-\H; is given by 

Shi = -ie {H; l+r \ J d^x J^x) A^{x) \ I) , (2.44) 

where | /) is the initial state corresponding to the two incoming nuclei, | H; is the 
final state which corresponds to a lepton pair plus the rest of the interacting system. 
The parameter e is the renormalized charge which couples the leptonic current with the 
virtual photon field A^. Since we assume that the lepton pair does not interact with 
the emitting system, the matrix element can be factorized as 

{H-m- I Jl{x)A^{x) I /) = {H I A^{x) I /)(/+/- I Ji(x) I 0), (2.45) 
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where | 0) is the vacuum state. Putting = ilj{x)'j^ip{x) where ip is the field operator 
for the leptons, one obtains the expectation value in terms of the Dirac spinors u{pi) 
and v(p2) as 

( i+r I jUx) I 0) - ^^^i )^^^(P^ ) ei(j>i+P2ya> (2.46) 

where denote the Dirac matrices, Ei — y^pf + m^, with i — 1,2 are the energies of 
the leptons and V is the volume of the system. Therefore, 

^ _^^^iPMP2) U^^^n-x^H I A^{x) I I). (2.47) 

The leptons of four-momenta pi and p2 are produced from a single virtual photon of 
four- momentum ? = (?o, ?) so that qQ — Ei-\- E2 and q — pi+P2- 

Assuming that a thermalized system is produced in the collision, the dilepton mul- 
tiplicity N is obtained by summing over the final states and averaging over the initial 
states with a weight factor Z{P)~^ e~^^^; 

^-zp)V¥' ' WW ^ 

where Ej is the total energy in the initial state, Z{P) is the partition function and 
(3 — l/T is the inverse temperature. After some algebra N can be written in a compact 
form as follows: 



(2.49) 



■(27r)3Ei {2nyE2' 
where Lf^^, is the leptonic tensor defined by 

^'^^ = I E U{p,)r V{P2)V{P2)Y U{pi) 

spins 

= p'iP2+P'2Pi-^9'\ (2.50) 
and H/j^i, is the photon tensor 

Using translational invariance we can write 

(/ I A^{y) \H) = {I\ e'^-yA^{0)e-'^-y \ H) = e'^'-^^^^ {I \ A^(0) | H) , (2.52) 
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where V is the four-momentum operator and pi and ph are the total four-momenta of 
the initial and final states respectively. We now have 



(2.53) 

Using the conservation of four-momentum, pi — pn + q and the completeness relation 

^|/)(/|=lwcgct 
I 

H,,^VLe-^^'DlM, (2-54) 

where O (= V.t) is the four-volume of the system and D^^ is the component i-D^J, of the 
exact photon propagator D^j, defined in the last Section. The time ordered propagator 
is the (1, 1) component of D^,^. In coordinate space it is defined as 

iD'^A^) ^ -^^T.{H\T{A,{x)AM}\H)e-^''- 

= e{xo)D>^{x) + 9{-xo)D<,{x). (2.55) 

where is the component iDj^, xq is the time component of the four vector, X/^ — 
{xo,x) and 9{xo) is the step function. Using the integral representation of the 9- 
functions, 

dz e-'y^ 



9{y) = I J 



277 z + ie 

and taking the Fourier transform we get [97] 



,^11 / ^ duo 



— uj + ie (f — Lij — ie 
Using the Kubo Martin Schwinger (KMS) relation in momentum space. 



(2.56) 



^^.(?o,g) = e^^°L><(?o,?), (2.57) 

we have 

^> (50, q) = - ^ ^l_^J mDll{qo, q). (2.58) 
The rate of dilepton production per unit volume (N/fl) is then obtained as 
dN _ 2e^ ruu.r,ii(^ d^pi d^p^ 
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Now, the spectral function of the (virtual) photon in the thermal bath is defined as 

pAlo^q) = / d'xe"^-^Y.{H \ [A,{x),A,m \ H)e-^^-, (2.60) 



27rZ{f3) 
so that, we have [82, 97] 



Dlli^o, q) = r du P^'^^^^^l - 2z7:fBE{qo)p,u{qo, q), (2.61) 

J-oo Qn — (jj -\- le 



-oo Qo— co + te 
where fsEiqa) = [e'^^" — 1]~^- This leads to 



ImDlliqo, q) = -7r[l + 2fBE{qo)]pAqo, q)- (2-62) 



In terms of the photon spectral function the dilepton emission rate is obtained as 

fi - 2..^L^V„.(*. ^/..(*). (2.63) 

This relation which expresses the dilepton emission rate in terms of the spectral function 
of the photon in the medium is an important result. Inserting, 



j dSS^^\pi+P2-q) 



and using the identity 



2m^, / 4m2 



the dilepton rate {dR = dN/d x) can be expressed as 



where m is the lepton mass and a denotes the fine structure constant. We will now 
proceed to evaluate the photon spectral function. 

As is well known, it is not the time-ordered propagator that has the required analytic 
properties in a heat bath, but rather the retarded one. We thus introduce the retarded 
propagator which will enable us to express the dilepton rate in terms of the retarded 
photon self energy. The retarded photon propagator in momentum space is defined as 

zD;^,{qo,q) = / d'x e'^-^e{xo)Y. (H \ [A,{x),A^{0)] \ H)e-^^'', (2.66) 



H 
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which leads to the relations 



lmDll = {l + 2fBE)lmDj^^ (2.67) 



and 



1. 



p^, = — ImD;,. (2.68) 

TT 



The above equation implies that in order to evaluate the spectral function at T we 
need to know the imaginary part of the retarded propagator. We note that the above 
expression for spectral function reduces to its vacuum value as P —>■ oo since the only 
state which enters in the spectral function is the vacuum [98] . 

Now, as mentioned before, the exact retarded photon propagator can be expressed 
in terms of the proper self energy through the Dyson-Schwinger equation: 

where, —i^^^, is the sum of all IPI (one particle irreducible) retarded photon self energy 
insertions which can be decomposed as 

nj,^ vn^ + s^.nf. (2.70) 

Here A'^^ and B^'' , as defined in the previous Section are the transverse and longitudinal 
projection tensors respectively and and Ilf are the transverse and longitudinal com- 
ponents of the retarded photon self energy. The presence of the parameter ^ indicates 
the gauge dependence of the propagator. Although the gauge dependence cancels out in 
the calculation of physical quantities, one should, however, be careful when extracting 
physical quantities from the propagator directly, especially in the non-abelian gauge 
theory. 

Inserting the imaginary part of the retarded photon propagator from Eq. (2.69) in 
Eq. (2.68) we get 

^ ^ (2.71) 



where 



n - ^ ^"""^'^ f2 72l 

~ TT (?2 + Re n«^)2 + (Im n«^)2 ■ ^ - > 
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Using, 



q^q^p,, = 0, (2.73) 



the dilepton rate is finally obtained as 



dR 9/, 2m, I 4m2 , ^^ , ^ 

A ^ "T^^^ ^1^1 '^^ "^ ^ PDfBEiqo). (2.74) 

This is the exact expression for the dilepton emission rate from a thermal medium 
of interacting particles. It has been argued by Weldon [99] that the electromagnetic 
plasma resonance occurring through the spectral function at — —KgYIt^l) could be 
a signal of the deconfinement phase transition provided the plasma life time is long 
enough for the establishment of the resonance. 

Since H^^, and H^j, are both proportional to a (the fine structure constant) they 
are small for all practical purposes. Neglecting Rell^^ and Imllji^ in the denominator 
of Eq. (2.72) one obtains 

1 [2Imnf + Imnf ] 1 imn^/^ 

2pT + PL^ — ^ ^ = j^. 2.75 

TT q^ TT q^ 

This corresponds to the free propagation of the virtual photon in the thermal bath. 
Using Eqs. (2.74) and (2.75) we get 



This is the familiar result most widely used for the dilepton emission rate [82]. It 
must be emphasized that this relation is valid only to O(e^) since it does not account 
for the possible reinter actions of the virtual photon on its way out of the bath. The 
possibility of emission of more than one photon has also been neglected here. However, 
the expression is true to all orders in strong interaction. 

2.2.3 Emission Rate and Current Correlation Function 

The emission rate of dileptons can also be obtained in terms of the electromagnetic 
current correlation function [50] . We denote the electromagnetic current of the strongly 
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interacting particles (quarks or hadrons) by the operator and the leptonic current by 
Jl- As before, e is the couphng between Jl and the virtual photon. For the present, the 
current is taken to contain the coupling constant. The matrix element for dilepton 
production is 

Shi = -ie{ H; l' \ j d''xd^yJl{x)D^'' {x - y) J^iy) \ I) (2.77) 

where D^^ is the free photon propagator. As in the earlier case the leptonic part of 
the current can be easily factored out and we get Eq. (2.46). The photon propagator is 
written in momentum space as 



D^'^ix -y)-J^, e-'^<^-yW^,'^(q) (2.78) 
to obtain 



Shi = -ie 



VV2^m'^''^^'^^ I rf^xe^^-(// I J!:ix) I I). (2.79) 
Squaring the matrix elements and using Eq. (2.48) one obtains the rate of dilepton 
production 

where VF^(g) is the Fourier transform of the electromagnetic current correlation func- 
tion defined as 

/p—pEn 
d'xe^'^-^{H I J';{x)fM I H) (2.81) 

Note that this definition is different from that of Mclerran and Toimela [50] where 

zW)' 

The correlation function is symmetric in // and v. One can use translational invariance 

and the KMS relation to show that 

In this case the dilepton rate is 

dR = LI''' 



H 



g4 (27r)3Ei (27r)3E2' 
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as in Ref. [50]. 

It is seen from Eq. (2.80) that from the measured dilepton and photon distribu- 
tions the full tensor structure of W^" can in principle be determined. This will yield 
considerable information about the thermal state of the strongly interacting system. 

Now, W^j^ is related to the retarded correlator by 

W> = 2e^'''>fBE{qo)l^W^. (2.82) 

where 

W^M) ^ ^ J d'xe'^-e{xo) E {H I [j;(x), fM] I H)^^. (2.83) 

Using the identity Eq. (2.64) and the transversality of the correlation function, the 
dilepton emission rate is obtained as 

5^ = -12^" + f)f^"^^ (2.84) 
We now define the improper photon self energy through the relation 

jjR, ^ j^R, a(5 ^ j^R,a^ pR j^R,./3 ^2.85) 

where —iP^i, is the sum of all self energy diagrams. The advantage is that can be 
defined in coordinate space as [100] 

iP^M^9{xo)^{H I [Ji{x),j':m I H)^. (2.86) 

Taking the Fourier transform and comparing with Eq. (2.83) we see that 

P^M = -KM ■ (2.87) 
Therefore, the dilepton rate can also be expressed as [101] 

It is important to realize that the analysis is essentially nonperturbative up to this 
point. To O(e^) we note that P reduces to the proper self energy 11 (= P DqD~^) 
and consequently Eq. (2.88) reduces to Eq. (2.76). This approximation is the same as 
implied in Eq. (2.75). 
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Let us now try to relate the approach discussed in this Section with the previous 
one. We know that the equation of motion of the photon field is given by the Maxwell 
equation 

which has the formal solution 



A^{x) = j d'yD^\x-y)Jt{y), 



where Dq is the free photon propagator. This can be used in Eq. (2.44) to get Eq. (2.77). 
Again, the connection between the electromagnetic current correlation function and the 

spectral function can be expressed in a straight forward way by substituting and 
from the Maxwell equation in Eq. (2.81) to obtain 

= (gV-(i-r')gM?a)^f (gV-(i-r')») 

= 27r (q^g^^ - (1 - r')g^fe) {q^Qp. - (1 - C"')?/??.) (1 + /b£;).(2.89) 

The gauge dependent terms will not contribute due to current conservation and we have 

W>M = 27rgVM-(l + fBE). (2.90) 

Substituting this in Eq. (2.80) we can recover Eq. (2.63). This establishes the connection 
between the approaches of Refs. [83] and [50] . 

Let us now consider real photon emission from a system in thermal equilibrium. The 
matrix element is given by 

SHi^-i{H-^\ I d'x j;:(x)A^,(x)\I) (2.91) 

where, J^(a;) is the electromagnetic current of the strongly interacting particles which 
produces the photon and ^o(^) is the free photon field. Since the photon escapes 
without re-interacting, the matrix element can be taken to factorize as, 

{H;^\j;:{x)A^,{x)\i) = {H\j;:{x)\i) {im^m. (2.92) 

In the case of a single photon of four-momentum p = {E,p) the free field can be 
expanded as 

A'oix) = e^^-, (2.93) 
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so as to obtain 

Shi = / d'xe^^-{H\Ji{x)\I). (2.94) 

The thermally averaged photon multiplicity is given by 

^,~i:ei5«.^-^. (2.95) 

The sum over the photon polarization is performed using '^e^e'^ — —g^". Proceeding 

pol 

as before we obtain the photon emission rate as 

Note that this expression also follows from the dilepton emission rate given by Eq. (2.80) 
with a few modifications. The factor L^^/q^ which arises from the lepton spin sum 
of the square modulus of the product of the electromagnetic vertex 7* — > the 
leptonic current involving Dirac spinors and the square of the photon propagator is to 
be replaced by the factor — gr^,^/2 from the polarization sum. The factor 1/2 follows 
from the normalization of the photon field (Eq. (2.93)). Finally the phase space factor 
for the lepton pair, dP'pi/[{2Ti)'^ Ei] (i^p2/[(27r)^i?2] is replaced by d^p/[{2'KYE]. 

Using Eqs. (2.82) and (2.87) and the fact that to lowest order in e the improper and 
proper photon self energies are equal, we get 

^^=(l^'""-^-^^^)- ^'-''^ 

In this form, the real photon emission rate is correct up to order in electromagnetic 
interaction but exact, in principle, to all orders in strong interaction. However, for all 
practical purposes one is able to evaluate up to a finite order of loop expansion. It 
is clear from the above that in order to deduce the photon and dilepton emission rate 
from a thermal system we need to evaluate the imaginary part of the photon self energy. 
The Cutkosky rules at finite temperature or the thermal cutting rules [84, 102, 103] give 
a systematic procedure to calculate the imaginary part of a Feynman diagram. The 
Cutkosky rule expresses the imaginary part of the n-loop amplitude in terms of physical 
amphtude of lower order [n — l loop or lower). This is shown schematically in Fig. (2.1). 
When the imaginary part of the self energy is calculated up to and including L order 




y 
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Figure 2.1: Optical Theorem in Quantum Field Theory 
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Y 




loops where L satisfies x + y < L + 1, then one obtains the photon emission rate for 
the reaction x particles — > y particles + 7 and the above formalism becomes equivalent 
to the relativistic kinetic theory formalism [104]. 

For a reaction 1 + 2 — >^ 3 + 7 the photon emission rate is given by [105] 



where 



a 
b 



-^Imin 
-E'2min 



E. 



E 



dR 



(is 



dt 



16(27r)''i? 7(mi+m2)^ Jt______ 

f{E,)f{E,)[l + f{Es)] 



\M\'I 



dEi 



X 



/ 



dEo 



^JaE^ + 2bE2 + c 



(2.98) 



/ I J. 2 2\2 

— [s + t — — mgj 



Ei{s + t — ml — ml){ml —t)+ E[{s + t — ml — ml){s — ml — ml) 
-2mj{ml-t)] 

-Eliml - tf - 2E'i_E[2ml{s + t-ml-ml)- {ml -t){s-ml- ml) 



-E^[{s 



m-i 



777.2) ~ 4777i7772] — (s + t — 7772 — m^)(m2 — t) 



x{s — m^ — 7772) + 777.2(5 + t — 777.2 ~ ''^3) + 'mi{m2 — t) 



{s + t 
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ml) 



AE 
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Eml 
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m^ — t 
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t 



AE 
— ac 



M is the overall degeneracy of the particles 1 and 2, Al is the invariant amplitude of 
the reaction (summed over final states and averaged over initial states), / denotes the 



40 



distribution functions and s, t, u are the usual Mandelstam variables. Now, the rapidity 
of a particle is defined as 

y = o In ■ 



2 E-p, 

where E and Pz are the energy and longitudinal momentum of the particle respectively. 
For massless particles the transverse momentum is 

PT = {E^-Plf/'. 

We then have E — pt cosh y, pz = prsinhy and d^p/E — d^prdy in the case of real 
photons. 

The dilepton emission rate derived above in terms of the photon self-energy can be 
connected by the optical theorem to the kinetic theory rate for a reaction a a — > l'^ l" 
which is given by 

dR f d^Pa f d^Pa J., s f d?Pi f d^p-2 



dK ^ f d-'pa f d^a f d^pi r _(£P2 

d^q J 2£;„(27r)3-^^^"V 2E-a{2TTf^^^~^' J 2E^{2Tif J 2^2(2? 

I-M IL^W- {'^T^Y5''^\Pa+p-a-pi-p2)S^'\q-Pa-Pn). (2.99) 



where f{po) is the appropriate occupation probability for bosons or fermions. The Pauli 
blocking of the Icpton pair in the final state has been neglected in the above equation. 
The transverse mass of the lepton pair is defined as M|. = q^ + M'^ = ~ where 
q^ — M^; M being the invariant mass of the lepton pair. Using the above definition of 
rapidity, we have q^ — M^coshy, q^ — M^sinhy and d^q — MdMd'^Mxdy. In terms 
of the cross-section for the production of a lepton pair the dilepton production rate can 
be expressed as 

^ =/ ^^^^'^^ / ^^^^^^ <e1(Pa,P.)<f;-(Pa,P.) S^'\q - p^ - p,) (2.100) 

where 
and 

{Pa ■ Pa - miy/'^ 



„,aa _ 



rel IT' IT' 

is the relative velocity of the colliding particles a and a. 
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2.3 Emission Rates from Quark Matter 

In this Section wc will discuss the rates of photon and dilepton emission from a thermal 
system composed of quarks, antiquarks and gluons at a temperature T due to various 
processes which are known to contribute substantially to the yield from QGP. Since 
the constituents of the system are massless for all practical purposes, we will encounter 
infrared singularities in the evaluation of the rates. We will discuss how these can be 
screened by summing the Hard Thermal Loops (HTLs) [58, 59] in the theory. The yield 
due to the primary scattering of partons embedded in the colliding hadrons will also be 
discussed since they are expected to constitute the principal background to the thermal 
photons (dileptons) in the region of large transverse momentum (invariant mass) . 

2.3.1 Thermal Photon Emission Rates from QGP 

Naively, one expects that the properties of QGP at high temperature (T >> Tc) can be 
studied by applying perturbation theory due to the small value of the strong coupling 
constant, as{T) which is given by the parametrized form [106] 

"^^^^ = (33-2iV,)ln(8T/Tj- ^^'^^'^ 

However, QCD perturbation theory at high temperature is plagued by infra-red prob- 
lems and gauge dependence of physical quantities, e.g. the gluon damping rate [90, 
107, 108]. The gauge dependence of the gluon damping rate was cured by Braaten and 
Pisarski [58] by an effective expansion in terms of hard thermal loops - i. e. including 
all the relevant loop effects in a given order of the coupling constant in a systematic 
way. The idea of HTL is based on the observation that at non-zero temperature there 
are two energy scales - one associated with the temperature T, referred to as the hard 
scale and the other connected with the fermionic mass ~ gsT {gg << 1), induced by 
the temperature, known as the soft scale. A momentum p^^ appearing in the self energy 
diagram of photon would be called soft (hard) if both the temporal and the spatial com- 
ponents are ~ gsT (any component is ~ T) . If any physical quantity is sensitive to the 
soft scale then HTL resummation becomes essential, i.e. in such cases the correlation 
function has to be expanded in terms of the effective vertices and propagators, where 
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(a) (b) (c) 

Figure 2.2: Lowest order diagrams for photon and dilepton production from QGP. 

the effective quantities are the corresponding bare quantities plus the high temperature 
hmit of one loop corrections. 

However, the problem of infra-red divergences in QCD is not solved completely 
by the HTL framework. The quantities which are quadratically divergent in naive 
perturbation theory such as the damping rate of fast moving fermions in QGP becomes 
logarithmically divergent in effective perturbation theory. On the other hand quantities 
which are logarithmically divergent in the naive perturbation theory turns out to be 
finite if one applies HTL resummation method. The hard photon {E > T) emission rate 
which falls in the second category, is the relevant quantity for the present discussions. 

The thermal photon emission rate from QGP is governed by the following Lagrangian 
density: 

^QGP — ■^QCD + ^jq, (2.102) 



where 



18 Nf _ .a 

^ a=l /=1 ^ 



^-r, = -jF^.F^'^-J^efi^frA^^l^f. (2.103) 
^ /=i 

In the above, is the non-abehan field tensor for the gluon field of color a, ipf is 
the Dirac field for the quark of flavour /, gs is the color charge, e/ is the (fractional) 
electric charge of quark flavor /, A"'s are the Gell-Mann matrices, F^j^ is the electro- 
magnetic fleld tensor and A'^ is the photon fleld. As mentioned in the introduction, the 
dominant processes for photon production from QGP are the annihilation {qq g^) 
and the Compton processes {q{q) — > q{q)^) as shown in Figs. (2.2a) and (2.2b) respec- 
tively. However, the production rate from these processes diverges due to the exchange 
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of massless particles. This is a well-known problem in thermal perturbative expan- 
sion of non-abelian gauge theory which suffers from infra-red divergences. One type 
of the divergences could be cured by taking into account the 'electric type' screening 
through the HTL approximation [58] . The non-abelian gauge theory also contains 'mag- 
netic type' divergences, which can be eliminated if there is a screening of the magnetic 
field [109, 110, 111]. This is in sharp contrast to Quantum Electrodynamics, which is 
free from screening of static magnetic field. However, the study of magnetic screening 
is beyond the scope of HTL approximation as the transverse component of the gluon 
self energy vanishes in the static limit in this framework. Magnetic screening is relevant 
if any physical quantity is sensitive to the scale g^T, where all the loop contributions 
are of the same order [112] and hence the perturbation theory breaks down [113]. The 
production of soft photons {E < QgT) from QGP is non-perturbative because it is 
sensitive to the magnetic screening mass of the gluons [114] and consequently the soft 
photon emission rate is poorly known. The production of hard photons {E > T) is 
insensitive to the scale g^T and hence infra-red divergences can be eliminated within 
the framework of HTL as discussed below. 

Let us try to understand the notion of HTL using massless 0^ theory described by 
the Lagrangian density 

C^\{d<pf-g^\ (2.104) 

The thermal mass (self energy) resulting from the one loop tadpole diagram in this 
model is m\-^ ~ g^T"^- At soft momentum scale (p'^ ~ gT) the inverse of the bare 
propagator goes as ~ g'^T'^- Thus, the one loop (tadpole) correction is as large as the 
tree amplitude. Therefore, this tadpole is a HTL by definition. Braaten and Pisarski [58] 
have argued that these HTL contributions should be taken into account consistently 
by re-ordering the perturbation series in terms of effective vertices and propagators. 
Therefore, according to their prescription we have 

^ = \{d<t>f - 9^' - ImU' + ^ml<l>' = Aff + At, (2.105) 

where Cct — "^th0^/2 is the counter term which should be treated in the same foot- 
ing as the 0^ term. Cct has been introduced in order to avoid thermal corrections 
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at higher order which has already been included in the tree level. With the counter 
term the Lagrangian remains unchanged, so the effective theory is a mere re-ordering 
of the perturbative expansion. A similar exercise has to be carried out in gauge the- 
ory keeping in mind that an addition and subtraction of local mass terms will violate 
gauge invariance. The effective action for hot gauge theories have been derived in 
Refs. [115, 116, 117, 118, 119], whereas the authors of Refs. [120, 121] follow the clas- 
sical kinetic theory approach for the derivation of the HTL contributions. It has been 
shown in Ref. [122] that the contribution of HTL to the energy of the QGP is positive. 
The counter term required to avoid double counting in evaluating the virtual photon 
production from QGP in the two-loop approximation has been derived in [123] recently. 

The photon emission from Compton and annihilation processes can be calculated 
from the imaginary parts of the first two diagrams in Fig. (2.3). Since these processes 
involve exchange of massless quarks in the t/u channels the rate becomes infrared 
divergent. One then obtains the hard contribution by introducing a lower cut-off to 
render the integrals finite. In doing so, some part of the phase space is left out and 
the rate becomes cut-off dependent. The photon rate from this (soft) part of the phase 
space is then handled using HTL resummation technique. The apphcation of HTL 
to hard photon emission rate was first performed in Refs. [55, 56]. For hard photon 
emission, one of the (soft) quark propagators in the photon self energy diagram should 
be replaced by effective quark propagators (third diagram in Fig. (2.3)), which consists 
of the bare propagator and the high temperature limit of one loop corrections [124, 125]. 
When the hard and the soft contributions are added, the emission rate becomes finite 
because of the Landau damping of the exchanged quark in the thermal bath and the 
cut-off scale is cancelled. The rate of hard photon emission is then obtained as [55] 

e'^-^ = l^T^e-^'^H2m2E/glT). (2.106) 

where as is the strong coupling constant. Recently, the bremsstrahlung contribution 
to photon emission rate has been computed [57] by evaluating the photon self energy 
in two loop HTL approximation. The physical processes arising from two loop contri- 
bution (Fig. (2.4)) are the bremsstrahlung of quarks, antiquarks and quark anti-quark 
annihilation with scattering in the thermal bath. The rate of photon production due to 
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Figure 2.3: Two loop contribution to the photon self energy. A diagram interchanging 
the blob in the internal line of the third diagram should also be considered. 





Figure 2.4: Two loop photon diagram relevant for bremsstrahlung processes. The blob 
on the gluon (spiral line) indicates effective gluon propagator. The circle on the vertices 
represent those required to evaluate the imaginary part of the photon self energy in the 
framework of thermal cutting rules (see Refs. ([57]) and also ([103])). 

bremsstrahlung process for a two-flavour thermal system with E > T is given by [57] 

_dR9^P 40 



E- 



aa 



(2.107) 



and the rate due to q — q annihilation with scattering in the thermal bath is given by, 

dRS^^ 40 



E- 



■ aa. 



..ETe-''/'^{JT-JL), 



(2.108) 



d^p 277r5 

where Jt ~ 4.45 and Jl ~ —4.26. The most important implication of this work is that 
the magnitude of the two loop contribution comes out to be of the same order as those 
evaluated at one loop [55, 56] due to the larger size of the available phase space. In case 
of soft thermal photon {E ~ QsT) emission, all the vertices and the propagators have to 
be replaced by the corresponding effective quantities. It has been shown [126, 127] that 
the result is divergent due to the exchange of massless quarks introduced through the 
HTL effective vertices itself. However, such coUinear singularities for light-like external 
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Figure 2.5: Feynman diagram for soft photon production in the SPA. 

momentum could be removed with an improved action [119]. It is also shown that such 
infrared singularities could be removed through KLN (Kinoshita - Lee - Naucnberg) 
theorem [128, 129] by including appropriate diagrams and summing over all degenerate 
initial and final states [130, 131], but the rate is non-perturbative because it is sensitive 
to the scale glT [114]. 

The emission rate of hard photons is well under control within the framework of 
HTL resummation. However, there are important issues in hot gauge theories which 
cannot be addressed within the HTL resummation method [107, 132]. For example, (i) 
HTL resummation is based on the weak coupling limit {qs « 1) to distinguish between 
hard (T) and soft momentum scale {QsT) but such a limit may not be achievable in 
URHIC even for the highest energy to be available at the CERN LHC in the near future. 
Extrapolation of results obtained in HTL approximation to higher values of coupling 
constant will be demonstrated in Section 3.2 of Chapter 3 through the photon spectra, 
(ii) It cannot cure the infra-red divergence problem that arises in the damping rate of 
fast fermions, (iii) it cannot remove the mass shell singularities in the soft photon (real) 
emission rate, (iv) the next to leading order correction to the Debye mass diverges unless 
one includes magnetic screening, which is beyond the scope of HTL approximation and 
finally (v) HTL works for a system in equilibrium ; extension of the formalism to non- 
equilibrium processes is still in the early stages of development. Results from other 
methods such as ladder approximation [133], renormalization group equation [8] etc. 
will be very important in these cases. 

We have studied soft (low transverse momentum) photon production due to parton 
bremsstrahlung in a thermal system within the soft photon approximation (SPA) [134]. 
These photons are produced from one of the external legs in a parton scattering diagram 
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(Fig. (2.5)) i.e. due to processes like ah — > air) where a is a quark or antiquark and 
h can be a quark, antiquark or a gluon. These are hence 0{aa1) processes compared 
to the 0{aas) processes described earher. Naively, the soft photon approximation 
provides that the cross section for the process factorizes into a scattering part and a 
photon production part. The photon can be emitted from any of the external lines in 
Fig. (2.5). The emission of photons from the interior of the scattering vertex (the central 
blob in the figure) is neglected because in the limit of very low energy of the emitted 
photon its contribution is very small. The dependence of the photon momentum p is 
neglected both in the strong part of the matrix element M.q &s well as in the phase 
space delta function. The latter is however taken care of through a correction factor. 
The matrix element is then written as 

A4 = eMo^% (2.109) 
and the rate of production of soft thermal photons at temperature T is given by 



x$(s, S2, ml ml)K,{z) E—, (2.110) 

where z^m = {ma + mb)/T, z = ^/s/T and gab is the colour and spin degeneracy. The 
cross-section for the process ab — > cd'j is given by 

da' _ a a{s) 

(2-111) 



with 

ais) 



-\{s,ml,ml)/s 

Here J is the real photon current and 



r dt^^^ (E'\e-J\l^J. (2.112) 



where A(x, z) — — 2{y + z)x + {y — z^. The strong interaction differential cross- 
section daqq/ dt and daqg/ dt for scattering of quarks and gluons are obtained from semi- 
phenomenological expressions used earlier by several authors for this purpose. The 
emission rates of hard and soft thermal photons are shown in Fig. (2.6). One observes 
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Figure 2.6: Thermal photon emission rates from QGP at T— 200 MeV. 

that the 2-loop corrections to the photon rate namely, the bremsstrahlung contribution 
given by Eq. (2.107) (the long-dashed curve) and the contribution from qq annihilation 
with scattering as given by Eq. (2.108) (dot-dashed line) are of the same order as the 
1-loop result (dotted line). Soft photons calculated in the SPA is plotted only up to 
E —1 GeV for consistency [134] regarding the Landau-Pomeranchuk effect [135]. 

2.3.2 Hard QCD Photons 

As remarked earlier, the large transverse momentum window for the photons will have 
a contribution due to partonic interactions which occur when nuclei start to interpen- 
etrate at very high collision energy This contribution can be understood in terms of 
perturbative QCD, and can provide reliable information about the partonic distribu- 
tion in the nuclei, as well as a means of providing near-absolute normalization of the 
photon measurements at larger transverse momenta. Since these photons are emitted 
at the very early stages of the collision with large momenta they are usually referred 
to as hard or prompt QCD photons. The cross section to leading order of perturbative 
QCD for the production of a photon in a hadronic collision, is obtained by convoluting 
the cross-section for the elementary processes e.g. Compton and annihilation with the 
gluon or quark contents of the participating hadrons [136]. Neglecting the correction 
due to neutron-proton asymmetry the nucleus- nucleus {A — B) collision is then built up 
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as an incoherent sum of independent nucleon-nucleon {N — N) collisions. The N — N 
cross section for the QCD Compton process at y = is given by 

da'^°^^ aag dxaXaXb 



/: 



X 



mXa) + q[Xa)\g[Xb) — h [Xa ^ Xb) 



xlxf 



(2.114) 



where, q{x) and g{x) are the quark and gluon structure functions of the nucleon respec- 
tively and 

_ XaXT 
2Xa - Xt 

XT = 2pr/ 

XT 



Xr] 



For the annihilation process we have 



da 



ann 
N 



d^PTdy 



Xt 



dxa 



Jxmin Xa - {xt/2) 



X 



X 



Xi, 



q{Xa)q{Xb) 3 + {Xa ^ Xb) 



/y>3 <T^3 

•^a-^b 



(2.115) 



A large contribution to the prompt photons also comes from the fragmentation 
of partons in the final state in a parton-parton scattering, ah cd. These are 
bremsstrahlung processes and though of higher order {0{aal)) compared to Comp- 
ton and annihilation {0{aas)), is found to contribute in the same order of magnitude 
at all values of the transverse momenta. Using an effective structure function the cross 
section ior N — N collisions is obtained as [136] 



da''^""' 
d^PT dy 



27rs2 A2 XT J XT {yr/'^.y 



2/t/(2-j/t) Xa - yTf^ 



X 



F2(Xa)[G{Xb) + -Q{Xb)] 



+ iyT/2f 



XI 



{Xa ^ Xb) 



(2.116) 



where Xb = XayT/ {^Xa - yT) and 



F2{x) ^ x^e^q[q{x) + q{x)], Q{x) ^ x^[q{x) ^- q{x)] and G{x) ^ xg{x). 
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The strong coupling is given by 

127r 

= (33-2iV,)lnWVA^c„)' 
with g2 ^ p2^^ ATy = 4 and Aqci? = 0.23 GeV. 

Neglecting the effects of shadowing, the photon distribution in a coUision of nuclei 
A and B at an impact parameter b due to hard scattering of partons is given by 

E^W=TMt)E^ (2.118) 

where Tab{^) is the nuclear overlap integral. It is defined as 

TAB{h)^ I d''sTA{s)TB{\b-s\). (2.119) 

Ta = J dz pa{z, s) is the nuclear thickness function and pA is the nuclear number density 
normalized to the mass number A. Also, Tab is normalized such that 

J d%TAB{h)^ AB. (2.120) 

Let us now consider central [h — 0) collisions of two identical nuclei of mass number A. 
We will assume a constant nuclear number density po; so that 

4 , 

-'kR'aPo = A 



and 



Therefore, 



Ta{s) = 2po{Rl - s'fl\ 



Taa{Q) = / d^sTA{s)TA{s) 
Jo 

_ 9 A^ 

Writing Taa(O) ~ A'^/ttR^ we obtain 

dN^ = ^da^_ 
d^Prdy t^R\ d^pTdy 
In Fig. (2.7) we have plotted the photon yield due to the three processes described 

above using the MRSD-/ [137] set of structure functions. One observes that the yield 

from bremsstrahlung is of the same order as the sum of the Compton and annihilation 

processes. We have used K — 2 in order to take into account the higher order QCD 

corrections. 
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Figure 2.7: QCD (prompt) photons from Pb-Pb collisions at 200 A GeV. 

2.3.3 Thermal Dilepton Emission Rates from QGP 

The most dominant contribution to thermal dilepton production from QGP comes 
from quark anti-quark annihilation as shown if Fig. (2.2c). To obtain the differential 
number of lepton pairs emitted per unit four- volume we start from Eq. (2.100). Since 

quantum effects in the QGP is known to have a negligible effect we will use Boltzmann 
distribution for the quarks. The phase-space integrals are performed using the delta 
function to obtain 
dR 



1 4m^ 

Integrating over Mr we arrive at 

dR 1 . ^mf" 



(2.122) 



dW^dy 4(27r)4' M2 

and on integration over y we get 
dR 1 



ox ,..^M'T' / Mcosh|/\ , _ , 
1 - ^) a,,-(M)— 1 + — ^ exp(-Mcoshy/r) 



cosh y 



T J 



m p,).^'r(l-3^).«(M)if.(M/T). 



(2.123) 



(2.124) 



The cross section for two light quark flavours can be evaluated using the diagram 
Fig. (2.2c) to obtain [22] 



SOtt c? / 4m2y/' / 



2m2 



(2.125) 
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where m is the lepton mass and rUq is the quark mass. The thermal mass of the quarks 
will be used in the calculations. It is worth mentioning that the ct^ corrections to this 
rate comes from the Compton and annihilation diagrams i.e. the diagrams labelled a 
and h in Fig. (2.2) with a virtual photon in place of the real one [138]. One encounters 
mass (collincar) singularities in the evaluation of these processes which can be cured 
through the use of KLN theorem. 

A significant contribution to the dilepton yield in the lower invariant mass region 
comes from soft (bremsstrahlung) processes. The rate of production of soft thermal 
dileptons can be evaluated analogously as the soft photons using the SPA. The photon 
emitted from the external lines in Fig. (2.5) in this case is a virtual one which eventually 
produces a lepton pair. The rate is given by [139] 

,,.0 , = -rr^ / dz / ^-KMTdMT 

dM^ay lOTT^ J{ma+mt+M)/T JM 

^ «' '^-^s,S2,ml,ml)K,{z) ^^^^^^^^ (2.126) 



where 



r4 V > a> 0/ .V ^dM^ffMrdy 

da^^^~ o? d{s) 



dM-^d^MTdy 127r3M2 Mlcosh^y' {^-V^^) 
d{s) is as given in Eq. (2.112) with J appropriately modified to account for virtual pho- 
ton emission. The rates for thermal dilepton production from qq annihilation (0(q;^)) 
is plotted in Fig. (2.8) along with the soft 0{p?o?^ contributions discussed above. 

2.3.4 Drell-Yan Dileptons 

In the region of large invariant mass the dilepton yield is dominated by the Drell-Yan 
process A ^ A ^ + X and thus forms the principal background to dileptons 
emitted from the QGP. Here, a quark from one of the nucleons in nucleus A annihilates 
an antiquark from one of the nucleons in the other nucleus to produce a virtual photon 
which subsequently decays into a lepton pair. The differential yield of such lepton pairs 
produced is A — A collisions is obtained by an incoherent sum of the contributions from 
independent nucleon-nucleon coUisions. For central collisions, the Drell-Yan yield is 
given by 
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Figure 2.8: Thermal e'^e and 77 emission rates from QGP at T—200 MeV. 

where as in the case of prompt photons, q{x) and q{x) are the quark distribution 
functions of a nucleon. At y = 0, Xa = Xf, = M/ ^/s where s is the square of the center 
of mass energy of the coUiding nucleons. 



2.3.5 Diphotons 

We will briefly discuss the production of diphotons produced due to the annihilation of 
quarks and antiquarks. The importance of thermal diphoton emission lies in the fact 

that an experimental detection of large mass diphotons can possibly provide a vahiable 
confirmation of the results obtained from the measurement of dileptons [140]. The 
diphoton cross section is given by: 



27ra=^ iVe i2S + if V \ , 
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1+ " 



M2 









1/2- 






1 + 








1 2ml 


M2 







4m^ 
M2 



1 _ ^ 
M2 



1/2- 



(2.129) 



In the above Nc — 3, S = 1/2 and Cg is the charge of the quark and we have used 
TTiq — mth — (27rQ;s/3) T. The emission rate is obtained by using this cross-section 
in Eq. (2.123) or Eq. (2.124). The rate of emission of thermal 77 pairs from QGP at 
r=200 MeV is shown in Fig. (2.8). 
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Figure 2.9: QCD dielectrons and diphotons from qq annihilation. 

Diphotons with large invariant mass are also produced from the hard QCD anni- 
hilation of quarks and antiquarks in the colliding nuclei. The yield can be calculated 
analogously as the Drell-Yan pairs to get [140] 



Ki y — Xa — xi, — M/y/s. An arbitrary cut-off on the momentum transfer Pc{— 2 
GeV) has been introduced so that perturbative QCD remains valid in this case. In 
Fig. (2.9) we have shown the dilepton yield due to the Drell-Yan process for Pb-Pb 
collisions at 200 A GeV. Also shown is an estimate of the diphoton yield in such a 
collision. As before we have used the MRSD-/ set of nucleon structure functions. 




X [QMq{Xb) + {Xa ^ Xb)] . 



(2.130) 



Chapter 3 

Medium Effects and Emission Rates 
from Hot Hadronic Matter 

In this Chapter we will consider photon and dilepton emission from a thermal system 
of interacting hadrons. We have seen that the photon and dilepton emission rates are 
related to the imaginary part of the photon self energy in the medium. As a result 
the rates will depend on the in-medium modifications of the hadrons appearing in the 
internal loop of the photon self energy diagram. Here the hadronic medium consists 
of mesons and baryons at a finite temperature. Due to the interactions with real and 
virtual excitations, the properties of these hadrons are expected to get modified. As 
a result the propagators appearing in the photon self energy undergo modifications. 
These are studied in the framework of Thermal Field Theory. 

Many of the hadrons in a hot hadronic gas are electrically charged and hence couple 
to the electromagnetic field. Pions and p mesons form the most important constituents 
of such a system. This is because pions are light and the p mesons have large spin- 
isospin degeneracy. Lagrangian densities constructed with the tt, p, a;, r] and ai fields 
have been used to calculate the amplitudes for photon production. Of these, the medium 
modifications of the p and ui mesons are known to affect the photon and dilepton spectra 
significantly. Though we have included the presence of nucleons and antinucleons for the 
evaluation of the medium effects of vector mesons, we have neglected their contribution 
to the production of photons and dileptons. Also, we will assume the net baryon number 
to be zero. In Section 3.1 we will study the in-medium modifications of hadronic masses 
and decay widths using well known models. Thereafter, in Sees. 3.2 and 3.3 we shall 
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discuss the emission rates of photons and dileptons where these medium modifications 
are taken into account. 

3.1 Hadronic Properties at Finite Temperature 

In this Section we will consider the in-medium modifications of the mass and decay 
widths of the p and uj mesons. We will discuss the Quantum Hadrodynamic (QHD) 
model, the gauged linear and non-linear sigma models, and the hidden local symmetry 
approach. We will also show how the QCD sum rules can be used to constrain the 
spectral functions of the vector mesons in the medium. 

The change in the hadronic mass in the medium can be understood from the fol- 
lowing phenomenological arguments [141]. Let us consider the propagation of a vector 
meson in a nuclear medium. The attenuation of the amplitude at a distance z, in a 
Fermi gas approximation, is given by e""*^^, where n is the density of nucleons and a 
is the meson-nucleon interaction cross section. The optical theorem relates a to the 
imaginary part of the forward scattering amplitude; a — A7rlmT{E)/k. It then follows 
that the meson wave function ip ~ exp[2mnzT{E) / k], where the imaginary part of T 
accounts for the attenuation and the real part modifies the dispersion relation of the 
propagating particle. In terms of an effective mass {rrics = m + Am), the propagation 
can also be described by -0 ~ exp[i^ E'^ — mlf[z]. Comparing the arguments of the 
exponential we get 

^^^_27m^ (3.1) 
m 

This relation clearly shows that the enhancement or reduction of hadronic masses de- 
pends on the sign of ReJ^{E). 

3.1.1 Quantum Hadro dynamics 

In the Quantum Hadrodynamic model [71, 72] of nuclear matter the vector meson 
properties are modified due to coupling with nucleonic excitations. The discussion has 
two parts. We will first study how the properties of nucleons are modified in matter 
at finite temperature. The nucleons interact through the exchange of scalar a and the 
vector cu mesons and their mass is modified due to the scalar condensate. Thereafter, 



57 



we will consider the changes in the p and uj meson masses due to coupling with these 
modified nucleonic excitations. 

a) The Nucleon Mass 

The interaction in QHD is described by the Lagrangian 

C>QHD = -9.NN A^7m ^ ^'^ + 9cNN NaN, (3.2) 

where N{x), 0"(x), and u}{x) are the nucleon, a, and ui meson fields respectively. The 
cr{u!) field couples to the scalar (vector) current of the nucleon with the coupling constant 
gaNN{gu}NN) which wiU be specified later. 

As discussed in the previous Chapter (Section 2.1), the free nucleon propagator at 
finite temperature and density in general has four components. The time-ordered i.e. 
the (ll)-componcnt is physically relevant for our purpose and we will denote this as 
G°(p) where p denotes the four- momentum of the nucleon. So we have 



+ 2niS{p'^ - M%)r]{p.u) 



_p2 -M% + ie 

= G%{p) + G'j,{p), (3.3) 

where the first term (G'p) describes the free propagation of nucleon-antinucleon pairs 
and the second term {G%) allows for the on-shell propagation of particle- hole pairs. Mjv 
in the above equation is the free nucleon mass. 

The effective mass of the nucleon in matter at finite temperature in presence of 
the interaction described by Eq. (3.2) will appear as a pole of the effective nucleon 
propagator. In the Relativistic Hartree Approximation (RHA) [71, 72] one obtains the 
effective propagator by summing up scalar and vector tadpole diagrams self-consistently 
i.e. by using the interacting propagators to determine the self energy. The effective 
propagator referred to as the Hartree propagator is given by 

G^(p) = G'ip) + G°(p)E^(p)G^(p). (3.4) 

This is pictorially shown in Fig. (3.1). Here S^(p) is the nucleon self energy which 
contains contributions from both scalar (S^) and vector (E^) tadpole diagrams [71, 72] 
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Figure 3.1: Diagrammatic representation of Dyson-Schwinger equation for nucleons in 
RHA 



and is given by 



where 



and 



ml J (27r) 
~ m2 J (27r)4 



Tr[7^G^(p)]. 



(3.5) 



(3.6) 



(3.7) 



Here, (?7iu) is the mass of the neutral scalar (vector) meson. The solution of Eq. (3.4) 
now reads, 

1 



G^ip) = (^+M* 



Nl 



p2 - + ie 
^ Cf (p) + G'g(p) 



- + 27ri5{p' - M*^')r}{p.u) 



(3.8) 



One observes that the pole structure of the effective nucleon propagator in RHA re- 
sembles that of the non-interacting propagator with shifted mass and four-momentum 
i.e. p = p + Sf^ and = Mm + , is the effective mass. Using G% in place of the 
full Hartree propagator in Eqs. (3.6) and (3.7) defines the Mean Field Theory (MFT) 
values of the self energies. This is equivalent to solving the meson field equations with 
the replacement of the meson field operators by their expectation values which become 
classical fields i.e. cr — > {a) and u — > {u). This yields 



(3.9) 
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which indicate that the nuclear ground state contains scalar and vector meson conden- 
sates generated by baryon sources, ps is the baryon density of the medium and Ps 
is the (Lorentz) scalar density. The spatial part of the uj condensate vanishes due to 
rotational symmetry in infinite nuclear medium. These condensates are related to the 

scalar and vector self energies generated by summing tadpole diagrams in QHD as 



The mean field approximation is thus to neglect the fluctuations in the meson flelds 
which themselves are generated by the nucleons. 

RHA is obtained when one includes the vacuum fluctuation corrections to the MFT 
results. This amounts to the inclusion of the Dirac part of the propagator Gp in the 
calculation of the self energies. Summing over the vacuum tadpoles results in a sum 
over all occupied states in the negative energy sea of nucleons. Vacuum (or quantum) 
fluctuations, as these are called, form an essential ingredient in a relativistic theory of 
many particle systems. Since there are inflnite number of negative energy states in the 
vacuum one expects that the vacuum contribution to the self energy is infinite. 

Let us now find the Hartree self energy of the nucleon with the full nucleon prop- 
agator consisting of a medium as well as a vacuum part. The vector part of the self 
energy is obtained from Eq. (3.7) as 



The first term of this equation appears to be divergent. The usual procedure is to 
regularize the integral in n dimensions by dimensional regularization to render the 
integral finite. One can then shift the integration variable from p to p. The resulting 
integral vanishes by symmetric integration. The vector self energy then reduces to 



9aNN{cr) 




(3.10) 




(3.11) 



^v" = -glNN^^'PBlml 



(3.12) 



and gives rise to a an effective chemical potential. 



(3.13) 
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The scalar part of the self energy follows from Eq. (3.6): 



X 



ml J (27r)4 p2 _ ^ ^2 J (27r)3 E 

fMl^*,T)+fMl^*,T) 



(3.14) 



where 



^"'^^^^'^^ ~ exp[(^* -\*)/T] + 1 
fFD{fi*,T) = exp[(£;* + /i*)/T] + l 



(3.15) 

The baryon density of the medium is given by 

PB-J^sj d'p[fMl^\T) - hoil^^T)]. (3.16) 

The first term in Eq. (3.14), to be denoted by E^^\ represents the contribution to 
the scalar self energy from the filled Dirac sea and is ultraviolet divergent. We will 
now proceed to renormalize this divergent contribution. The first step is to isolate the 
divergent part through dimensional regularization. This gives 

= _ ^giVivr(2-n/2) 
* ml 27r2 

= _^£(lz^(M|, + 3M^Ef + 3M^,Ef + Ef) (3.17) 

smce Mj^ = Mjv + Ef . The divergence in E^ now appears as the pole of the F-function 
for physical dimension n = 4. The counter terms needed to remove the divergent 
contributions from the loop corrections to the measurable amplitudes are 

>CcT = E"nCTVn! (3.18) 

n=l 

Including the contributions from the counter terms the renormalized self energy becomes 
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where 



^CTC 



,1 2/1 ) (3-20) 

The coefficients (aj) are fixed by defining a set of renormalization conditions. Since the 
scalar density ps {={iIj4')) is not a conserved quantity the tadpole diagrams appear in 
the self energy. The tadpole contribution must vanish in normal vacuum (free space) i.e. 
((t)o = 0. This is ensured by the term aia in jCct- Oi2cr^ is the meson mass counter term 
which ensures that is the physical (measured) mass. Since the original Lagrangian 
of QHD [71] does not contain and o"'^ terms, three and four point meson amplitudes 
must vanish at the tree level. The last two counter terms in Eq. (3.18) are chosen to 
maintain this condition at zero external momenta for the a meson when nucleon loop 
corrections are included. We thus have 



«n = -i{-gaNNnn - 1)! / -^Tr[G^(p)1. 



Consequently the effective nucleon mass reads 



(3.21) 



NN 



N 



+ 



(27r)3 E* 



M, 



N 



- MUM*^ - Mjv) 



- ^M^(M^-M^)2-H(M^_M^)3 



(3.22) 



The solution of this equation (with glj^j^ = 54.3 and ma = 458 MeV) gives the effective 
nucleon mass as a function of temperature and baryon density. At zero baryon 
density it can be parametrized as [142] 

^ 8.94n 



N 



1 - 0.0264 



0.16 



(3.23) 



where T is in GeV. We thus observe that in nuclear matter, scalar((T) and vector(c(j) 
mean fields induced by nucleon sources give back-reactions to the nucleon propagation 
itself and modify its self energy. This is the origin of < in QHD. 

In Fig. (3.2) we depict the variation of nucleon mass with temperature for a set of 
baryon densities. We observe that the nucleon mass falls steadily with density for a 
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Figure 3.2: Variation of nucleon mass with temperature for different values of baryon 
densities. The normal nuclear matter density p^=0.1484 fm~^. 




Figure 3.3: Same as Fig. (3.2) with (sohd) and without( dashed) vacuum fluctuation 
corrections to MFT. 
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Figure 3.4: The vector meson self energy due to NN polarization. The double lines in 
the loop indicate nucleon propagators in RHA. 

fixed temperature. However, the variation with temperature for given values of baryon 
densities shows interesting features. At zero baryon density the nucleon mass decreases 
monotonically as a function of temperature, but for finite densities it increases slightly 
before falling. This trend is similar to that obtained by Li et al [143], and may be 
attributed to the modification of the Fermi-sea at finite temperature and density. Our 
calculation shows a 35% reduction of the effective nucleon mass at T = 160 MeV and 
two times normal nuclear matter density compared to its free mass. In order to highlight 
the effect of vacuum fluctuation corrections we compare the MFT results with those 
obtained using RHA. This is plotted in Fig. (3.3). We observe that the effect of vacuum 
fluctuation is substantial for higher values of the baryon density. The contribution of 
the antinucleons from the Dirac sea is responsible for such an effect. 

b) The Vector Meson Meiss 

In a medium, meson properties get modifled due to its coupling to nuclear excitations 
as shown in Fig. (3.4). This modiflcation is contained in the meson self energy which 
appears in the Dyson-Schwinger equation for the effective propagator in the medium. 
The interaction vertices are provided by the Lagrangian 

= 9VNN (a^T.t^A^KT - ^^^/^-^"^^'^«'^) ' (3.24) 
where V^'^ — {u!^,p^}, Mjv is the free nucleon mass, is the nucleon field and Ta — 
{1, r}, T being the Pauli matrices. 

The lowest order contribution to the vector meson self energy is expressed in terms 
of the self-consistent nucleon propagator described in Eq. (3.8). This is given by 

W^ik) = -2ig'y^^ I ^Tr [r^{k) G^{p)n-k) G^{p + k)] , (3.25) 
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where F'^ represents the meson-nucleon vertex function obtained from Eq. (3.24) and is 
given by 

ri'ik) = for u 

r>^(k) = ^^^ + ^-Jl^a'^-K■, for p (3.26) 

where cr'^" = |[7'*,7"]- The vector meson self energy can be written as a sum of two 
parts 

n'^-(A;) = U''/{k) + Il%'{k), (3.27) 



where 



U^/{k) = -2ig'yj,^ J -0^Jr[r^{k)G^{p)r'^{-k)G^{p + k) 



U^^ik) ^ -2t9'yj,^ I ^Jr[T'^{k)G^{p)T^{-k)G^{p + k) 
+r^^ik)G^{p)r''i-k)G^ip + k) 

+ r^{k)G^{p)T^{-k)G^{p + k)] . (3.28) 

n^'^ is the vacuum polarization. This is a bilinear function of Gp and hence describes 
the correction to the meson propagators due to coupling to NN excitations. The NN 
pairs can be excited only if the four-momentum carried by the mesons is in the time- 
like region (A;^ > 0). Hence the shift in the mass of the vector mesons due to vacuum 
polarization is caused by processes like V — > NN — > V where N represents nucleons in 
the modified Dirac sea having an effective mass M^, smaller than what it would be in 
free space. From Eqs. (3.28) and (3.8) we have 



(27r)4 {p2-M*^^)[ip + k)^-M\ 



*2] 
N . 



Prom naive power counting it can be seen that this part of the self energy is ultraviolet 
divergent and has to be renormalized. A few comments about renormalizability of the 
interaction given by Eq. (3.24) is in order here. At very large momenta the propagator 
for massless bosons ~ 0{k~'^), whereas for massive vector bosons it goes as ~ 0(1). 
This poses severe problems to the renormalizability of the theory with massive vector 
bosons. However, in a gauge theory with spontaneous symmetry breaking the vector 
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gauge bosons acquire mass in such a way that the renormahzabihty of the theory is 
always preserved. The theory which involves neutral massive vector bosons coupled to 
a conserved current is also renormalizable. This is because in a physical process the 
propagator D(^^ — {—g^'^ + k^k" /m'^)/{k'^ — + ie) appears between two conserved 
currents and J^, and the offending term k^k^ /m^ docs not contribute because of 
current conservation {k^j^J^ = 0), making the theory renormalizable. This is the case 
for the cu meson [144, 145] which we shall consider first. The counter term required in 
this case is 

JC^m^-IcV^^V,^. (3.30) 

We use dimensional regularization to separate the divergent and the finite parts. The 
divergences now appear as a pole in the F-function at the physical dimension n — A. 
The renormalized vacuum polarization tensor for the cu is then given by 



U'P'ik) = {g^'' - k^'k''/k^)W^''{k^), (3.31) 



where 



n'^^fc") = iimL[Y{2-n/2) C dzzil-z) 

TT'' I Jo 

- £ dzz{l - z)lii[M;^ - k'^z{l -z)]^-C (3.32) 
in which the counter term contribution 

has been included. ( is now determined by the renormalization condition 

nr(A:')|Mj^-.M. = 0. (3.34) 

Finally, we arrive at 

n-(^2) = \Re{wrr, 



^%v^2 dzz{l-z)ln 



_ j^2^^i _ ^) 



(3.35) 



- k^z{l - z) 

Because of the tensor interaction in Eq. (3.24) the vacuum self energy for the p meson 
is not renormalizable. We employ a phenomenological subtraction procedure [146, 147] 
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to extract the finite part using the condition; 







(3.36) 



M*-*Mn 



with (n = 0,1, 2,.. ..cx)). Using dimensional regularization and the above subtraction 
procedure we arrive at the following expressions: 



9pNN 



/i + M; 



. 1 



N 



2M, 



-I2 + 



N 



where 



Jo 



dz z(l — z) In 



2 '2Mjv 
M*^ - e z{l - z) 



f 

Jo 



dz In 



Ml - z{l - z) 
M*^ - k'^ z{l - z) 



(3.37) 



(3.38) 



(3.39) 



Ml-k'^z{l-z) 

The medium dependent part of the polarization, n^J", describes the coupling of the 
vector mesons to particle-hole excitations. It contains at least one on-shell nucleon 
propagator which provides a natural ultraviolet cut-off in the loop momenta. This 
part of the self energy leads to an increased effective mass of the vector mesons in the 
medium. 

We recall that in a hot and dense medium because of Lorentz invariance and current 
conservation the general structure of the polarization tensor takes the form 

n'^'^ = IlT{ko, k)A'"' + UL{ko, k)B'"' 

where the two Lorentz invariant functions 11^ and Ul are obtained by contraction as 

k^ 



1 



ti^ being the four velocity if the thermal bath. 

In the case of a vector meson of four-momentum k interacting with real particle-hole 
excitations in the nuclear medium these are obtained as 



I 



Tr 



^VNN I ^27r)4 

(n^'^ + n^'''* + n^'*)^^ 



r^{k)Gp{p)r,{-k)GD{p + k) + {F^D) 



(3.40) 
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with 



9vNN J_ f Pdp 
- 8|p1|^| 



UJr, 



{k' + 2M*jI) In 



(A;2 + 2|p1 1^1)2 -4fcoV 



2,. ,2 
P 



{k''-2\p\\k\Y-Aklujl 



(3.41) 



and 



7r2 ^ \2MnJ \k\ J ujp 



{e + 2\p\\k\? - 

(A;2-2|p1 1^1)2 -4A;gcu2 



X 



fMf^*,T) + fFD{f^*,T) 



9VNN 

47r2 
X In 

X 





2 

) 




" pdp 






1^1 


Up . 


+ 2 






1 7.2 ,2 



(A;2-2|p]|A;|)2-4A;ga;2 



-4|p1|^| 



The longitudinal component of the polarization tensor is given by 



(3.42) 



(3.43) 



(3.44) 



with 



f yjvjv k'^ f pdp 



47r2 



/ 



UJr, 



{(fco - 2a;,)^ - |fc|2^i, fe^-2feoc^P + 2|p1|fe| 

A;2-2A;ou;p-2|p1|A;| 



X 



A;2 + 2A;oa;p-2|p1|A;| 
fFDil^*,T)+fMl^*,T)\ , 



(3.45) 



9VNN 



TT 



\ k 



X 



2M;v/ |A;| 
/F^(/i*,T) + /^^(/.*,T) 



(F-2|p1|A;|)2-4A;gu;2 



and finally 



9VNN 



27r2 V2M 



'AT 



2 /• pdp 



2 A;^ {k"^ -2kQUJpf' 



2iPr-y- 



2|^F 



(3.46) 
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Figure 3.5: Transverse and longitudinal dispersion relations of the p and uj mesons. The 
solid and dashed curves pertains to the transverse and longitudinal modes, respectively. 



X In 



k'^ -2kQUjp + 2\p\\k\ 



X In 



X 



k"^ -2kQUp-2\p\\k\ 
k"^ + 2kQUjp + 2\p\\k\ _ A\^kl 
k'^ + 2koijjp-2\p\\k\ \k\ 

fFD{l^*,T) + fMl^*,T) 



{k^ + 2^00;^)^ 



(3.47) 



where Up = fr^ + and fpD stands for Fermi-Dirac distribution functions for the 
nucleons. In the above the superscripts 'v\ 'vV and 't' represent the vector- vector, 
vector-tensor and tensor-tensor components respectively arising from the product of 
vector and tensor terms in Eq. (3.26). The dispersion relation for the longitudinal 
(transverse) mode now reads 



kl - -ml + Renf(2.)(A;o, k) + Rell^{k^) = 0. 



(3.48) 



It is important to point out that the self-energy functions of the nucleon as well as that 
of the vector mesons calculated in this section are actually the 11-component of the 
2x2 self-energy matrix in the real time formulation of thermal field theory. Since here 
we are concerned about the real(dispersive) part of the self-energy in our discussion on 
effective masses we do not make any distiction between these and the scalar self energy 
function following Eq. (2.19) of Chapter 2. 
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In Fig. (3.5) we plot the dispersion relations for p and uo mesons at T = 180 MeV 
and twice normal nuclear matter density. This is obtained by solving Eq. (3.48) for 
the p and uj mesons. One observes a small difference between the longitudinal(L) 
and transverse(T) modes in case of the to meson but in case of the p this sphtting is 
negligible (attributable to the smaller vector coupling constant). We have observed 
that the quantity — along the dispersion curve remains almost constant ~ my 
which is defined as the value olkl oXk — Q on the mass hyperbola. This means that a 
simple pole approximation of the p and uj propagator at k'^ — my is good enough for 
our calculations. The splitting between the transverse and longitudinal components of 
the self energy of vector mesons with both vector and tensor interactions can be shown 
to be [148], 



fpD + fpD 



X 



(/ cos'^ — r cos + tr 



(3.49) 



where 



C + 8poko\p\ \k\ cos - 4|p]2 \k\2 cos^ 0_ 

u = ?,ki\p\^-\km 

V = 4/coPo|p1 1^1 

w = 2pl\kf + \kf\p\' - kl\p\' 

C — 1^1^ — k^p^. 

The following values of the coupling constants and masses [144, 146] have been used in 
our calculations so as to reproduce the nuclear saturation density: Kp — 6.1, fi'pjvAr — 
6.91, mp = 770 MeV, Mn = 939 MeV, = 0, and ^^^^v = 102. 

Thus the physical mass (my) is defined as the lowest zero of Eq. (3.48) in the limit 
k —>■ 0. In this limit — n£ — U^, and we have. 

In;; = n = + (3.50) 



where 



n^(^o,^-0) = -^^ / p'dpF{\p\,M*^)[fMf^*,T)+fMf^*,T)] (3.51) 
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Figure 3.6: Variation of p meson mass with temperature for various baryon densities. 



with 



n\p\,M*r,) 



1 



+ 



3'2M, 



N 



(3.52) 



and a;^ = + M*^ . The effective mass of the vector meson is then obtained by solving 
the equation: 



kl-m\^ Ren = 0. 
The effective masses take the following parametrized forms [142] : 



(3.53) 



1-0.127 



1 - 0.0438 



^T(GcV)' 

. 0.16 ; 

'T{GeV) 
0.16 



5.24 



7.09 



(3.54) 



The effective mass of the ai meson has been estimated from that of the p mass 
using Weinberg's sum rule [149]. One finds reference to two other kinds of masses in 
the hterature. The invariant mass is defined as the lowest order zero of Eq. (3.48) with 
neglected. Again, the screening mass of a vector meson is obtained from the pure 
imaginary zero of the quantity on the left hand side of the same equation with ko — 0. 
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Figure 3.7: Variation of cu meson mass with temperature for various baryon densities. 



These two masses are different because of the non-analyticity of the polarization tensor 
at the origin i.e. at ipo,p) = (0, 0). 

In Fig. (3.6) the effective mass of the p meson is plotted against temperature for 
various values of baryon density. We observe that the variation of the p mass follows 
quahtatively the same trend as that of the nucleon. In this case, the p mass decreases 
by 45% at T= 160 MeV and two times normal nuclear matter density compared to its 
free space value. This is due to the fact that the large decrease of the modified Dirac 
sea contribution to the p self energy dominates over the in-medium contribution which 
is seen to increase with temperature. We then evaluate the effective cu mass with the 
values of the coupling constants mentioned above. The results are plotted in Fig. (3.7). 
The quantitative difference in the p and uj meson masses is due to the different numerical 
values of the coupling constants e.g. the tensor interaction is absent in case of the cu 
meson and quite significant for the p meson. 

Before we proceed further a few comments on the QHD model calculations are in 
order. In this model the major contribution to the change in the masses of the p and u 
mesons arises from the nucleon-loop diagram. For the dressing of internal lines in matter 
we restrict ourselves to the Mean Field Theory (MFT) to avoid a plethora of diagrams 
and to maintain internal consistency. It has been shown [101, 105] that the change 
in the p mass due to pinr interaction is negligibly small at non-zero temperature and 
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Figure 3.8: Decay of p meson. 

zero baryon density. Therefore the change in the p meson mass due to pmr interaction 
is neglected here. At finite baryon density, the dynamics is more involved due to 
the medium effects on the pmr vertex, the pion propagator coupled with delta-hole 
excitation, and the coupling of the p meson with A^*-hole excitations [79, 150, 151, 152, 
153, 154, 155, 156]. The effect of such medium modifications is to broaden the p-peak 
as well as to produce complicated structure around the peak. Here we have restricted 
our calculations within the realm of MFT, i.e the internal nucleon loop in the p and cu 
self energy is modified due to tadpole diagram only. 

c) The Vector Meson Width 

The physical decay width of an unstable particle is related to the imaginary part of its 
self energy. For a particle at rest the width is obtained from the relation 

Imn(A;o) = koT{ko) (3.55) 

where ko is the energy of the decaying particle. There are various ways to calculate 
the imaginary part. The Cutkosky rules at finite temperature [84, 102, 103] provides 
a simple and systematic way to calculate the imaginary part of the self-energy. Again, 
use can be made of the fact that the self-energy function develops cuts along the real 
axis when the particles in the internal loop become on-mass shell. The discontinuity 
across these cuts is pure imaginary for real ko so that we have 

Discn(A;o) = [n(A;o + ie) - U{ko - ie)] = 2iImn(A;o). (3.56) 

Let us consider the p meson width. The imaginary part of the p self energy is totally 
dominated by the pion loop which is borne out by the fact that the two pion decay 
mode of the p has a branching ratio of ~ 100 %. In contrast, the real part of the p self 
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energy which is responsible for the mass modification has a neghgible role to play as far 
as the pion loop is concerned. It is found to cause a small positive shift of the p pole. 
The p — TT interaction is described by the Lagrangian 

£}^,^-g^,^- {7^x8^7^). (3.57) 

The 11-component of the self energy of the p meson due to pion loop (Fig. (3.8))is given 

by 

-iKl{k) = -g%^ / (2P - kY^^^.M (2P - kyi^l.iP - k), (3.58) 
where the pion propagator is 

A?i(g) = ^-^---2m6{q'-ml)fBE{\qo\) 

l + fBEi\qo\) fBEi\qo\) 



q2 _ ^2 _|_ q2 _ ^2 _ 
After integration over po, the imaginary part of the self-energy function is obtained as 



(27r) 



(2p - kY{2p - kf 



I (1 + fBE{0Jp) + fBE{uJp-k)Y{ko -LUp- UJp^k) + 
{fBE{i^p-k) - fBEM)Siko -UJp+ UJp^k) I + ^ I 



(3.59) 



where use has been made of Eq. (2.20). Using the relations given in Eq. (2.29) the 
longitudinal and transverse components can now be worked out. The terms involving 
the thermal distribution functions in the above equation can be interpreted in terms of 
pion absorption from and emission into the medium. The first and second ^-functions 
correspond to time-like and space-like regions of k respectively. Restricting to the time- 
like region, we define the spin-averaged quantity, 

a'^" 1 1 

^[A^arnUr + B^JmUL] = -[2ImnT + ImH^] = -Imn^^ 



487r 
X In 



2T 

1 + ^ 



\k\W{k^) 



l-exp[-f(fco + |A:|M^(A:^))] 
l-exp[-l{ko-\k\W{k^))] 



(3.60) 
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Figure 3.9: The p — > tttt decay width as a function of temperature for different values 
of baryon densities. 



where W{k^) — y^l — Aml/k'^. In the rest frame of the p {k = Q) this reduces to 



Imn;^(A;o) _9%n.,2 . 2^3/2 



Using Eq. (3.55) the p decay width is obtained as 



^fBEif) + 1 



(3.61) 



^p{ko 



487r 



^■2 



(3.62) 



with fBE{x) — [e^ — l]"^. It is interesting to note that the phase space factor (2/be + 
1) when written in this form clearly shows that the in-medium width is actually the 
difference between the rates of decay and formation of the resonance. 

In Fig. (3.9) we demonstrate the in-medium effect on the decay width of p meson. 
The observed enhancement of the decay width with temperature at non-zero values of 
the baryon density is solely due to the stimulated emission of pions in the medium; a 
consequence of the (1 + Jbe) terms in the decay width. This is just a manifestation 
of the well known Bose enhancement (BE) effect which is more clearly observed in 
Fig. (3.10). 

The effect of the modifications in the mass and decay width of an unstable particle 
is conveniently illustrated through the dimensionless quantity Ay which is basically the 
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Figure 3.11: Spectral function of p meson in QHD. Solid (long dashed) line corresponds 
to T =180 MeV (T =150 MeV). The spectral function in vacuum is shown by the 
dotted line. 
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Figure 3.12: Diagrams contributing to the width of the u meson, 
spectral function apart from some factors. It is defined as 



(3.63) 



where M is the invariant mass of a lepton pair and niy — my — ReU. The p spectral 
function Ap, in units of e is plotted in Fig. (3.11). The shifts in both the spectral 
functions towards the lower invariant mass region correspond to the reduction of their 
masses due to thermal interactions. 

Let us now consider the cu meson. The vacuum width (8.5 MeV) of the lu is known 
to be dominated by the a; — > 37r mode. A substantial contribution to the cu width 
also comes from the process con ^ tttt in a thermal bath [148, 157]. The Feynman 
diagrams are shown in Fig. (3.12). We shall use the Gell-Mann Sharp Wagner (GSW) 
interaction [158] given by 

Cgsw = ^ e^^apd^uj^d'^p^TT (3.64) 

for the ujpn vertex and the Lagrangian given by Eq. (3.57) for the pTTir vertex. The 
cu ^ St: width is obtained as 

ra;-»37r(^o) dz dx \ F \^ S (3.65) 

where S is the phase space factor, given by 

^ = [(1 + fBE{E,)){l + /s£;(^2))(l + fBEiEs)) - f be{Ei) f be{E2) f be{E^)] 

and 
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Figure 3.13: The u; width as a function of temperature. 
The hmits of integration are 

Zmax = {kl - 3ml)/2ko, 

Xmax = \/0-5ko{z - z.^^){z'^ -ml)/{2kQZ - kl - ml), 




= z, 

E2 = x+{ko-z)/2, 

E3 = -x + {ko-z)/2, 

I I = \/Ei - ml, 

where Pi is the pion 3-momentum. The amphtude for the process is 

\FMp,\'\p,ni-Z',)H 

where 

^ _ A;2 + m2 - 2koiEi + E2) + 2E1E2 
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i=l 

1 

dl2 + miri 
1 

1 

2((ii2(ii3 + m^r^)/ii/i2 

2(^13^23 + m^rj)/l2/i3 

2(^12(^23 + m^r^)/ii/i3 
(Ei + E2r-p1-mJ 

(£;i + £;3r-p1-mj 

(E2 + E3)=^-p?-m2. 

The contribution from the reaction cuir — > tttt to the decay width of the cu is calculated 
analogously. 

In Fig. (3.13) we have shown how the decay width of the cu meson increases with 
temperature. The narrow peak of the u in vacuum is broadened substantially due 
to interactions with the thermal pions. Both the modes discussed above are seen to 
contribute almost equally to a ten-fold broadening of the cu in the medium. This 
is also observed in the spectral function of the u; shown in Fig. (3.14) where F^^ = 

3.1.2 Models with Chiral Symmetry 

In this Section we will discuss finite temperature effects on the vector meson properties 
considering a few models which respect chiral symmetry. The effects of in-medium prop- 
erties of vector mesons on the electromagnetic probes will be presented in Section 3.2. 
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Figure 3.14: Spectral function of u; meson in QHD. Solid (long dashed) line corresponds 
to T —ISO MeV {T —150 MeV). The spectral function in vacuum is shown by the dotted 
line. 

a) The Gauged Lineeir Sigma Model 

The hnear sigma model (LSM) is a convenient tool to describe the low energy dynamics 
of pions, because it shows explicitly how the spontaneous symmetry breaking (SSB) of 
global chiral symmetry (SU(2)0SU(2)) by the isosinglet sigma field generates pions, the 
Nambu-Goldstone (NG) bosons. However, there are reservations about the description 
of the a meson as a well-defined degree of freedom because of its large decay width 
which is comparable to its mass. But, it has been argued [44, 159] that in the limit 
of chiral symmetry restoration, the decay of a into two pions should be disallowed as 
a and tt become degenerate in mass in this limit. It has been explicitly shown that 
the width of a due to cr — > 27r decay vanishes as T ^ T^, where is the critical 
temperature for chiral transition. Though it is still not known whether or not the 
critical temperature for deconfinement and chiral transition are the same; presently we 
will make no distinction between them. 

The simplest version of LSM contains isosinglet a field and isotriplet pion field and 
respects the charge conjugation, parity and time reversal symmetry (CPT) [42]. The 
Lagrangian obeying these constraints is 

jClsm = tr I d^^ I' +11^ tr | $ [' +^Xtr{\ $ - htr{^) (3.66) 
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where $ is defined as 



$ = i(o- + ZTT • f) (3.67) 



with r being the Pauh matrices. The non-zero value of h ensures that the pions are 
massive and consequently PCAC (partially conserved axial current) relation is satisfied. 
Note that | $ a^+n^ is chirally invariant. Elimination of the a field by imposing the 
condition a'^ + — results in the Non-Linear Sigma Model, which will be discussed 
in the next Section. 

In the gauged LSM [42, 66, 160] one introduces the vectors and their chiral partners 
(axial vector) through left and right handed fields as 

Vr^{fr + a'^) • t+ {cu^ + /r) (3.68) 

VJ' = (p^ - a>') ■ t+ {u'' - /O (3.69) 

where Oi and /i are the chiral partners of the p and uj mesons respectively and t = t/2. 

The inclusion of axial vector mesons will increase the number of possible couplings 
and hence the number of arbitrary parameters become large. In order to include vector 
mesons in LSM with minimal coupling to the matter fields tt and a, one requires [161], 
that the Lagrangian and its chiral transformation properties are such that the current 
generated by the chiral transformation is proportional to the vector field itself. This 
leads to the field-current identities and eventually the idea of VMD [162]. The field- 
current identity is achieved by promoting the SU(2) global chiral symmetry of vector 
fields to a local gauge symmetry as was done by Yang and Mills [163] for the isospin 
symmetry. The Lagrangian for the vector field reads, 

U = hr I Fr r +ltr I Fr r +mltr[{Vn' + {V^^f] (3.70) 

where Fl^^ — V/^^ — V/^j. — ig V/^^, V{^^ . In the above Lagrangian the kinetic term 
for the gauge fields remains invariant under the transformation and the field-current 
identity is obtained through Gell-Mann Levy theorem from the mass term of the gauge 
fields as 



Jt,r - -'-fv,^ (3.71) 
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We note that chiral symmetry is a global one in QCD. Therefore, the local symmetry 
has to be broken and this is achieved by the mass term of the vector fields in the 
Lagrangian. Next, one has to introduce the interaction of matter fields (tt and a) 
and the gauge fields preserving the field-current identity. Noting that the ordinary 
derivatives occurring in Eq. (3.66) spoils the field-current identity, we introduce the 
required interactions consistent with the gauge principle i.e. by replacing the partial 
derivatives by covariant derivatives: 

D^'^ = - ig{Vi^^ - $ y/). (3.72) 

Finally the Lagrangian density for the gauged LSM is obtained from Eqs. (3.72), (3.66) 
and (3.70) as 

Cgism = tr I \^+ixHr I $ \^+^Xtr{\ $ \Y - htr{^) + U. (3.73) 

Expanding the kinetic term for the matter field one finds 

1 



tr I D^^ 



|2_ _ 

' ~ 2 



(3.74) 

The above equation indicates that (i) a shift in the a field (u — > do + a) gives rise to 
mixing between tt and fields (a term ~ gaodf^n ■ arises from the second term of 
the r.h.s. of Eq. (3.74)), which has to be eliminated by an appropriate shift in the Oi 
field, (ii)there is no interaction term involving cu, which can only be introduced through 
anomaly and (iii) the kinetic term for pion gets modified because of the shift in the Oi 
field. Thus to get back the canonical form of this term one has to renormalize the pion 
field TT l^fZ^., where = rr?plrv?^^. Zt,^ = 1/2 gives the Kawarabayashi - Suzuki - 
Riazuddin- Fayyazuddin (KSRF) relation [164, 165]. After some algebra one gets [66], 

= /i/(Z^(To), 
= /i/(7o + 2A(7o, 



/tt = \] Zt^gq. (3.75) 

Taking = 137 MeV, = 600 MeV, rup = 770 MeV and rua^ = 1260 MeV, we get 
(70 = 152 MeV, h = (102 MeV)^, = 412 MeV and A = 7.6. 
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With these inputs the thermal masses of p and oi mesons to lowest order in g at 
low temperatures are obtained as [66] 



ml{T) 



-3 



2 , 

45m2 



+ 



2m^ 



ai 



(3.76) 



(3.77) 



In the chiral limit ctq goes to zero and many of the couplings vanish. Assuming 
the validity of VMD in the medium Pisarski has showed that [166] p and oi become 
degenerate with a mass value ~ 962 MeV {i.e. p mass increases). On the other hand, if 
one adopts a scenario where vector meson dominance (VMD) is not valid in the medium 
then mp(T^) = maj^{T^) = 630 MeV {p mass decreases). However, it is important to 
mention at this point that chiral symmetry can also be realized via the Georgi limit [167] 
where the p meson becomes massless. Pisarski [42] has argued that the results obtained 
by Georgi in the non-linear sigma model can be translated in terms of the gauged linear 
sigma model without the validity of VMD, for which there is no unique prediction for 
the behaviour of the p mass at non-zero temperature. Thus the behaviour of in-medium 
p depends on the validity of VMD in the medium. 

b) The Gauged Non-Linear Sigma Model 

It is well-known that the global SU{2)i<^ SU {2)r symmetry of two-flavour QCD is 
expected to be spontaneously broken to the subgroup SU{2)v and the pions appear 
as the N-G bosons. The non-linear sigma model with SU{2)i^ SU{2)r/ SU{2)v is an 
effective theory of QCD for the description of pion dynamics. The in-medium properties 
of vector mesons have been studied by Song [67, 168] in the framework of gauged 
non-hnear sigma model (NLSM) [61]. We will discuss this model briefly because it is 
very similar to the gauged LSM; the main difference is that the a degree of freedom is 
eliminated in NLSM by the non-linear realization of chiral symmetry as mentioned in the 
previous Section. We start with the observation that a perfectly valid parametrization 
of the matter field $ could be 



U — exp 



27_^^ 



= exp 



2t 



-0 



(3.78) 
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where = (p'^r'^ / \/2 is the pseudoscalar field and = ^J^f-^- The Lagrangian for the 
NLSM based on the manifold SU{2)i SU{2\/SU{2)v is given by 

J TT , 



Co 



-tr 



(3.79) 



The vector and the axial vector fields can be introduced as the Yang-Mills gauge 
fields as before to minimize the number of arbitrary parameters in the model. The 
resulting Lagrangian is given by 



-NLSM 



tr 



1 

tr 

2 



1 2 



--tr I Fr 



+mitr[{vrr + (vm 



(3.80) 



where V^''' = (v^ ± af^)/2, and denote vector and axial vector fields. To improve 
the phenomenology of the model, the following higher dimensional terms can be added 
to the Lagrangian [67, 169] without spoiling the symmetry under consideration 



Cmm = -i^tr \d^UD,U^ F'^^"^ + D^U'' DM F^'^"" 



(3.81) 



where ^ is a constant determined form the decay of vector mesons [67]. The thermal 
shift of the p-mass evaluated with pion loop, pion tadpole and pion-ai loop resulting 
from the interaction given by Eqs. (3.80) and (3.81) shows negligible change in the 
p-mass from its vacuum value [67]. 

The effective masses of p, ai and uu at non-zero temperature have also been evalu- 
ated [168] with a SU (3); SU (3)^ symmetric Lagrangian: 



NLSM 



^tr IdMD^'U^ - -tr I I' --tr I F^ P 
+ml triiVrr + + -Jl tr [m{U + C/t - 2) 



itr \Df,UDM^F^''"'+ D^U^DMF 



+Ktr 



Fl^F^'^M^], 



(3.82) 



where U is defined as in Eq. (3.78) with the Pauli matrices r" replaced by the Gell-Mann 
matrices A". The two higher dimensional terms with co-efficients ^ and k are added 
to improve the phenomenology. It may be noted that although these terms retain the 
gauge invariance of the model the renormalizability of the model is spoiled. 
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The dynamics of the uj meson is governed by the anomalous interaction, also known 
as Wess-Zumino interaction given by 

ll'anamaly = e ^,^^d^ u'' tr{&' (3.83) 

This is very similar to the Gell-Mann Sharp Wagner [158] interaction already described 
before. 

The following values of the parameters consistent with the vacuum properties of the 
vector and axial vector mesons have been considered [168]: ((?, k, ^) = (10.30, 0.34, 0.45) 
and (6.45, —0.29,0.06), referred to as set 1 and 11 respectively. The calculation of the 
thermal mass shift of the vector and axial vector mesons with these inputs reveal that: 

(i) for parameter set I p and cu masses increase with different rate and oi mass decreases, 

(ii) for parameter set II the thermal mass shift of p and cu is negligibly small but Oi 
mass decreases shghtly. 

c) The Hidden Local Symmetry Approach 

In case of the two chiral models described above the vector mesons are introduced as 
Yang- Mills field and the mass term for the gauge bosons are put in by hand which is 
not entirely satisfactory. In the hidden local symmetry (HLS) approach the p meson 
is generated as a dynamical gauge boson of a hidden symmetry in the NLSM [170, 
171]. It has been explicitly shown that, in general, any NLSM corresponding to the 
manifold G/H is gauge equivalent to a "linear" model having G global ® Hiocai symmetry. 
Accordingly, the Lagrangian of Eq. (3.79) can be written in a form that exhibits, besides 
SU{2)i^ SU{2)r global, a local SU{2)v symmetry - the hidden symmetry and the p 
meson appears as a gauge boson corresponding to this symmetry. (The axial vector Oi is 
not included in the minimal version of HLS Lagrangian.) To make it more explicit, one 
introduces two SU{2) matrix- valued variables ^i{x) and $,r{x) with the transformation 
properties [170], 

ii^r{x) ^ h{x)ii^r{.x)gl (3.84) 

with 

U = eje. (3.85) 
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where h{x) e [SU{2)v]iocai and gi^r e [SU{2\r] global- ii,r is parametrized as 

= exp[iS(a;)//s T W/^] (3-86) 

where tt = tt" and S = i". The unwanted degrees of freedom E, which have entered 
through Eqs. (3.85) and (3.86) are known as the "compensators"- the would be N-G 
bosons which have to be "eaten up" by the hidden gauge boson, p. These extra degrees 
of freedom then reappear as the longitudinal polarization of the (massive) p. 
Now the covariant derivatives are defined as 

D^,ir = dnir - igV^.ir + ^Cr , (3.87) 

where /^(r^) is the external field corresponding to the gauging of SU{2)i SU{2)r and 
is the gauge field corresponding to the symmetry [SU{2)v]iocai- With these fields 
two types of [SU{2)i^ SU {2)r\giobai®[SU {2)v]iocai invariants can be constructed [170, 
171, 172] which are 



P 



D^ii ■ Ci - D^ir -il' . (3.88) 



A linear combination C, — La + aLy is equivalent to the original Lagrangian given 
in Eq. (3.79). By fixing the gauge = C,r = exp(i7r//7r) and hence eliminating the 
unphysical degrees of freedom, E) one can show that Ca = Lq while Cy vanishes when 
the equation of motion for is used. So far, has been treated as an auxiliary field. 
It is assumed that the kinetic term for this field is generated by quantum effects or by 
QCD dynamics. The full Lagrangian with the kinetic term is 

Chls ^jCa + aCv - Iq^uQ^" (3.89) 

where g^i, is the non-abehan field tensor for the p meson. The Lagrangian of Eq. (3.89) 
is then written as [170] , 

Chls = ^(9,7?)^ + ^ag^ ■ tt x d,7f + \g'afy^ - \e,.r + (3.90) 
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Figure 3.15: Shift in the pole position of the p spectral function for gauged Linear and 
Non-Linear Sigma Models and Hidden Local Symmetry Lagrangian at T = 150 MeV. 

The above equation implies that the mass of the p meson (m^ = ag^f^) is generated 
due to SSB via Higgs mechanism and the unphysical N-G modes S (and not tt) are 
"eaten-up" by the gauge boson i. e. the three extra degrees of freedom get converted to 
the three degrees of polarization appropriate for the massive gauge boson. For a — 2 
one recovers the KSRF relation. This value of a also results in universal couphng of p. 

Harada et al [68] have evaluated the finite temperature effects on the p-mass upto 
one loop order in the HLS approach due to the thermal pion and p meson interactions. 
Their results reveal that at high temperature the reduction in p mass due to pion loop is 
overwhelmed by the increase due to thermal p loop contribution, although the net shift 
is rather small. The contribution of thermal pions to the p self energy in this model 
is different from other calculations because in HLS approach there is no pion tadpole 
contribution. 

In Fig. (3.15) the shift in the pole position of the p-spectral function is depicted 
for the Linear Sigma Model (LSM), Non-Linear Sigma Model (NLSM), and Hiddden 
Local Symmetry (HLS) approach. For both the NLSM and HLS interactions the p mass 
increases by an amount 90 MeV and 10 MeV respectively. Due to the enhancement of p 
mass in the NLSM, a larger phase space is available for the decay process p — > tttt and 
consequently the p appears to be broader in this case compared to HLS interaction. On 
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the other hand, in the gauged LSM the p mass reduces by about 45 MeV at T = 150 
MeV. It may be noted that p mass decreases in gauged LSM for low temperatures and 
increases for temperatures in the vicinity of the chiral phase transition. 

3.1.3 Spectral Constraints at Finite T 

We will now discuss the medium modifications of hadronic properties in the QCD sum 
rule (QSR) approach. We will briefly discuss the basic principles of QSR in vacuum [76, 

173, 174, 175, 176, 177] and at finite temperature [178, 179, 180] and see how these ideas 
can be utilized to constrain the spectral function of the vector mesons p and u). 

a) QCD Sum Rules at T = 

Because of the complex infrared properties of QCD it is very difficult to extract infor- 
mation on hadronic spectra from the QCD Lagrangian. The idea of QSR formalism is 
to approach the bound state problem in QCD from the asymptotic freedom side i. e. to 
start at short distances and move to larger distances where confinement effects become 
important, asymptotic freedom starts to break down and resonances emerge as a reflec- 
tion of the fact that quarks and gluons are permanently confined within hadrons. The 
breakdown of asymptotic freedom is signalled by the emergence of power corrections due 
to non-perturbative effects of the QCD vacuum and are known to be more important 
than higher order ag corrections [173]. These are introduced via non- vanishing expec- 
tation values of quark and gluon condensate operators such as, (0|gg|0), (0|G^j,G'^''"|0) 
etc. where q{x) is the quark field and is the gluon field tensor. In standard per- 

turbation theory these matrix elements vanish after normal ordering. In the following 
we will discuss how the QCD sum rule approach connects the perturbative and non- 
perturbative domains and leads to the determination of hadronic resonance parameters 
like mass, coupling constant etc. in terms of the condensates. 

The QCD sum rule approach starts with the Wilson operator product expansion 
(OPE) of a product of suitable currents. The gluon and quark condensates appear 
as higher dimensional operators in the expansion. The coefficients of this expansion 
contain the short distance part and the long range part is contained in the vacuum 
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expectation values. The coefficients can be evaluated perturbatively in terms of the 
parameters (cts and the quark masses) of the Lagrangian used. In our discussions on 
QSR in vacuum, we shall consider the time-ordered or causal current correlator 

W^M = i I d'xe^'^-^ (0|r{J^(x) J.(0)}|0) (3.91) 

which has a tensor structure of the form 

KM = (^M- - (3-92) 

The source (electromagnetic) currents can be defined in terms of the quark fields (in 
units of e); 

2 1 - 1 

J,x = -u-f^u - -d-ff,d - -s7m«- (3-93) 

Defining the current in the p, uj and channels as 

= (l/2)(«7^«-J7^(i), 
= (l/2)(iZ7^ii + d7^d), 

Jt = ^"7m^, (3.94) 

one can express the electromagnetic current in terms of p, cu and fields as, 

J,-J', + lj:i-ljt- (3.95) 

Presently, we will confine our discussions to the p meson (J^'^ = 1 ,7 = 1) only. 
The analytic structure of the correlator (W), for spacelike = —q^, can be expressed 
through a dispersion relation: 

W{Q') = - I I"^^(^) ^" + (subtraction). (3.96) 

TT J S + Q 

The imaginary part of W is proportional to the spectral density which can be modelled 
as consisting of a conspicuous resonance and a continuum with a sharp threshold Wq, 

Imiy (s) = TT ^ 5{s - ml) + ^ (l + — ) e{s - u^) (3.97) 

Res 8^ ^ ^ ^ 

Qr is the resonance strength and the pole position is at m|j. 
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/ 



The theoretical side of the sum rule is derived from an operator product expansion 
for large — —q^ (deep Euclidean region) where asymptotic freedom is realized. Thus 
we write 

" d^xe''i-^T{J^{x)Mx)}) = Ci{q) + J2 Cn{q)On (3.98) 

n 

where / is the identity operator, C's are the Wilson coefficients, and OnS are the local 
gauge invariant operators constructed from the quark and gluon fields. The operators 
are ordered by their increasing dimensions and therefore, the coefficients fall off by 
corresponding powers of g^. On dimensional grounds one observes that the operators 
of dimension d > leads to power corrections. However, for large — —q^ a 
fewer number of power corrections {d — 6) is sufficient to converge the series. Taking 
the vacuum expectation value of Eq. (3.98) we obtain [173] 

W(Q^) = (^1 + -ij In ^ + — {0\m^uu + mjd\0) 

1 rv TTfy — 

TTCy — 

-^,mm,X^u + dj,yd) Yl Ql,>^''m- (3.99) 

^ q=u,d,s 

The left hand side (l.h.s.) is the well known two point Greens function which can 
be expressed in terms of the phenomenological parameters characterizing the strong 
interaction processes via the dispersion relations, consistent with the current under 
consideration. The right hand side (r.h.s.) has been evaluated by using OPE in the 
short distance (asymptotic freedom) region. The vacuum expectation value of the higher 
dimensional operator appears as a power correction to the asymptotic contribution (the 
first logarithmic term in the r.h.s. of the above equation). 

The sum rule therefore, becomes (modulo subtractions) 



TT J 



™^ = W(Q^). (3.100) 

In Eq. (3.100) the r.h.s. corresponds to large or small distance scale with fewer 
power corrections and l.h.s should be saturated by the lowest resonance, which is a long 
distance phenomenon. Therefore, in order to get a balance between the two sides we 
would like to have a weight function which enhances the low contribution relative to 
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the high one. This can be done by taking additional derivative with respect to Q^, 
and then taking and the number of derivatives n to infinity. This yields the Borel 
transformed sum rule. Borel transformation is equivalent to the following mathematical 
operation: 



Ml 



where 



lim 

Q2/n=M^=const. 



[n 



1)! 



Q 



2n 



d 



(3.101) 



(3.102) 



and Mb is the Borel mass. Applying Eq. (3.101) on the l.h.s. of Eq. (3.100) and 

Eq. (3.102) on the r.h.s. and expressing vacuum expectation value of four fermion 
operators in terms of two fermions, we obtain [173] 



B lmW{s) ds = 



1 

Stt 
+ 



M 



B 



1 + — + 



Stt 



TT 



j^{0\mqqq\0) 



3Mf 



(0|^G^,G^-|0) 

B TT 

448 TT^a, 2 



(3.103) 



Substituting the various values of the matrix elements as given in Ref. [173] we obtain 



/ 



Stt 



, c^. , 0-04 

^ TT ^ 



0.03 



M% 



(3.104) 



Differentiating with respect to l/M^ we obtain another sum rule: 



/ 



-^/^l lmW{s) sds = — M| 
Stt 



, a, 0.04 0.06 
1 J i \ 

TT Mi M% 



(3.105) 



In Eqs. (3.104) and (3.105) the terms and M^^ arise due to gluon and quark 

condensates respectively. Assuming that r.h.s. of Eq. (3.97) is saturated by the p 
resonance we get from Eqs. (3.104) and (3.105) 



Ml 



(1 + as/n) fl - (1 + cuo/M|) e-'^«/^B] - 0.04/M| + O.Oe/M^ 



1 + a,/7r) [1 - e-^^/^B] + 0.04/M4 - 0.03/M6 



(3.106) 



The above expression still depends on the Borel mass Mb and the continuum threshold 
ujq. The value of ujq can be inferred from the data of e^e" annihilation. The absolute 
value of p mass is then obtained by looking for the stability plateau i.e. choosing Ml 
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such that dmp{M^) / dM^ = 0. To determine the resonance strength for the p meson 
we keep only the p resonance in the sum of Eq. (3.97) and substitute it in Eq. (3.104) 
to obtain, after an elementary integration, 

a, 0.04 0.03 



47r6?.= 



2t: m?p 



(3.107) 



where Qp = l/g^p. 

Eqs. (3.106) and (3.107) indicate how the resonance parameters of vector mesons 
can be extracted by using QCD sum rules in vacuum. In the next Section we will briefly 
discuss the QCD sum rules at non-zero temperature. 

b) QCD Sum Rules at Finite T 

As mentioned earlier, it is not the causal (time-ordered) but the retarded current cor- 
relator has the required analytic properties in a thermal system. QCD sum rules for 
vector mesons in medium [179, 180] start with the retarded current correlation function, 

W^il(?o,g) =^ |ciW<^-^^(xo)([J^(x), J,(0)]) , (3.108) 

where = (go, O) is the four momentum and the currents are defined in Eq. (3.95). 
As discussed earlier there are two independent invariants in medium, the transverse 
(W^) and the longitudinal (W^) components of the polarization tensor both of which 
satisfy fixed q dispersion relations. These are defined through 

W^, = -q'{A^.Wi + B^^WE) (3.109) 

where A^^i, and Bf^i, are defined by Eqs. (2.26) and (2.27) respectively. In the limit 
g — > 0, as there is no spatial direction, we have 

= = VF^ = W;^J{-2,ql) (3.110) 



where the last relation follows from the trace of Eq. (3.109). In this limit, 

V? - ql 

Rel^^ can be calculated using perturbation theory with power corrections in the deep 
Euchdean region q^ — > — oo using OPE. For example, OPE for Rel^^(go), which is the 



1 /-oo ^rr\\Al^(n\ 

ReM/^(go) = -P / du^ J + (subtraction). (3.111) 

TT Jo — qn 
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same as the OPE for the causal (Feynman) correlator W^{qo), has a general form at 
ReW\ql ^ -oo) ^ -Colng^ + g ^""^""'^^Y^^^' ^^'^\ o^i^^^)), , (3.112) 

n=l ^ 

where /i is the renormalization point of the local operators which separates the hard 
scale \Q\ and soft scales such as Aqcd and T. C„ are the c- number Wilson coefficients 
which are T independent. All the medium effects are contained in the thermal average 
of the local operators On- Since {On)^ ~ T"^' • A.qqjj with l + m — n due to dimensional 
reasons, Eq. (3.112) is a valid asymptotic expansion as long as ^ and Aq^j^. 
The local operators On{lJ^^) in the vector meson sum rule arc essentially the same with 
those in the lepton-nucleon deep inelastic scattering (DIS) and can be characterized 
by their canonical dimension {d) and the twist (r=dimension-spin) . They are given in 
Ref. [180] up to dimension 6 operators and we will not recapitulate them here. For 
g — s> 0, Eq. (3.112) is an asymptotic series in l/Q"^ or equivalcntly an expansion with 
respect to d. The medium condensates {Onin'^)) j. may be evaluated by low energy 
theorems, the parton distribution of hadrons and lattice QCD simulations. 

Matching the left and right hand sides of Eq. (3.111) in the asymptotic region 
— > — oo is the essential part of QSR. This procedure gives constraints on the spectral 
integral and hence the hadronic properties in the medium as well as in the vacuum. 
There are two well-known procedures for this matching, namely the Borel sum rules 
(BSR) [177] and the finite energy sum rules (FESR) [181], which can be summarized as 



roo 

Jo 



dql V{ql) [W«(go)- W<f,,(go)]-0, (3.113) 





V{s) 



ql'^Oiuo-ql) (FESR), 
e-9o/A^| (BSR). 



Here ImlV^^^^(go) is a hypothetical imaginary part of obtained from OPE and Mb 
is the Borel mass. 

We have seen in the previous Section that in QSR in the vacuum, the spectral 
function [i.e. Iml^^ in Eq. (3.111)) is usually modelled with a resonance pole and the 
continuum to extract the mass and decay constant of hadrons. In the medium, such 
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a simple parametrization is not always justified because of the thermal broadening of 
the spectrum and also because of the new spectral structure due to Landau damping 
and the thermal mixing among mesons. Therefore, the model independent constraints 
obtained from QSR are only for the weighted spectral integral. 

For example, the first three finite energy sum rules at finite T read [180] 

POO 

h = [W^(go)-ImTy^^^^(go)]rfgo = 0, (3.114) 

h = [\mW^{qo)-lmW^^^{qo)]qldql = -C2m,, (3.115) 

h = [lv^W{qo)-luvW^^^{qo)]qUQl-C^{Oz)^. (3.116) 

Similar sum rules hold for the axial vector channel (in the chiral limit) except that one 
has a different operator for O^. One can also generalize the above sum rules to finite q 
[182, 183]. 

Explicit forms of Cn(Pn)j. have been calculated as [180] 

Co = -^(1 + ^), Ci = 0, (3.117) 
C2(02). = ^A'^G^), + \{Sqi^oDoq)r, (3-118) 

16 

C^{0^)^ = -(scalar 4- quark))^ + y(5gi7oL>oL>oL>o?)T- (3-119) 

Here we have neglected the terms proportional to the light quark masses and the quark- 
gluon mixed operators. S is used to make the operators symmetric and traceless. At 
low T, one may use the soft pion theorems and the parton distribution of the pion to 
estimate the r.h.s. of the above equations. When T is close to T^, one has to look for a 
totally different way of estimation; the simplest approach is to assume a resonance gas 
to evaluate the r.h.s., while the direct lattice simulations will be the most reliable way 
in the future. 

The sum rules U can be used to check the vahdity of the calculations of the spectral 
functions using effective theories of QCD. This is in fact quite useful for the spectral 
function at finite baryon density. At finite T, especially near the critical point, the 
behavior of the condensates with d > 4 is not known precisely. Therefore, it is rather 
difficult to make a strong argument on the spectral constraints near at present. The 
future lattice simulations of these condensates are highly called for. 
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c) Parametrization of the Spectral Function 

We will now introduce a parametrization of the correlator at finite T. The parametriza- 
tion should be consistent with the experimental data from e+e^ hadrons. It should 
also be consistent with the high energy behaviour known from perturbative QCD at 

As the vector mesons appear as resonances in the electromagnetic correlator, using 
Eqs. (3.108) and (3.95) we can write 

ImWj = ImW^'' + UmW^f + ^ImlVj^ (3.120) 

which shows that the contributions of uj and mesons to the electromagnetic probes 
are less by almost an order of magnitude compared to p mesons. In the limit g = 0, we 
have from Eq. (3.110) 

Im W;^^{q^) = -2,ql lmW\q^). (3.121) 

Our next task is to parametrize ImVK^(go) which is now a positive dimensionless quan- 
tity. We take a Breit-Wigner form with an energy-dependent width for the resonance 
along with a continuum: 

MH^^ta, f = 0) = ffj^^^^^f^^, + 3^(1 + ^)y^^s- (3.122) 

where indicates the in-medium values of the parameters. At T = 0, the above form 
reduces to a relativistic generalization of the parametrization used by Shuryak [75] to 
fit the experimental data of e'^e" — > hadrons. Here Imll^ is the imaginary part of 
the self-energy which should in principle contain all the channels which can destroy 
or create a p in the thermal bath. Hence ImH^ is the difference of the decay-width 
and the formation width and is given by Imll^ = qoV^qo). However, we have seen 
that for baryon free matter the most dominant contribution to ImH^ comes from the 
pion-loop [148]. The quantity ujq in Eq. (3.122) is the continuum threshold above which 
the asymptotic freedom is restored and fp is the coupling between the electromagnetic 
current and the p field in vacuum; 

(0 I p) = /,m,e^. (3.123) 
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Assuming vector dominance in the medium we obtain 

9p = mp/f,. (3.124) 

In vacuum, the standard parameters for the p spectral function are given by, nip — 0.77 
GeV, = 0.14 GeV , fp = 0.141 GeV, gp = 5.46, coq = 1-3 GeV, S = 0.2 GeV and 
as = 0.3. 

Let us now concentrate on the spectral function in the cu channel. We again take a 
Breit-Wigner form along with a continuum: 

TmTT^ 1 ry 1 

In vacuum is the coupling of the current with the a; meson defined as 

{0\J^\u;)^ Um^e^. (3.126) 

Note that /^^ here is defined as factor 3 larger than Shuryak's definition [75] . Imll^, 
which is the imaginary part of the uj self-energy is given by, 

imn^(?o) = go(ra,-.3,r + r^7r-.,r^). (3.127) 

In vacuum the standard parameters for uj are as follows, = 0.782 GeV, ttItt = 0.14 
GeV, = 0.138 GeV, ujq = 1.1 GeV, 5 = 0.2 GeV and a, = 0.3. 

Since not much is known about the critical behavior of the scalar and tensor con- 
densates at finite T in QCD sum rules we take a simple ansatz for in-medium quantities 
for their T-dependence. A possible parametrization of *-quantities at finite T is 

(3.128, 

mv fv ^0 V ^c/ 
where A is a sort of dynamical critical exponent and V stands for vector mesons (p and 
a;). It may be noted that there is no definite reason to believe that all the in-medium 
dynamical quantities are dictated by a single exponent A. Since the numerical value 
of A is not known precisely, we take two typical cases: A = 1/6 (BR scaling) and 1/2 
(Nambu scaling) [77] with the following remarks: 

(i) Eq. (3.128) for m* is not entirely consistent with the low temperature theorem [184], 
which says there should be no O(T^) correction to the mass. Therefore, one cannot 
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Figure 3.16: Variation of vector meson mass with temperature for BR (long-dashed 
hne) , Nambu (dot-dashed hne) scahng with Tc— 200 MeV and in the QHD model for p 
(solid hne) and u; (dotted hne). 

take the ansatz too seriously at low T. For our purposes, however, T < 100 MeV is 
not relevant in any way since it is below the freeze-out temperature (130 MeV) we have 
considered. 

(ii) Local duality constraint h in QCD sum rules imphes that (/*)^ = 87r^(l-|-Q;s/7r)(a;Q)^ 
-I- (scattering (Landau damping) term) [180] which is shghtly violated for /* as defined 
in Eq. (3.128). 

(iii) The vector dominance assumption in the medium together with Eq. (3.128) simply 
leads to g* = gp. 

Under these reservations, we will use the parametrized spectral functions with the 
BR scaling and Nambu scaling ansatz in the calculation of the lepton and photon 
production in the following Section. The principal qualitative difference between the 
spectral function in QHD and that described in this Section is the existence of the 
continuum and its medium modification at finite T. 

In Fig. (3.16) we depict the variation of vector meson masses as a function of temper- 
ature in the BR and Nambu scaling scenarios. Results in the QHD model is also shown 
for the sake of comparison. The mass variation in the QHD model and BR scahng is 
slower than the Nambu scaling scenario. At higher temperature (near Tc) the QHD 
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Figure 3.17: Spectral function for the isovector (p) channel extracted from e'^e~ col- 
lisions (dotted hne) as a function of invariant mass. The dashed (solid) line indicates 
the spectral function when nip and cuo vary according to BR (Nambu) scaling. 

and the BR scaling results tend to converge. However, such a small difference in the 
mass variation in the above two scenarios may not be visible through photon spectra. 
We also note at this point that in QHD unlike the scaling scenarios the p and uj masses 
show different rate of reduction [185] due to different values of their coupling constants 
with the nucleons. In Fig. (3.17) the spectral function (Stt times Eq. (3.122)) for the 
isovector (p) channel is plotted as a function of invariant mass at T = 150 MeV and 
Tc — 160 MeV. We find that both the peak and the continuum threshold of the spectral 
function move towards lower invariant mass. In the case of Nambu scaling scenario 
the shift is more compared to BR scaling. In the Nambu scaling scenario the peak of 
the spectral function and the continuum are not well separated; a merging of the two 
would take place at T = Tc. This could possibly indicate the onset of a deconfinement 
phase transition. Fig. (3.18) shows the spectral function at T = 180 MeV and — 200 
MeV. Due to a larger separation between and T compared to the previous case the 
peaks in the spectral function in all the cases are well separated from the continuum. In 
Figs. (3.19) and (3.20) the spectral functions for the isoscalar (a;) channel obtained by 
multiplying Eq. (3.125) by 87r are shown for T=150 and 180 MeV respectively. In both 
cases the peaks in the spectral function corresponding to the BR and Nambu scalings 
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Figure 3.18: Same as Fig. (3.17) at T = 180 MeV and % = 200 MeV. 

is distinctly visible. The larger width in the isoscalar channel is due to the combined 
processes uj —>■ St: and a; tt — > tt tt as discussed before. 

The spectral functions for the vector mesons both in the isoscalar and isovector 
channels are plotted in Fig. (3.21) at a temperature T ~ Tg. As expected from the 
scaling law, the Breit-Wigncr peak has vanished due to its overlap with the continuum 
(see Eqs. (3.122) and (3.125)). All the hadrons in the thermal bath have melted to 
their fundamental constituents - the quarks and gluons. Such a spectral function would 
indicate a transition from hot hadronic matter to QGP. This behaviour should, in 
principle, be reflected in the dilepton spectrum originating from these channels. 

3.2 Photon Emission Rates from Hot Hadronic 
Matter 

In the energy regime of our interest the most important sources of photon production 
from hadronic matter are the reactions vrp 717, tttt p^, mi 777, vrr] — > 717, 
and the decays p —> 77777 ^ ~^ "^7 [55, 105, 185]. Apart from these we have also 
included those reactions which produce photons via the intermediary axial vector Oi. 
A non-zero width of vector and axial vector mesons in the intermediate state has been 
taken into account. 
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Figure 3.19: Spectral function for the isoscalar (a;) channel extracted from e+e~ col- 
lisions (dotted line) as a function of invariant mass. The dashed (solid) line indicates 
the spectral function when and loq vary according to BR (Nambu) scaling. 
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Figure 3.20: Same as Fig. (3.19) at T=180 MeV and Te= 200 MeV. 
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Figure 3.21: Spectral functions for isovector (p) and isoscalar (uj) channels at Tc . 

The relevant vertices for the reactions tttt ^ pj and tt p — > tt 7 and the decay 
p — > TT TT 7 are obtained from the following Lagrangian: 

>C = -gp^nfT ■ (7? X d^TT) - eJ^A, + ^F^- {p, X p;)3, (3.129) 

where F^^, — d^A,, — di^A^ is the field tensor for electromagnetic field and is the 
hadronic part of the electromagnetic current given by 

J"" = {p. X ^>^), + (7? X {d'^Tf + g^^Tf X p^))3, (3.130) 

with Qfj^^ = dfj_pi, - d^p^ - gpnnip,x X Pu)- The coupling strength of the pTrn vertex, gp^^, 
is fixed from the observed decay p tttt. 

Photon emission rates due to the reactions nr] — > 7r7, vrvr — > 777 and the decay 
a; — > 7r7 have been evaluated using the following interaction [158]: 

e^^^f^d^aj'-dy^TT + ^ e.^^^d^p'-dy^v + "^A.ff (3.131) 
nirj gp 

where e^j.^^ is the totally antisymmetric Levi-Civita tensor. The second term is con- 
structed analogously as the first term which is the familiar GSW Lagrangian. The 
last term is written down on the basis of Vector Meson Dominance (VMD) [162]. The 
values of gipp^ and g^jim are fixed from the observed decays, p ^ r]^ and a; — > 717 
respectively [185]. The constant gp is determined from the decay, p° — > e'^e" . 
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Figure 3.22: Thermal photon spectra at T = 160 MeV. Thick dots (squares) indi- 
cate photon emission rate from QGP including one and two loop contributions for 
Qs — 2(0.8). Dot-dash line represents photon spectrum from hot hadronic gas without 
medium effects. The result with the in-medium effects in the QHD model is shown by 
the long-dashed line. Dotted (solid) line indicates photon spectrum with BR (Nambu) 
scaling mass variation scenario. 

The importance of the role of oi as an intermediary meson in the process 7rp — > 
7r7 was emphasized in Refs. [186, 187]. Recently it has been shown [188] that this 
contribution is not so large. We use the following interaction Lagrangian for the vrpai 
and 7rai7 vertices [189, 190]: 



(2c + Z^)7r • X a'^ + ^^7? • {d„pu - d^p^) x {d'^t 



d'^TT ■ {di^p„ - d^pi^) X a" 



egpKef^ 
mi 



(3.132) 



where U/j, corresponds to the oi field, is the renormalization constant for pion fields 
and /7r(= 93 MeV) is the pion decay constant. The interaction terms with coefficients 
c and Kg are introduced to improve the phenomenology of the model. The following 
values of the parameters, — 1260 MeV, gp = 5.04, c = —0.12, = 0.17 and 
K,Q — 1.25 [190, 189] are considered which reproduce the width of the Oi meson in 
vacuum. 

The invariant amplitude for all the reactions mentioned above are listed in the 
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Figure 3.23: Thermal photon spectra at T = 180 MeV and Qs — 2. 

Appendix along with the corresponding Feynman diagrams. These are used in Eq. (2.98) 
to obtain the rates of photon production incorporating the effective masses and decay 
widths of the participating hadrons. While evaluating the photons from QGP we have 
considered both one loop and two loop contributions (as shown in Figs. (2.3) and (2.4)) 
to the photon self energy. 

The total photon emission rate from QGP and hadronic matter at T = 160 MeV 
is plotted in Fig. (3.22) as a function of the energy of the emitted photon for different 
values of strong charge Qs in the QGP phase and for various mass variation scenarios in 
the hadronic sector. The photon production rate from QGP has been evaluated in the 
HTL approximation, which is valid if the hard and soft scales are well separated, i.e. for 
Qs « 1 (which corresponds to ctg << 0.08). However, lattice QCD calculations [106] 
suggest that cts ~ 0.2 — 0.3 at the temperatures achievable in URHICs. We emphasize 
that the extrapolation of results obtained under HTL approximation to higher values of 
Qs (or as) may be dubious. We have evaluated the photon spectra for two values of the 
strong coupling constants = 0.8 (thick squares) and 2 (thick dots) to demonstrate 
the sensitivity of the photon spectra to the value of the strong charge and to show the 
uncertainties involved in the problem. In the hadronic sector the photon yield is seen 
to be enhanced in comparison to the case when the effects of the thermal interaction on 
the hadronic properties are neglected. This is true for almost the entire energy range of 
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the emitted photon under consideration. As a result of the similar mass shift in QHD 
and BR scaling the photon spectra in these two scenarios (long-dashed and dotted lines 
respectively) have a negligible difference, whereas the enhancement in the spectrum 
due to hadronic mass shift according to Nambu scahng is clearly visible (solid line) . In 
Fig. (3.23) wc show the photon emission rate at T = 180 McV. Photon spectra from 
hadronic matter with mass variation according to the Nambu scaling scenario overshine 
the photons from QGP even for a larger value of Qs (~ 2). 

The photon spectra at T = 150 MeV with in-medium masses calculated in the 
framework of gauged LSM, NLSM and HLS approaches is shown in Fig. (3.24). An 
increase of p mass in the NLSM reduces its number due to Boltzmann suppression 
which leads to a suppression in the photon emission rate (dotted line) . The production 
rate is enhanced due to a reduction in the p mass (solid line) in LSM. (We recollect that 
the p mass decreases in gauged LSM for low temperature and increases for temperatures 
close to the chiral transition temperature T^. Therefore, for T ~ wc will observe 
a reduction in photon emission rate and the net yield would be a superposition of all 
temperatures, from initial to freeze-out.) The change in the mass of p is so small in the 
HLS approach (short-dashed hue) that the production rate is almost indistinguishable 
from the spectra with vacuum masses of the hadrons. We have also demonstrated in 
Fig. (3.24) how the photon spectra is modified for a drastic change in the width of the 
p meson (Fp ~ 400 MeV) without any appreciable change in the pole mass {rup ~ 770 
MeV). Such a large value of the width have been proposed in Refs. [78, 191]. We observe 
that the effects of such modifications in the properties of p on the photon spectra is 
rather negligible (long-dashed line). 

So, from the above discussions we can infer that the relative photon yields from 
QGP and hot hadron gas depends on: (i) the value of the strong coupling constant, (ii) 
the degree of hotness of the medium and (iii) how adversely the hadrons are affected in 
the medium. 

We will now focus on a few other aspects that might affect the thermal photon 
emission from hot hadronic matter. 

The first is concerned with the fact that most of the photon producing reactions 
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Figure 3.24: The change in the photon spectra due to the finite temperature effects on 
the hadronic masses in Linear, Non-Linear Sigma Model and Hidden Local Symmetry 
approach at T= 150 MeV. 

which we have considered contain unstable particles (p and lu) in the external lines. 
Now, we know that the density of a stable hadron of mass m in a thermal bath is 
completely determined by the temperature, chemical potential and the statistics obeyed 
by the species through the distribution function 

dN M 1 



d?x d?k ds (27r)3 exp(A;o - n)/T±l 



5(s — m ) 



(3.133) 



where J\f is the degeneracy and ko — + s is the energy of the particle in the rest 
frame of the thermal bath. To account for "particles" like the p and uo mesons whose 
life-times are such that they can decay within the thermal system, the above expression 
is replaced by [185] 

dN _ M 1 
d^xd^kds ~ (2^ exp(A;o-/i)/T± 1 ^ 

where P{s) is the spectral function of the species under consideration and is given by 



P{s) = -- 



Imn 



(3.135) 



7r(s-m2 + Ren)2 + (Imn)2' 

[Note that P{s) reduces to the Dirac delta function in the hmit Imll — > 0.] The above 
exercise has been carried out wherever the p or uj appears in the photon producing 
reactions in the external lines. In the case these particles appear in the internal lines. 
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Figure 3.25: Effect of spectral function of vector mesons on photon emission rates at 
T — 160 MeV. Solid (dot-dashed) hne shows results with (without) the inclusion of 
spectral function. 



the propagators have been suitably modified, as discussed in the previous Section. The 
limits of the s integration in the case of different reaction channels have been determined 
from kinematical considerations. The results for the reactions tttt p7, tttt — > tt^, 
p ^ nir'-f and uj ^ iT'-f have been shown in Fig. (3.25). The difference caused by the 
inclusion of the spectral function is observed to be ~ a few percent. This is reflected in 
the total photon spectra as shown in Fig. (3.26). The effective mass and decay widths 
have been considered in QHD for the sake of illustration. 

The second point to consider is that though we have been treating the hadrons (p, 
u etc.) as point particles, in reality they are composite objects and may need form 
factors at high momentum transfer. We have investigated how the inclusion of form 
factors influence the results. Using the reaction tttt — > frf this has been demonstrated 
in Fig. (3.27). We have taken the same monopole form factor for both nnp and nn^ 
vertices [55] to suppress the contribution from very high momentum region where the 
quark structure of the hadrons could be relevant. The Ward-Takahashi identity has 
been used to obtain the dressed propagator. The in-medium mass of the p meson has 
been taken from QHD. The form factor effects for the above reaction reduces the photon 
production rate by about 10-15% at T —180 MeV. Once the space-time evolution is 
carried out such an effect will turn out to be neghgible. We have therefore neglected it 
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Figure 3.26: Total photon emission rate at T = 160 MeV. 



in our discussions. 
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Figure 3.27: The effect of the monopole form factor on the photon emission rate from 
the reaction tttt — > p^. 



Decay Photons 

Photons from -k^ and r] decays constitute almost the entire background to the spectra 
of thermal photons in the low transverse momentum region. We will briefly discuss how 
the transverse momentum spectra of decay photons can be estimated. One has to first 
obtain the spectra of the 7r° and r] mesons at the time of freeze-out. For the case of no 
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transverse flow the hadronic spectra is given by 
dN„ 



(3.136) 



n=l 



where, M is the degeneracy of the particles, t/ is the freeze-out time of the hadron 
and the plus and minus signs are for bosons and fermions respectively. Additionally, 
rriT = ^^Vt + m^, where m is mass of the hadron and pt is its transverse momentum. It 
is worthwhile to mention that the transverse expansion of the system may considerably 
enhance the transverse momentum of the the particles at the time of freeze-out. In any 
case, the decay photons have to be subtracted out by an invariant mass analysis, for 
the thermal photons to be identified. The decay photons from the process tt'^ — >■ 77, 
for example are then obtained as [192], 

= UL-k- y„f\ , (3.137) 

d?kT dy J E \ d^p J TT V 2 J ' ^ ' 

where p is the four-momentum of the 7r°, and k is the momentum of the detected photon. 
After some algebra, this reduces to, 

dN 1 



fOO 1-277 A]\J 

a'^KT dy TT k Jeq Jo a''p 



(3.138) 



where. 



Pt 



~k 



m 



^(E- Eo) sin^ 
k 



k\\ — kx sinhy , 
E 

Eq 



rriT coshy^r , 
2 

m 

k + 



k 



(3.139) 



3.3 Dilepton Emission Rates from Hot Hadronic 
Matter 



In order to discuss dilepton emission from hadronic matter we will use Eq. (2.84) which 
expresses the emission rate in terms of the retarded current correlation function. 

dR a 2m? f Am 



IvciW^^fBEiqo). 
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Putting d^q = 27r MdM MtcLMt dy this is written as 
dR ' 



dM 



a 1 



1 + 



2m^ 



'1 



4m^ 
M2" 



/ 



gMT cosh y/T _ ]^ W 



(3.140) 



The parametrized form of the electromagnetic current correlation function in the p and 
u} channels as discussed in Section 3.1.3c can be used in the above equation to obtain 
the dilepton emission rate from a thermal hadronic medium. However, instead of using 
the current correlation function directly in the above equation one can use vector meson 
dominance (VMD) to obtain the dilepton yield from tt'+'tt" — > which is known to be 
the most dominant source of dilepton production. In order to make a comparative study 
we state briefly how the dilepton emission rate from pion annihilation can be derived 
using Eq. (3.140). Recall that VMD relates the hadronic electromagnetic current to the 
vector meson field through the field current identity as 



(3.141) 



V 9v 

where, V = p,LJ,(f) and my stands for the bare masses. We shall keep only the p meson 
in the following. The electromagnetic current correlator can then be expressed in terms 
of the propagator of the vector particle in the following way: 



where 



1X1/ 



'■/J,!/ 



1 4 

^ImL)^-^ 



imn^^ 



(3.142) 



(g2 - m2 + Ren^^)2 + [Imn^^]2 



imn 



(3.143) 



(?2 - m2 + Renf )2 + [Imnf ]2 

The complete expression for the dilepton emission rate due to pion annihilation is then 
obtained from Eq. (3.140) as 

2m2\ \ 



dR 
dM 



2a'^mj 1 



1 + 



Aw? f MrdMrdy 



M2 



/ 



J^2 J gMT cosh y/T _ ]^ 



X 



+ 



imnf 



(M2 -ml + ReUfy + [ImUf]-^ {M^ -ml + ReUfy + [ImUf 



(3.144) 
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The quantities Imllj^ can be obtained from Eq. (3.59). It is however found that 
the difference between the longitudinal and transverse polarizations is very small upto 
reasonably high temperatures [101] for the interaction considered here. We hence make 
the approximation Uj^^ — 11^^ = 11''^, so that 



ImVl^^ = ^ 

yp 



(3.145) 



(g2 -ml + ReUpR)^ + [ImllpRY 
Writing Imll''^ = Mrf{M) and neglecting the thermal factor in the p width one obtains 

in the Boltzmann approximation. Ki is the modified Bessel function and a{M) is the 
cross-section for pion annihilation given by 



a{M) = Y^l - 4m2/M2 ^1 - 4m VM^ (1 + jM^) |F^(M)|^ (3.147) 



3M2 

with 



4 

ml 



which is known as the pion form factor. Apart from the thermal modification of the 
p meson properties, the rate given by Eq. (3.146) is equivalently derivable from the 
kinetic theory result given by Eq. (2.99). 

One should also add the contributions to the thermal dilepton yield coming from 
the decay of vector mesons. In a thermal medium the production of an off-shell vector 
meson (y) of four momentum q (where — M^) and its subsequent decay into a lepton 
pair leads to the dilepton emission rate [193] 

dR = T^p!. P'" fBE(qo)rr-.i+i- d% (3.149) 



(27r) 



where P/j^i, — J2poi ^ij. — ~9ij.v + Q.ij.(Iv/ is the projection operator as applicable for 
an unstable vector meson V, and 1™^^+;- is the partial decay width for the process 



V — ^ /+/ in vacuum, is the spectral function of the off-shell vector meson and is 



I Hi VclL.UUlil. p^ 

given by 



V 1 Imll^ 
^'^'^ " ^(g2-mf. + Ren«)2 + (Imn^)2 ^'^'^^ ^^"^^^^ 
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in the approximation Ily^ = 11™ = 11^. Using the relation P^'^ P^^ = (2 J + 1), we get 
the dilepton emission rate due to the decay of an unstable vector meson of spin J as 



dR ^ (2J+ 1) 



^vac 



imn^ 



(3.151) 



TT (g2 -ml + Ren^)2 + (Imn^)^ 
Note that for a particle which does not decay in the collision volume i.e. for which 
the total width Ftot = Imll-^/M is small, the spectral function in the above equation 
becomes 5{(f' — my), as it should be for a stable particle. As discussed before, in 
a medium the width Ftot of particle V should be calculated with all the processes 
involving the creation and annihilation of y, i.e. Ftot — Ty^aii — Taii^y [193]. 

Using similar approximations as in the case of pion annihilation, the invariant mass 
distribution of lepton pairs from the vector meson decays is obtained from Eq. (3.151) 
as. 



dM n 



2 



with 



(M2 -ml + Re°r4)2 + M^F^ ^^^'"^ ^" ' ^^'^^^^ 



Fr_z+.- = ^wJ^ - 4"^VM2 (1 + 2mVM2) (3.153) 



where m is the mass of the lepton. 

We will now discuss the results of static dilepton emission rates. The rates for BR 
and Nambu scahng scenarios are obtained by putting times Eqs. (3.122) and (3.125) 
in Eq. (3.120) which is then inserted in Eq. (3.140) using Eq. (3.121). The in-medium 
masses in this case is obtained from Eq. (3.128) with A=l/6(l/2) for BR(Nambu) 
scahng. The dilepton rate in the QHD model is obtained by adding the contributions 
from pion annihilation (Eq. (3.146)) and from u decay using Eq. (3.152). The effective 
masses of the p and ui in this case are as given in Eq. (3.54). It may be mentioned 
here that since we are considering dilepton production from pion annihilation with a 
p dominated form factor, adding the contribution from p decay to lepton pairs (which 
also contains the width due to p <-> tttt) separately will induce a double counting. 

In Fig. (3.28) we display the invariant mass distribution of e'^e~ pairs. The dilepton 
yield from qq annihilation is denoted by thick dots. The dotted line indicates the result 
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Figure 3.28: Thermal dilepton spectra at T = 150 MeV. Thick dots indicate dilepton 
emission rate from QGP. Thin dotted hne represents dilepton yield from hot hadronic 
gas without medium effects. The result with the in-medium effects in the QHD model 
is shown by the dot-dashed line. Long-dashed (solid) line indicates dilepton spectrum 
with BR (Nambu) scaling mass variation scenario. 

obtained from the parametrization of the electromagnetic current correlation function 
in the 'p' and 'a;' channels, when the medium effects are ignored. A large shift towards 
the lower invariant mass region of the p peak is seen in the Nambu scaling (solid line) as 
compared to the BR scaling (dashed line) consistent with the relative shift in the spectral 
functions in the two cases as discussed before. In the QHD model calculations (dot-dash 
line) the two peaks corresponding to p and uj masses are visible in the spectra [148, 157]. 
The separation between the two peaks is due to different mass shift of the p and ui. 
Measurement of such separation in hadronic masses (Am = m* — m*) would signal 
the in-medium effects. The validity of such results could be tested in URHICs by 
the CERES [194, 195] collaboration in future. Similar shift at zero temperature but 
finite baryon density could be detected by HADES [196] and CEBAF [64]. Effects of 
the continuum contribution in the spectral function on the dilepton spectra is clearly 
visible for M > 1 GeV since the value of the continuum threshold in vacuum is 1.3 
GeV. Due to the continuum contribution the dilepton rates from hadronic matter and 
QGP shine equally brightly in the mass range M > 1 GeV. The lepton pair spectra at 
T = 180 MeV is shown in Fig. (3.29). Since the effective mass of the p in QHD and BR 
scahng scenario is almost same in this case (see Fig. (3.16)), the corresponding rates 
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Figure 3.29: Thermal dilepton spectra at T = 180 MeV. 
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Figure 3.30: Thermal dilepton spectra at T = T^.. The vector meson peaks have merged 
with the continuum and the rates from QGP and hadronic matter have almost become 
equal. 



are very similar near the p peak. 

The dilepton invariant mass distribution at T = Tc is shown in Fig. (3.30). All 
the peaks in the spectrum have disappeared as expected. The rates obtained from the 
electromagnetic current correlator is close to the rate from qq annihilation, indicating 
that the qq interaction in the vector channel has become very weak i.e. signalling the 
onset of deconfinement [197,198]. 

In Fig. (3.31) we compare the dilepton emission rate at T — 150 MeV for vacuum 
mass of p with the rates in which the effective masses are obtained in the framework 
of gauged LSM, NLSM and HLS approach. The positive shift of p mass in NLSM is 
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Figure 3.31: The change in the invariant mass distribution of lepton pairs due to the 
finite temperature effects on the hadronic masses in Gauged Linear, Non-Linear Sigma 
Model and Hidden Local Symmetry approach at T= 150 MeV. The long-dashed fine 
indicates dilepton spectra for Tp — 400 MeV and rUp — 770 MeV (see text) . 

reflected in the peak position of the spectra towards larger value of M (dotted line). 
A very small change in the p mass in HLS approach does not cause any visible change 
in the invariant mass distribution of dilcptons. The long-dashed line indicates lepton 
pair distribution for nip — 770 MeV and Tp = 400 MeV. The large width of p leads to 
the disappearance of the p peak from the spectra, indicating that p ceases to exist as a 
quasi-particle. However, as mentioned before, it is interesting to note that the photon 
spectra is insensitive to such drastic broadening of the p meson. 



Chapter 4 



Electromagnetic Spectra with 
Space-Time Evolution 



Matter produced in highly relativistic coUisions of heavy ions can be either in the form 
of a hot hadronic gas or a quark gluon plasma. So far we have discussed the rate of 
photon and dilepton emission per unit time from unit volume from a thermal system 
made up of quark matter or hadronic matter at a fixed temperature. However, the 
highly excited state of matter produced at a high temperature will expand and cool 
emitting photons and dileptons in the process. This process will continue as long as 
the mean free paths of the constituents become comparable to the size of the system 
and freeze-out occurs. So, in order to compare the yield with experiments we must 
integrate the static rates over the space time volume of the collision from formation till 
freeze-out. The total yield is hence 



where dV stands for invariant phase space elements: d^q/E for photons and d^g for 
dileptons. T{x) is the local temperature, assuming that the produced matter is in local 
thermal equilibrium, x is the space-time coordinate and d'^x — dVdt where V is the 
three- volume. Since we are dealing with relativistic situations one must account for the 
fact that the thermal rates are evaluated in the rest frame of the emitting matter and 
hence the momenta of the emitted photons or dileptons are expressed in that frame. 



freeze — out 



dN 




(4.1) 



formation 
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Therefore, the invariant rate is a function of E*, the energy of the photon or lepton 
pair in the rest frame of the emitting matter which in this case is the fluid element. In 
a fixed frame like the laboratory or the centre of mass frame, where the 4-momentum 
of the photon or lepton pair is — {E , q) and the emitting matter element d^x moves 
with a velocity -u^ = 7(1, -u), wc have E* = u^q^. As is evident from Eq. (4.1) the key 
ingredients that we now need in order to carry through our program are firstly, the 
initial conditions viz. the initial time and the initial temperature, and secondly, the 
temperature as a function of space-time or, in other words the equation of state (EOS). 
These are handled by the relativistic hydrodynamic equations which are basically the 
equations of conservation of energy-momentum and particle current. We have seen in 
Chapter 3 that the hadrons redress themselves in the thermal medium thereby changing 
their vacuum masses. This feature is taken into account in the evolution dynamics 
through the EOS. As a result the coohng law in the hadronic sector will be quite 
different from the QGP which is made up of weakly interacting quarks and gluons 
which are essentially massless. We will discuss these aspects in Sees. 4.1 and 4.2. We 
then go on to evaluate the photon and dilepton yield in typical collision scenarios at 
the CERN SPS and BNL RHIC energies in Section 4.3. Our principal contention in the 
whole exercise is to comment on the initial probable state of matter, QGP or hadronic 
matter on the basis of our calculations in comparison with available data. We will 
compare the results of our evaluation of the thermal photon spectra with the data from 
the WA80 and WA98 experiments. The dilepton data will be compared with the results 
of the CERES experiment. 

4.1 Relativistic Hydrodynamics 

A basic ingredient of the hydrodynamic description of the collision volume is the exis- 
tence of a strong interaction time scale, 

n ~ -r-^^ ~ Ifm/C ~ Tformation (4.2) 

^^QCD 

In any hadronic collision the produced fragments can only interact after a proper time 
Ti has elapsed after their collisions. Thus, there is another time scale in the problem. 



116 




Region Z=0 

Central Rapidity Region 



Projectile 
Fragmentation 



Region 



(a) 







dN/dy / 











' ► 



Target 



(b) 



y 

Proj. 



Figure 4.1: Schematic diagram of (a) the coUision geometry in the central coUisions of 
identical nuclei and (b) the central plateau in the rapidity density. 



the so called transit time , which is defined as 



"^transit 



Tcm 



(4.3) 



Ra is the nuclear radius, is the Lorentz factor. If the value of 7cm (which is a function 
of the collision energy) is such that Ttransit < ^formation then most of the secondaries are 
formed after the nuclei pass through each other. Consequently these secondaries will 
not contribute to the energy density of the fluid in the central region (see Fig. (4.1a)). 
Such a scenario may be realized in the highly relativistic collisions of heavy ions. This 
particular feature has been taken into account in Bjorken's hydrodynamic model [29]. 

It has been observed experimentally that the particle spectra for the secondaries 
produced va. N — N collisions exhibit a central plateau in the rapidity space. This kind 
of behaviour is due to the frame independence symmetry of the hydrodynamic expansion 
of the system [28]. Bjorken assumed that the same kind of plateau will also be observed 
in nucleus- nucleus collisions [29] . In terms of the initial condition this means that the 
energy density, pressure etc. (all the thermodynamic quantities) will be a function of 
the initial thermalization (proper) time only and will not depend on the space time 
rapidity r] (deflned later). This initial symmetry of the thermodynamic quantities is 
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preserved throughout the evolution scenario. If the particle rapidity density is flat or 
invariant under Lorentz boosts then the entropy density (s) will be independent of the 
rapidity. 

The evolution of the fluid is governed by the energy-momentum conservation equa- 
tion 

df,T^"' = (4.4) 

where 

Ti"' = (e + P)m'^m" + g'^'P (4.5) 

is the energy-momentum tensor for ideal fluid; e and P are the energy density and 
pressure of the fluid element. For an isentropic flow the entropy conservation reads 

d^s^'^O (4.6) 

where s'^ — su'^ is the entropy current. Considering fluid motion along the z (beam) 
direction, u^^ — (coshy/, 0, 0, sinhy/) where the fluid rapidity y/ is given by 

yf = -\n^^, (4.7) 

Vz being the fluid velocity. Let us consider the situation in a two dimensional sub-space 
{t — z plane). We will make a change of variables from {t, z) to (r, 77) where 

r^VW^; ^ = ^ln^, (4.8) 

T] being the space-time rapidity. As a result Eqs. (4.4) and (4.6) become 

d d 

— {sT cosh(y/ - ^)) + ^ smh(y/ - ry)) = (4.9) 
d d 

— {Trsmhivf + ^ (rcosh(y^ - 77)) = 0. (4.10) 

In order to solve these equations we need to know how the fluid velocity Vz depends 
on the space-time coordinates, i.e. we need yf as a function of r]. We will make use 
of the observation that the amount of particle production is intimately connected with 
the stopping of the colliding nuclei. At high collision energies the nuclei would slow 
down considerably, but due to the large rapidity gap between the colliding nuclei, the 
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produced matter will still span a large rapidity interval as seen in Fig. (4.1b). The 
rapidity distribution of secondaries is likely to exhibit a plateau-like structure in the 
central region as conjectured by Bjorken. In this situation it is reasonable to assume 
that in the central rapidity region the longitudinal velocity of produced matter exhibits 

scaling behaviour, 

V, = z/t. (4.11) 

We then have yf = 7] and r becomes the proper time of the fluid frame which is related 
to the centre of mass frame by a Lorcntz transformation along the 2;-axis with velocity 
z/t. Putting Uf — rj in Eqs. (4.9) and (4.10) we get 

^ (sr) = (4.12) 



and 



dT 

-^0. (4.13) 



Equation (4.13) implies that T is independent of r] and so are all the thermodynamic 
quantities; they depend only on r. This gives rise to the frame independence or boost- 
invariant expansion. Prom Eq. (4.12) we have 

ST — const. (4-14) 

which is the Bjorken's scahng solution. The resulting space-time picture of the collision 
is shown in Pig. (1.2) of Chapter 1. It may be noted that the above results were 
obtained without any specific input from the equation of state. It is thus a general 
result that one dimensional similarity fiow is necessarily isentropic even if there is a 
phase transition. Por a relativistic massless gas with statistical degeneracy gk, s and T 
are related through the equation of state: 

' = "^"^SkT". (4.15) 

Putting this expression for entropy density in the Bjorken scaling solution we get the 
cooling law 

TV = const. (4.16) 
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which is routinely used to evaluate the signals of QGP. The initial temperature of the 
system is determined by observing that the variation of temperature from its initial 
value Tj to final value Tf (freeze-out temperature) with proper time (r) is governed by 
the entropy conservation (Eq. (4.12)) 

s{T)t = s{Ti)n . (4.17) 

The entropy density is then expressed in terms of the observed particle (pion) multi- 
phcity. Using Eqs. (4.15) and (4.17) one gets the initial temperature as 



' " 45C(3)7ri?i4afcri dy ^ ' ' 



where dNT^/dy is the total pion multiplicity, Ra is the radius of the system, Tj is the 
initial thermalization time, and C(3) is the Reimann zeta function, — (7r^/90) Q}^ is 
the degeneracy of the produced system and k stands for either QGP or a hot hadronic 
gas. The rapidity density for the secondaries is obtained as [199], 

(f)..-"'(f)„ 

where is known as the rcscattering parameter. dN/dy can be parametrized to 
fit the experimental data in the central region as 

=o.81n(^). (4.20) 

The assumption of a central plateau in the rapidity distribution is not experimentally 
observed in nucleus-nucleus collisions at the presently available energies. Hence the 
boost invariant hydrodynamics may not be a valid concept at these energies. The con- 
cept of complete stopping in Landau model [200] is not valid either at these energies. 
The physical situation may be in-between the boost invariant model of Bjorken and 
the Landau model of complete stopping, which means that there may be an overlap 
between the formation zone and the collision zone. However, for its simplicities the 
Bjorken model will be used to describe the space time evolution of matter formed in 
URHICs. Appropriate generalization has been made to take into account the temper- 
ature dependent hadronic masses in the evolution. 
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It is worth mentioning that though significant stopping of baryons have been ob- 
served in relativistic collisions of heavy ions in the recently concluded experiments at 
the CERN SPS we will assume the central region in rapidity space to be approximately 
free of baryons. With strong stopping the net baryon density in the central region 
can be appreciable but the matter is still not baryon rich because strong stopping also 
causes large secondary particle production and baryons form only a small fraction of all 
particles. In other words, the ratio of the baryon chemical potential to the temperature 
is small. However, at RHIC and LHC the net baryon number can certainly be neglected 
in the central rapidity region. 



4.2 Initial Conditions and Equation of State 

The set of hydrodynamic equations is not closed by itself; the number of unknown vari- 
ables exceeds the number of equations by one. One thus needs to postulate a functional 
relation between any two variables so that the system becomes deterministic. The most 
natural course is to look for such a relation between the pressure P and the energy 
density e. Under the assumption of local thermal equilibrium, this functional relation 
between P and e is the EOS. Obviously, different EOS's will govern the hydrodynamic 
flow quite differently [201] and as far as the search for QGP is concerned, the goal is to 
look for distinctions in the observables due to the different EOS's (corresponding to the 
novel state of QGP vis-a-vis that for the usual hadronic matter). It is thus imperative 
to understand in what respects the two EOS's differ and how they affect the evolution 
in space and time. Recently, the sensitivity of the photon emission rate on various 
evolution scenarios has been studied in Ref. [202]. 

A physically intuitive way of understanding the role of the EOS in governing the 
hydrodynamic flow lies in the fact that the velocity of sound = {dP/de)s sets an 
intrinsic scale in the hydrodynamic evolution. One can thus write a simple parametric 
form for the EOS: P — c1{T)€. For an ideal gas of massless constituents = 1/3. 
Inclusion of interactions, however, may drastically alter the value of [203]. In our 
calculation we assume the MIT bag model equation of state for the QGP where the 
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energy density and pressure are given by 

^Q = 9Q^ + B, (4.21) 

and 

Pq-9q^T'-B. (4.22) 

The effective degrees of freedom in QGP, qq — 37 for two flavours. The entropy 
density sq is given by sq = 2c/Q(7r^/45)T^. Putting Uk = ag ^ (7r^/90)g'Q the initial 
temperature for a system produced as QGP can be determined from Eq. (4.18). 

In the hadronic phase we have to be more careful about the presence of heavier 
particles and the change in their masses due to finite temperature effects. The ideal limit 
of treating the hot hadronic matter as a gas of pions originated from the expectation 
that in the framework of local thermalization the system would be dominated by the 
lowest mass hadrons while the higher mass resonances would be Boltzmann suppressed. 
Indirect justification of this assumption comes from the experimental observation in high 
energy collisions that most of the secondaries are pions. Nevertheless, the temperature 
of the system is higher than during a major part of the evolution and at these 
temperatures the suppression of the higher mass resonances may not be complete. It 
may therefore be more reahstic to include higher mass resonances in the hadronic sector, 
their relative abundances being governed by the condition of (assumed) thermodynamic 
equilibrium. The hadronic phase is taken to consist of tt, p, cu, r), Oi mesons and nucleons. 
The nucleons and heavier mesons may play an important role in the EOS in a scenario 
where mass of the hadrons decreases with temperature. 

The energy density and pressure for such a system of mesons and nucleons are given 

by 

Yl [ d'pEh fBE{Eh, T) + / d?pEM fFD{EN, T) (4.23) 

/i=mesons V^^^ l^^^ 

and 

Y. / iS^h.(E,, T) + ^^l d^J^JME., T) (4.24) 

/i=mesons ^ ' "" ^ ' 
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where the sum is over aU the mesons under consideration and N stands for nucleons 
and Efi — ^Jp'^ + m\. The entropy density is then 

SH = = 4aeff(r) = A'^g^^{wr{nT)T' (4.25) 

where g^e is the effective statistical degeneracy. Thus, we can visuahze the finite mass of 
the hadrons having an effective degeneracy ges{m*{T),T). Because of the temperature 
dependence of the effective degeneracy Eq. (4.18) has to be solved self-consistent ly in 

order to calculate the initial temperature of the system initially produced as a hot 
hadronic gas. We thus solve the equation 

^-^nR\4a,,m)TS (4.26) 

where aeff (Tj) = (vr^/QO) gef[{m*{Ti),Ti) . The change in the expansion dynamics as well 
as the value of the initial temperature due to medium effects enters the calculation of 
the photon emission rate through the effective statistical degeneracy. 

If the energy /entropy density in the fireball immediately after the so-called "forma- 
tion time" Tj is sufficiently high, then the matter exists in the form of a QGP. As the 
hydrodynamic expansion starts, the system begins to cool until the critical temperature 
Tc is reached at a time tq. At this instant, the phase transition to the hadronic matter 
starts. Assuming that the phase transition is a first order one, the released latent heat 
maintains the temperature of the system at the critical temperature Tc, even though 
the system continues to expand; the cooling due to expansion is compensated by the 
latent heat liberated during the process. Together with the possible explosive events, 
we are neglecting the scenarios of supercooling or superheating. This process continues 
until all the matter has converted to the hadronic phase at a time th, still at T = Tc; 
from then on, the system continues to expand, governed by the EOS of the hot hadronic 
matter till the freeze-out temperature Tf is reached at the proper time rj. Thus the 
appearance of the so called mixed phase at T = Tc, when QGP and hadronic matter 
co-exist, is a direct consequence of the first order phase transition. The possibility of the 
mixed phase affects the bulk features of the evolution process and plays an important 
role in QGP diagnostics. 
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In the mixed phase, the relative proportion of QGP and hadronic matter must be a 
function of time; initially the system consists entirely of QGP and at the end, entirely of 
hot hadronic matter. If we denote the fraction of the QGP by /q(t) then the hadronic 
fraction in the mixed phase is /i/(T) — 1 — /q(t) so that fqiTq) — 1 and fHiTn) — 1- 
Then the entropy density in the mixed phase is given by 

SmUr) ^ fqir) s'^Q + Mr) s% (4.27) 

where Sq {sjj) denotes the entropy density of QGP (hadronic) phase at Tc. The life-time 
of the mixed phase is given by 

rmr' =rH- tq. (4.28) 

The scaling law governing the variation of s with r must continue to hold also in the 
mixed phase; substituting Eq. (4.27) in Eq. (4.12) we obtain for Tj > 

/«« = ^(^?-)-,^(t-) (-^) 

where r (= gg/ges) is the ratio of the degeneracy of QGP phase and the effective 
degeneracy in the hadronic phase. In the above equation we have used the relation 
th — ttq, valid for (1 + 1) dimensional isentropic expansion. 

Since the entropy density s(t) in any phase can be expressed as in Eq. (4.27) with 
suitable values of /q(t), the volume fractions can also be expressed as 



fQ{r) = 



S[T) - Sh 



M-) = (1.30) 

If Tj = Tc, i.e. if the system is formed in the mixed phase with a fraction Jq of the 
QGP phase then 

(4.31) 



fQir) ^ 



:i + (r-l)/o)^-l 



r — 1 

The mixed phase in this case ends at a proper time = (1 + (r — l)/o)Tj. For Sj < s^, 
the value of /ij(T) is always unity. 

Now that we have all the necessary requisites, we will evaluate the photon and 
dilepton yield in heavy ion coUisions at the SPS and RHIC. We will consider two 
possibilities: 
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A+A^(QGP)^(Mixed Phase) ^Hadronic Phase 

or 

A + A — >Hadronic Phase. 

The former (latter) case where the initial state is formed in QGP (hadronic) phase 
will be called the 'QGP scenario' ('no phase transition scenario'). For the QGP sec- 
tor we use a simple bag model equation of state (EOS) with two flavour degrees of 
freedom. The temperature in the QGP phase evolves according to Bjorkcn scaling law 
T^T — TfTi- The cooling law in the hadronic sector is quite different from that of 
the QGP because of the presence of massive hadrons. These hadrons redress them- 
selves in the medium thereby changing their vacuum masses. This is accounted for by 
introducing temperature dependence in the statistical degeneracy which takes care of 
the mass varying with temperature. In Table 1 we quote the values of the initial tem- 
peratures obtained by assuming various mass variation scenarios. The value of initial 
thermalization time has been assumed as 1 fm/c both for SPS and RHIC energies. The 
multiphcity density {dN / dy) in the two cases are given as 600 and 1735 respectively. 
For a two-flavour QGP the effective degeneracy is 37. The freezeout time Tj is taken 
as 130 MeV [143] in all the calculations, d and 5 dictate the variation of temperature 
with proper time for the hadronic matter according to the cooling law T — -d/r^ . The 
values of 6 indicates a slower coohng in the hadronic phase as compared to that of QGP 
phase where T oc l/r°-^^. 
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=4.6 fm 




(MeV) 
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5 


Th (fm) 


^? 


S 


bare masses 


270 


10.8 


0.267 


0.215 


31.9 


0.337 


0.215 


QHD 


220 


9.4 


0.247 


0.194 


27.6 


0.305 


0.194 


BR 


195 


8.2 


0.236 


0.184 


23.9 


0.288 


0.185 


Nambu 


195 


4.7 


0.203 


0.151 


13.9 


0.239 


0.152 



Table 1 : Values of initial temperatures and evolution parameters for SPS and RHIC. 
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Figure 4.2: Variation of effective degeneracy as a function of temperature. 

In Fig. (4.2) we depict the variation of effective degeneracy as a function of temper- 
ature with and without medium effects on the hadronic masses for various scenarios. 
We observe that for T > 140 MeV the effective degeneracy becomes larger due to 
the reduction in temperature dependent masses compared to the free hadronic masses. 
Physically this means that the number of hadrons in a thermal bath at a temperature 
T is more when in-medium mass reduction is taken into account. Eq. (4.26) implies 
that for a given pion multiplicity the initial temperature of the system will be lower 
(higher) when medium effects on hadronic masses are considered (ignored). This is 
clearly demonstrated in Fig. (4.3) where we show the variation of temperature with 
proper time for different initial conditions. The thick dots indicate the yield when 
QGP is formed initially at Tj = 185 MeV and cools down according to Bjorken law up 
to a temperature Tc at proper time tq, at which a phase transition takes place; it re- 
mains constant at up to a time th — 9.4 fm/c after which the temperature decreases 
as T = 0.247/r°'^^^ (when medium effects arc taken from QHD) to a temperature Tj. 
If the system is considered to be formed in the hadronic phase then the initial temper- 
ature is obtained as Tj = 220 MeV (270 MeV) when in-medium effects on the hadronic 
masses from the QHD model is taken into account (ignored) . The corresponding cooling 
laws are displayed in Fig. (4.3). The above parametrizations of the cooling law in the 
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Figure 4.3: Variation of temperature as a function of proper time in Pb-Pb collisions 
at the SPS. 

hadronic phase have been obtained by solving Eq. (4.17) self-consistcntly. An initial 
state with the vanishing meson masses at Tj — 195 MeV (tj = 1 fm/c) could be realized 
in the case of BR and Nambu scahng scenarios with pion multiplicity dN/dy — 600. 

At RHIC a scenario of a pure hot hadronic system within the format of the model 
used here, appears to be unrealistic. The initial temperature considering bare hadronic 
masses turns out to be ~ 340 MeV whereas for the other extreme case of massless 
hadrons it is ~ 290 MeV. With temperature dependent masses the initial temperature 
will lie somewhere between these two values. For such high temperatures, clearly a 
hot dense hadronic system cannot be a reality, the hadrons would have melted away 
even for lower temperatures. Thus, for RHIC we have treated the case of a QGP 
initial state only. The temperature profile for RHIC is depicted in Fig. (4.4) where 
we observe that the length of the plateau, which indicates the life time of the mixed 
phase r^{^ = th — tq, depends on the masses of the hadrons in the hadronic phase. 
The effective degeneracy plays an important role here. At the transition point there is 
a large decrease in the entropy density. This decrease has to be compensated by the 
expansion (increasing the volume) to keep the total entropy constant. Since we are 
considering (1+1) dimensional isentropic expansion, this change in the entropy density 
will be compensated by increasing r which is a measure of the volume, so that the total 
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Figure 4.4: Variation of temperature as a function of proper time in Pb-Pb collisions at 
RHIC. The initial temperature has been determined by assuming 'QGP scenario'. The 
initial temperature Ti — 265 MeV for = 1 fm/c and dN/dy — 1735 

entropy remains constant. We have seen earher (Fig. 4.2) that the effective degeneracy 
in the hadronic phase is the largest for the Nambu scahng and smallest for the bare 
mass scenario, resulting in smallest (largest) discontinuity in the entropy density for 
the former (latter) case. Consequently the time taken for the system to compensate the 
decrease of the entropy density in the Nambu scaling scenario is smaller compared to 
the case where bare masses are considered. Hence the life time of the mixed phase for 
the Nambu scaling case is smaller than all other cases. 



Our ultimate goal in URHICs is to distinguish between the two possibilities described 
above, that of a hot hadronic gas initial state or a QGP which undergoes a phase tran- 
sition to hadronic matter via a mixed phase of coexisting quark and hadronic matter. 
In the following we will compare the photon and dilepton spectra originating from these 
two scenarios. As discussed before, the experimentally observed photon and dilepton 
spectra originating from an expanding QGP or hadronic matter is obtained by convo- 
luting the static (fixed temperature) rate with expansion dynamics. In a first order 
phase transition scenario the photon and dilepton spectra from a (1 + 1) dimensionally 
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expanding system is obtained as 



/ 



( 



dR 



QGP 



Q{s - 4) 




) 



HG "^Q 



HG 



QGP 



Q{Sff — s) Tdrdr) 




— s 



H 



9(4 - s)e{s - sj,) 



(4.32) 



where Ra is the radius of the nuclei and © functions are introduced to get the con- 
tribution from the QGP, mixed and hadronic gas (HG) phases. For (1+1) dimen- 
sional expansion, d^x — TrR\ dzdt — TrR\ Tdrdr). In the case of photons dV stands for 
d'^p/E{= d'^pTdy); Pt and y being the transverse momentum and rapidity of the emitted 
photon. For dileptons, dV = d^q = MdMd'^qrdy where, M, qt and y are the invariant 
mass, the transverse momentum and rapidity of the lepton pair respectively. 



The thermal photon spectra at SPS energies is shown in Fig. (4.5). In order to interpret 
the spectra corresponding to the different scenarios let us recall that the effective degen- 
eracy in the hadronic phase Qes is obtained as a function of T by solving Eq.(4.25). A 
smaller (larger) value of g^s is obtained in the free (effective) mass scenario. As a result 
wc get a larger (smaller) initial temperature by solving Eq.(4.26) in the free (dropping) 
mass scenario for a given multiplicity. Naively we expect that at a given temperature 
if a meson mass drops its Boltzmann factor will be enhanced and more of those mesons 
will be produced leading to more photons [105, 204]. However, a larger drop in the 
hadronic masses results in smaller initial temperature, implying that the space time 
integrated spectra crucially depends on these two competitive factors. Therefore, with 
(without) medium effects one integrates an enhanced (depleted) static rate over smaller 
(larger) temperature range for a fixed freeze-out temperature (T/ = 130 MeV). In the 
present calculation the enhancement in the photon emission due to the higher initial 
temperature in the free mass scenario (where static rate is smaller) overwhelms the 
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Figure 4.5: Total thermal photon yield corresponding to dN/dy — 600 and Ti — 1 fm/c. 
The solid (long-dash) line indicates photon spectra when hadronic matter formed in the 
initial state at Tj = 195 McV (Tj = 220 MeV) and the medium effects are taken from 
Nambu scaling (QHD). The dotted line represents the photon spectra without medium 
effects with Tj = 270 MeV. The solid (dotted) line with thick dots represent the yield for 
the 'QGP scenario' when the hadronic mass variations are taken from Nambu scahng 
(free mass). 

enhancement of the rate due to negative shift in the vector meson masses (where the 
initial temperature is smaller). Accordingly, in the case of bare mass ( Nambu scaling) 
scenario the photon yield is the highest (lowest). In case of the QHD model, the photon 
yield hes between the above two hmits. In the 'QGP scenario' the photon yield with 
in-medium mass is lower than the case where bare masses of hadrons are considered. 
However, the difference is considerably less than the 'no phase transition scenario'. This 
is because, in this case the initial temperature is determined by the quark and gluon 
degrees of freedom and the only difference between the bare and effective mass scenarios 
is due to the different lifetimes of the mixed phase. In Fig. (4.5), the photon spectra in 
the 'no phase transition scenario' overshines the ones from the 'QGP scenario'. 

We now compare the above set of curves to the preliminary single photon data 
obtained by the WA98 Collaboration [205] in Fig. (4.6). We find that for the entire 
range of pt the photon spectra in the 'no phase transition scenario' with medium effects 
evaluated in the QHD as well as Nambu scaling scenario explain the data reasonably 
well. However, the yield with the QHD model parameters seems to overpredict in the 
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Figure 4.6: Same as Fig. (4.5) with data from WA98 experiment. 

low region by a small amount. A QGP initial state as well as a hadronic gas initial 
state with vacuum properties of hadrons seem to be improbable. The thermal photon 
spectra for Pb-Pb coUisions at SPS have also been evaluated in Ref. [202] for a variety 
of evolution scenarios without introducing medium effects. 

We have also compared the thermal photon spectra evaluated [142] in the QGP 
and hadronic gas scenarios with the data obtained by the WA80 [206] collaboration. 
This is shown in Fig. (4.7). The experimental data in this case is the upper bound 
of the thermal photons in S-Au collisions at 200 A GeV at the CERN SPS. The pion 
multiplicity is given as dN/dy —22b. Our conclusions are very similar to the Pb-Pb 
case. The photons evaluated in the hadronic gas scenario with bare masses (short- 
dashed line) can be ruled out. Also, in such a scenario when medium effects from QHD 
are incorporated, the yield (solid line) just crosses the upper bound. However, thermal 
photons in the 'no phase transition scenario' with medium effects according to Nambu 
scaling (dot-dashed line) turns out to be a viable description. Quite a few attempts have 
been made to explain the WA80 upper bound using various models. Using relativistic 
hadron transport model Li et al [207] have shown that the photon spectra with either 
free or in- medium meson masses in dense matter do not exceed the upper limit. With 
in-medium modifications of mesons calculated in the hidden local symmetry approach, 
Halasz et al [188] have reached similar conclusions. However, using (3+1) dimensional 
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Figure 4.7: Total thermal photon yield in S-Au collisions at 200 A GeV at CERN SPS. 
The long-dashed line shows the results for a QGP initial state with Tj=190 MeV. The 
solid (dot-dash) line indicates photon spectra in the hadronic gas scenario at Tj = 220 
MeV (Ti = 195 MeV) with medium effects calculated from QHD (Nambu scahng). The 
short-dashed line represents the photon spectra from a hadronic gas with bare masses 
and Tj = 270 MeV. The initial time in all these cases have been taken as 1.2 fm/c. 
The data correspond to the upper bound obtained by WA80 Collaboration. 

expansion as well as various probable equation of states SoUfrank et al [208] does not 
claim evidence for a phase transition. The WA80 data have also been explained by Hui 
et al [209] in a hadron and string cascade model. These results are in sharp contrast 
to the calculations performed by Srivastava and Sinha [210], Shuryak et al [211, 212], 
Dumitru et al [213], Arbcx et al [214], Neumann et al [215] and Song [204] who have 
concluded that the WA80 upper limits are satisfied only if a QGP is assumed to have 
been formed in the initial stages. 

4.3.2 Dilepton Spectra at SPS 

The space time integrated dilepton spectra for the 'QGP scenario' and 'no phase transi- 
tion scenario' with different models of mass variation are shown in Fig. (4.8). The shifts 
in the invariant mass distribution of the spectra due to the reduction in the hadronic 
masses, particularly the p meson mass, according to the different models are distinctly 
visible. Their is no appreciable difference between the QGP and 'no phase transition' 
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Figure 4.8: Total thermal dilepton yield corresponding to dN/dy ~ 600 and Ti — 1 fm/c. 
The solid (long-dash) line indicates dilepton spectra when hadronic matter formed in 
the initial state at Tj = 195 McV {Ti = 220 MeV) and the medium effects are taken 
from Nambu scaling (QHD). The dotted line represents the yield from hadronic gas 
scenario with bare masses and Tj = 270 MeV. The thick dots represent the yield for 
the 'QGP scenario' with mass modifications included. 



scenarios which is expected because one is looking at low values of the invariant mass 
where the dominant contribution in either case will be from hadronic matter. Dileptons 
originating from the high temperature QGP phase will show up at higher values of M. 

Let us now make a comparison of our predictions with the dielectron spectrum in 
Pb-Au collisions at 158 A GeV as measured by the CERES Experiment [216]. In this 
case one has to incorporate [208] the kinematical cuts in the momenta of the electron 
and the positron relevant for the CERES detector. In Fig. (4.9) we have plotted the 
dilepton yield normalized to dNch/dy within the rapidity interval 2.1 to 2.65 in the 'no 
phase transition scenario' with the spectral functions in the QHD and Nambu scaling 
scenarios. We observe that dileptons from an evolving hot hadronic gas with medium 
modifications calculated in the Nambu scaling scenario clearly fits the data beyond 
M ~ 500 MeV. The QHD model calculations under similar conditions overshoots the 
data in the region of the shifted p peak between 500 and 700 MeV. 
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Figure 4.9: Our calculations compared with the inclusive e^e spectra measured by 
CERES in 158 A GeV Pb-Au collisions in the '96 run. 



Though not a direct signal of QGP, the observed enhancement of low-mass dileptons 
compared to the yield from hadronic decays at freeze-out (shown by the dotted curve 
in Fig. (4.9)) has triggered a host of theoretical activity. Though the pion annihilation 
channel tt+tt" via the p meson accounts for a large fraction of this enhancement, 

it turns out that a quantitative explanation of the data requires the incorporation of 
medium modifications of the vector mesons. Li, Ko and Brown [143] were the first 
to achieve an excellent agreement with the CERES data using a decreased p mass in 
the dense fireball. Rapp et al [78] have used a large broadening of the p meson line 
shape and hence a shorter life time due to scattering off baryons in order to explain the 
data. However, both approaches rely on a high baryon density for the dropping mass 
or the enlarged width of the p meson but the role of baryons is still a debatable issue. 
It is worth mentioning that we have assumed zero baryon density in our calculations. 
Koch [217] also finds very little effect due to baryons and are able to obtain a reasonable 
explanation of the data. It is worth emphasizing here that as yet it has not been possible 
to explain [208, 218] the observed low-mass enhancements of dileptons measured in the 
Pb-Au coUisions as well as in the S-Au collisions at the CERN SPS in a scenario which 
does not incorporate medium induced dropping of the vector meson mass, the p meson 
mass in particular [212]. 
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Figure 4.10: Thermal photon spectra at RHIC energies. A first order phase transition 
scenario has been considered. 

4.3.3 Photon and Dilepton Spectra at RHIC 

Finally we study the electromagnetic probes for RHIC energies. As discussed earher, 
in this case a scenario of a pure hot hadronic system appears to be unrealistic. So we 
have treated the case of a QGP initial state only. The baryonic chemical potential has 
been taken to be zero, which in this case is certainly a valid assumption as far as the 
central rapidity region is concerned. The thermal photon yield for RHIC is displayed 
in Fig. (4.10). The sohd hne represents the total thermal photon yield originating from 
initial QGP state, mixed phase and the pure hadronic phase. The short dash line 
indicates photons from quark matter (QM) (= pure QGP phase + QGP part of the 
mixed phase) and the long dash line represents photons from hadronic matter (HM) (= 
hadronic part of the mixed phase + pure hadronic phase) . In all these cases the effective 
masses of the hadrons have been taken from Nambu scaling. For pt > 2 GeV photons 
from QM overshines those from HM since most of these high photons originate from 
the high temperature QGP phase. The dotted and the dot-dash lines indicate photon 
yields from QM and HM respectively with bare masses in the hadronic sector. The 
HM contribution for the bare mass is larger than the effective mass (Nambu) scenario 
because of the larger value of the hfe time of the mixed phase in the earher case (see 
Table 1). It is important to note that for pt > 2 GeV, the difference in the QM and 
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Figure 4.11: Thermal dilepton spectra at RHIC energies in a phase transition scenario. 

HM contribution in the effective mass scenario is more that the bare mass scenario. 

Thermal dilepton yield at RHIC energies for QGP initial state for different mass 
variation scenarios are shown in Fig. (4.11). The shape of the peak in the dilepton 
spectra in case of QHD is slightly different (broader) from the other cases because of 
the larger mass separation between p and u; mesons in this case (see Fig. (3.16)). The 
dilepton yield beyond the vector meson peak is larger in the bare mass scenario because 
of larger initial temperature as seen in Table 1. 




Chapter 5 

Summary and Discussions 



As emphasized in Chapter 1, electromagnetic signals are the only direct probes of evolv- 
ing strongly interacting systems and hence the study of their spectra is of tremendous 
importance as far as the detection of QGP in URHICs is concerned. In this thesis 

we have studied photon and dilcpton production from a thermal source composed of 
strongly interacting matter, QGP or a hot hadronic gas, likely to be produced in the 
ultra-relativistic collisions of heavy ions. 

In Chapter 2 we have reviewed the formulation of photon and lepton pair production 
from a thermal medium in the framework of Thermal Field Theory. We have seen that 
the production rate per unit time per unit volume is proportional to the electromagnetic 
current correlation function or, more generally, to the spectral function of the photon 
in the medium. Neglecting the possible reinteractions of the virtual photon on its way 
out of the medium, the emission rate is found to be proportional to the retarded self 
energy of the photon. In Chapter 2 we have also reviewed how Hard Thermal Loop 
(HTL) resummation technique has been used to cancel the infra-red divergences so as 
to obtain finite rates of photon emission due to quark and gluon interactions in the 
QGP. 

As mentioned in the introduction, our main emphasis has been to study the effect of 
the medium modifications of hadronic properties, particularly of the p and uj mesons, on 
the spectra of photons and dileptons obtained in URHICs. In Chapter 3 we have made 
a detailed study of the static (fixed temperature) rates of photon and dilepton emission 
from hadronic matter incorporating medium modifications of the masses and decay 
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widths of vector mesons evaluated within various well-known models. In the Quantum 
Hadrodynamic (QHD) model the p and uj masses are found to decrease differently. 
The disentanglement of the p and uj peaks in the dilepton spectrum resulting from 
URHICs would be an excellent evidence of in-medium mass shift of vector mesons [194] 
and/or the validity of such model calculations for the situation under consideration. 
We have studied the modification of the spectral functions of the p and ui mesons using 
QCD sum rules at finite temperature. Here, the hadronic spectral function in vacuum 
have been parametrized using the experimental data on e^e~ — > hadrons in terms of 
a resonance and a continuum. Due to the lack of our understanding of the critical 
behaviour of scalar and tensor condensates the vector meson masses and continuum 
threshold as a function of temperature were taken to vary according to Brown-Rho 
and Nambu scaling. The interesting possibility of the disappearance of the p and cu 
peaks in the dilepton spectra due to temperature induced lowering of the continuum 
threshold was investigated. We have also made a brief study of the p mass modifications 
within the gauged linear and non-linear sigma models and the hidden local symmetry 
approach. The changes in these models are observed to be small. The static rates of 
photon and dilepton production from interacting hadronic matter have been evaluated. 
It is observed that the in-medium effects on the dilepton and the photon spectra are 
prominently visible for QHD model and Brown-Rho and Nambu scaling scenarios. It is 
interesting to note that the dilepton spectra spectra is affected by both the changes in 
the decay width as well as in the mass of the vector mesons whereas the photon spectra 
is affected only by the change in the mass of the vector mesons and is rather insensitive 
to the change in width. The effects of the continuum on the dilepton spectra are seen 
to be substantial. 

In Chapter 4 the photon and dilepton yield from a longitudinally expanding system 
likely to be produced in URHICs at the SPS was contrasted between a QCP initial state 
and a hot hadronic gas initial state scenarios. It was observed that the incorporation 
of medium modifications resulted in substantial changes in the evaluation of the initial 
temperature as well as the cooling profile. When compared with the preliminary Pb- 
Pb data obtained by the WA98 Collaboration we find that the photon yield in the 
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hot hadronic gas scenario with medium effects calcuiated from QHD as well as Nambu 
scaling scenarios explain the data reasonably well. Similar conclusions were reached 
when we compared the thermal photon yield in S-Au coUisions at 200 AGeV to the 
upper limits obtained by the WA80 Collaboration. The thermal photon spectra from a 
hadronic gas with spectral modifications evaluated in the Nambu scaling scenario falls 
within the experimental upper limits. The corresponding QHD model calculations does 
not seem to qualify as a viable possibility in this case. 

We have also compared the results of the dielectron yield to the Pb-Au data obtained 
by the CERES Collaboration. Here again the hot hadronic gas scenario with medium 
effects in the Nambu scaling scenario seem to fit the data fairly well. This is certainly 
due to the increased production of low mass electron pairs originating mainly from the 
decay of p mesons with decreased mass and enhanced width in the medium. This is 
however not a direct signal of QGP because the contribution from the QGP phase is 
very small in the low invariant mass region. 

It is interesting to note that so far the observed low mass enhancement in the 
dilepton data could not be explained without the consideration of p mesons with a 
depleted mass. In this connection it is worth mentioning that CERES is planning to 
substantially improve the mass resolution to values of the order of the natural line width 
of the u! meson (8.5 MeV). It should then be possible to directly measure the yield of the 
vector mesons p, uj and even (r<^=4.5 MeV) including any possible changes in their 
mass or decay width [195]. Also, the mass shifts of vector mesons at zero temperature 
but finite baryon density could be detected by HADES [196] and CEBAF [64]. 

In the following we will discuss some general aspects regarding photon and dilepton 
emission from evolving matter with reference to the calculations presented in this thesis: 

The photon production from QGP has been evaluated using HTL resummation 
based on the assumption « 1, which is impossible to meet in URHICs even at 
the highest energy to be available at the CERN LHC in future. The strong coupling 
constant is hkely to attain a value g's ~ 2 at RHIC/LHC. Evaluation of the photon 
emission rates at such high values of the strong coupling is a formidable task. In this 
respect the development of methods suitable for addressing non-perturbative effects 
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near and above the QCD phase transition point is of paramount importance. Extension 
of the self-consistent resummation scheme developed in 0^ theory [219] to non-abelian 
gauge theory [220, 221] would be a very important step towards the understanding of 
the phenomena near the QCD phase transition. 

The exact value of the critical temperature (T,.) for deconfinement phase transition is 
still uncertain. However, recent lattice simulation [13] for two flavour QCD indicates a 
value of Tc for chiral transition ~ 130-160 MeV. We have taken = 160 MeV, although 
it is not known whether the values of for the chiral and deconfinement transition are 
the same or not. The value of the initial thermalization time is unfortunately also 
an unknown quantity. We take = 1 fm/c as a canonical value following Bjorken [29]. 
A similar value of Tj has been considered in the hterature [51, 204, 208]. The freeze-out 
temperature T/ is another parameter which is not precisely known. This is important 
because the yield from the hadronic phase depends substantially on the value of Tf. We 
have also considered a baryon-frcc system based on the assumption that the ratio of the 
baryon chemical potential to the temperature is reasonably small. This might not be 
totally justified at SPS energies, but at RHIC a baryon-free central region is certainly 
a possibility. 

We have assumed thermodynamic equilibrium in our calculations. Such a situation 
may not be realized practically [31, 32, 222, 223, 224, 225]. Also, emissions from the 
pre-equilibrium era may pollute the kinematical domains where one looks for signals of 
QGP [32]. Unfortunately, although considerable progress has been made [226, 227, 228], 
the general techniques for solving non-equilibrium quantum field theoretical problems 
is still in the early stages of development [229]. 

We have assumed that the produced system behaves as an ideal, non-dissipative 
fluid. Incorporation of viscosity breaks the time-reversality of the evolution. The ensu- 
ing generation of entropy during the temporal evolution invalidates the role of dN/dy 
as a handy constant of motion and consequently affects the estimation of the initial 
temperature [230]. 

In this thesis we have neglected transverse expansion of the evolving matter. In- 
corporation of medium modifications in a 3-1-1 dimensionally expanding system is a 
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non-trivial task. We intend to take up this project in the near future. 

We conclude with a few comments based on the experimental results from URHICs 
at the CERN SPS. Prom the measured hadron yields and momentum spectra including 
transverse and directed flow eflects there is a strong evidence of production of a col- 
lectively behaving strongly interacting matter with a large volume and finite lifetime. 
The thermal photon spectra as well as the pattern of low-mass enhancement of e+e" 
pairs observed in S-Au, Pb-Pb and Pb-Au collisions also lend support to the picture 
of a collective thermal system of hadrons. On the other hand, the observation of a 
significant enhancement in multistrange baryon production as well as the anomalous 
suppression of J/^ in Pb-Pb collisions cannot be explained without invoking partonic 
degrees of freedom and deconfined matter. At this critical and interesting juncture we 
look forward to the experiments at RHIC where larger and hotter systems are likely 
to be produced. Along with quantitative gains in all the signals there is a distinct 
possibility of direct detection of electromagnetic radiation from the QGP phase which 
will pave the way for an unambiguous conclusion regarding the formation of this novel 
form of matter. 
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Invariant Amplitudes for Photon Production 

We list the invariant amplitudes for photon production from hadronic matter consisting 
of TT, p, uj, 7] and ai mesons which have been considered in the evaluation of the photon 
yield from hot hadronic matter. The medium modifications of the p and cu mesons have 
been taken into account in the respective propagators. 
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Figure A.l: Feynman diagrams for tt'^tt"" — > p^^ 
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Figure A.2: Feynman diagrams for n^n^ — > 7 
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Figure A. 3: Feynman diagrams for n^p^ 
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Figure A. 4: Feynman diagrams for ir^p^ — > tt*^ 7 
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Figure A. 6: Feynman diagrams for tttt — > 777 and nr] — > 777. 
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and r(p ^ 777) = (57±10.5) keV. 
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Figure A. 7: Feynman diagram for tt p — > Oi — > tt 7. 
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